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the  Grafenberg  array  (6  to  13  stations)  were  averaged  for  each  event  and  comparer 
to  the  single  station  (Al)  results*  Additionally,  transverse  component  data 
were  available  for  three  of  the  GRF  stations.  At  MAIO  only  vertical  component 
Lg  data  were  available.  For  the  propagation  path  from  Semipalatinsk  to  GRF,  the 
narrow-band  filtering  method  and  the  noise-corrected  power  spectra  of  both 
vertical  and  transverse  components  give  an  average  Q  for  Lg  of  approximately  1000 
For  the  path  from  Azgir  to  GRF  the  average  Q  is  around  1100.  For  the  path  from 
Semipalatinsk  to  MAIO  the  average  Q  is  approximately  625.  There  appear  to  be  no 
significant  differences  in  the  results  among  individual  events  for  each  path, 
indicating  that  the  strong  tectonic  release  effects  in  evidence  at  long  periods 
for  some  of  the  Semipalatinsk  events  do  not  influence  the  Lg  spectra  of  those 
events  appreciably  in  the  frequency  band  .3  to  1.2  Hz,  and  the  best  band  to  use 
to  establish  an  vs  nt  relationship  is  centered  at  .6  to  .7  Hz. 

To  overcome  long-period  non-isotropic  source  effects  for  Shagan  River  event! 
several  very  low  tectonic  release  events  are  used  to  obtain  an  empirical  "pure" 
explosion"  Rayleigh  signature  (shape)  which  then  is  incrementally  scaled  and  sub¬ 
tracted  from  the  total  Rayleigh  wave  signature  of  high  tectonic  release  events 
until  the  best  double-couple  solution  that  also  fits  the  Love  wave  signals  is 
found  (i.e.  when  the  trace  of  the  moment  tensor  is  closest  to  zero);  the  scaled 
signal  is  then  the  appropriate  explosion  contribution.  The  remaining  "pure 
earthquake"  Rayleigh  and  Love  wave  signals  for  each  event  are  interpreted  in 
terms  of. a  best-fitting  equivalent  earthquake  model. 

To  determine  the  approximate  threshold  for  obtaining  source  parameters  from 
long-period  surface  waves  in  an  intraplate  setting,  the  phase  matched  filtering 
method  is  used  to  analyze  the  1982  New  Brunswick  aftershock  sequence  with  the 
>  ^elC-recprded  mainshock  seismograms  as  the  reference  signals.  The  source  excita¬ 
tion  function  of  each  aftershock  is  extracted  from  the  observed  signal  spectrum 
•  by  correcting  for  the  path  transfer  function  and  the  main  shock  source  excitata- 
t'ion.  Preliminary  results  using  this  approach  on  recordings  at  GDSN  stations  in 
North  America  indicate  that  source  parameters  (including  depth)  can  be  extracted 
for  events  as  ^small  as  about  m.  4  to  4.5.  (This  will  be  part  of  a  M.S.  thesis  by 
B.  Yan).  D 

The  1982  New  Brunswick  earthquake  series,  including  the  m^  5.7  mainshock  and 
several  high  magnitude  aftershocks  were  investigated  to  determine  their  source 
mechanisms  and  focal  depths.  Band-pass  filters  applied  over  a  broad  frequency 
range  provided  group  velocity  and  spectral  amplitude  data  for  Rayleigh  and  Love 
modes  for  the  path  to  SCP.  Synthetic  amplitude  spectra  were  generated  assuming 
structure  parameters  common  to  both  the  source  area  and  the  propagation  path. 

Key  aspects  of  the  source  mechanism  could  be  inferred  from  the  fundamental  mode 
observations,  particularly  for  the  mainshock  on  January  9,  1982,  and  the  largest 
magnitude  aftershock  which  occurred  on  January  11.  Although  there  are  signifi¬ 
cant  differences  in  modal  excitation  at  the  higher  frequencies  from  event  to 
event,  only  qualitative  support  for  source  mechanism  could  be  inferred  because 
higher  order  modes  could  not  be  identified  at  the  higher  frequencies.  Addition¬ 
ally,  published  NEIS  teleseismic  nr  ,  network  m,  ^  ,  and  Mg  values  were  compared 
with  m^L  and  M  values  computed  with  the  SCP  aafa.  Thele  results  indicated  that 
the  sing?e-statfon  m,  values  for  the  four  New  Brunswick  events  studied  in 
detail  agreed  remarkably  well  with  the  network  average  magnitudes,  whereas  Mg 
values  commonly  disagreed.  In  spite  of  the  observed  variable  excitation  of  the 
individual  higher  order  modes,  consistent  with  the  theoretical  excitation  appro¬ 
priate  for  the  different  mechanisms  for  these  events,  the  ensemble  of  modes  that 
comprises  the  Lg  signal  gives  every  indication  of  providing  a  stable  measure  of 
source  strength  when  nr  is  measured  according  to  Nuttli's  magnitude  relationship. 

The  fundamental  moat  Love-to-Rayleigh  surface  wave  spectral  ratios  for  periods 
greater  than  ten  seconds  provided  estimates  of  source  mechanism  and  focal  depth — 
for  the  mainshock  and  its  larger  aftershocks.  The  accuracy  of  focal  depth 
determinations  is  approximately  ±  1  km. _ _ 
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SUMMARY 

Grant  No.  AFOSR-82-0054 


The  research  completed  on  this  project  Is  summarized  as  follows: 

1.  A  statistical  study  was  made  to  estimate  source,  path,  and  receiver 
terms  from  two  Lg  data  sets  for  North  America:  (a)  LRSM  measurements  for 
nearly  100  Nevada  Test  Site  (NTS),  and  (b)  WWSSN  and  Canadian  recordings  of  17 
eastern  North  American  earthquakes.  Decay  rates  of  r~2  are  acceptable  for 
both  data  sets,  with  a  tendency  to  a  lesser  rate  over  paths  in  the  Canadian 
Shield,  where  a  mean  crustal  Q  of  approximately  2000  is  indicated.  For  the 
NTS  events  the  source  terms  correlate  very  well  with  yield.  No  dependence  of 
Lg  excitation  on  non-isotropic  source  effects  was  found;  this  includes 
transverse  component  Lg.  This  analysis  method  has  been  extended  to  obtain 
similar  decomposition  of  source,  path  and  station  effects  as  a  function  of 
frequency  for  Lg  signatures  of  NTS  and  Eurasian  events. 

2.  Previous  studies  have  shown  that  excitation  of  fundamental  mode  Love 
and  Rayleigh  waves  corresponding  to  earthquake-like  source  contributions  is 
highly  variable  among  nearby  explosions  at  NTS  and  at  Semipalatinsk.  The 
associated  focal  mechanisms  typically  are  consistent  with  regional  tectonic 
stress  orientations.  As  discussed  above  we  have  found  there  is  not  a 
corresponding  excitation  of  Lg  waves  for  NTS  explosions,  and  consequently 
magnitude-yield  estimates  using  Lg  will  not  be  biased  significantly  by  these 
non-isotropic  source  effects.  For  Soviet  explosions  in  thrust-type  tectonic 
settings,  we  have  found  that  above  about  0.3  Hz  the  vertical-  and 
transverse-component  Lg  signals  also  are  not  significantly  biased,  but  unlike 
NTS,  excitation  of  somewhat  lower  frequency  vertical-component  Lg  in  some 
cases  appears  to  be  correlated  with  non-isotropic  source  contributions 
observed  in  the  long-period  data  from  the  Shagan  River  test  sites.  Lg  data 
above  about  0.3  Hz  seem  to  provide  a  reliable  additional  means  of  estimating 
yield  for  events  at  both  test  sites  that  is  comparable  to  P-wave  yield 
estimates  when  propagation  and  station  effects  are  properly  taken  into 
account. 

3.  Several  approaches,  applied  separately  or  in  combination,  have  been 
used  in  an  attempt  to  overcome  the  biasing  effects  of  non-isotropic  source 
excitation  on  seismic  yield  estimates  from  long-period  observations.  They 
fall  into  two  general  categories:  (a)  decomposition  of  observed  seismic 
signatures,  especially  long-period  Love  and  Rayleigh  wave  data,  using  a  moment 
tensor  (or  equivalent)  inversion  formalism  to  estimate  the  explosion  signal, 
and  (b)  use  of  carefully  selected  seismic  phases,  frequency  bands,  azimuths, 
and  distances  that  are  least  affected  by  non-isotropic  effects.  The  success 
in  either  case  depends  on  the  nature  and  strength  of  the  non-isotropic 
contributions.  Different  strategies  are  required  for  strike-slip  tectonic 
settings  (NTS)  and  thrust-type  tectonic  settings  (Semipalatinsk).  An 
empirical  correction  procedure  was  developed  that  uses  a  suite  of 
low-tectonic-release  events  to  model  Che  explosion  contribution  in  the 
high-tectonic  release  long-period  signals. 

4.  We  have  found  that  for  similar  epicentral  distances  Lg  excitation 
(relative  to  m^)  is  comparable  among  granite  explosions  at  Semipalatinsk  and 
Hoggar  in  north  Africa.  Preliminary  work  was  directed  to  the  important 


question  of  whether  NTS  granite  explosions  have  comparable  or  larger  Lg 
excitation  for  a  given  m^  (yield)  than  Hoggar  or  Oegelen  explosions. 

5.  Several  experiments  were  carried  out  to  separate  near-receiver 
effects  from  source  effects  on  body  waves,  to  use  mode-converted  phases  to 
estimate  source  and  receiver  structure,  and  to  put  bounds  on  Q  at  source  and 
receiver  locations.  In  one  of  these  studies  (a  M.S.  project)  a 
frequency-domain  phase-difference  polarization  filter  has  been  implemented  and 
used  as  an  aid  in  discrimination  of  wave  types  present  in  three-component 
seismic  signals.  Crustal  structure  beneath  five  three-component  seismic 
stations  in  eastern  North  America  has  been  investigated  using  P-to-SV 
conversions  identified  for  teleseismic  earthquake  and  underground  nuclear 
explosion  sources.  These  results  are  contained  in  an  M.S.  thesis  (D.  O'Neill) 
completed  in  late  1983. 

6.  Methods  to  increase  the  dynamic  range  of  DWWSSN  long-period  and 
short-period  digital  recordings  were  developed  and  implemented  using  station 
SCP  as  an  example.  A  paper  describing  this  approach  has  been  prepared  for 
publication  and  is  included  in  the  Semi-Annual  Technical  Report  for 

March  1,  1983  -  August  31,  1983. 

7.  Fundamental  and  higher  mode  surface  waves  from  the  January  9,  1982 
New  Brunswick  m^,  5.7  mainshock  and  several  aftershocks  are  being  used  to 
investigate  source  mechanisms  and  focal  depth  for  these  events.  They  provide 
an  opportunity  for  comparing  source  parameters  obtained  from  regional 
observations  with  those  from  teleseismic  observations  of  the  same  events;  this 
source  region  is  also  a  good  analog  of  Semipalatinsk  and  the  earthquakes  are 
in  a  granite  pluton.  A  M.S.  paper  that  was  completed  on  this  topic  (C. 
Nichols)  appearts  at  the  end  of  this  report. 

8.  Spatial  wavenumber  analysis  of  multi-spectral  scanner  Landsat  images 
has  been  used  to  study  the  orientation  and  location  of  structural  lineaments 
in  the  region  of  the  Hoggar  massif  in  southern  Algeria,  Africa.  All  results 
are  compared  to  a  previous  photointerpretation  of  the  imagery  and  to  available 
seismic  observations  on  the  mechanism  of  "tectonic  release"  associated  with 
nuclear  explosions  detonated  in  the  massif.  Major  results  of  the  wavenumber 
analysis,  photo-interpretation  and  seismic  observations  are  in  agreement  with 
one  another.  In  particular,  the  N38°E  strike-slip  fault  inferred  by  Harkrider 
(1976)  from  analysis  of  surface  wave  radiation  patterns  is  consistent  with  the 
prominent  northeast-striking  features  observed  in  both  the  wavenumber  spectra 
and  the  original  imagery.  From  this  study  and  previous  work  relevant  to  the 
area,  the  maximum  principal  compressive  stress  axis  is  inferred  to  have  been 
in  the  east-west  direction  in  the  past,  but  at  present  is  oriented  in  the 
north-south  direction.  This  orientation  is  consistent  with  the  present-day 
direction  of  Af rican-Eurasian  plate  motion  and  with  focal  mechanisms  of  recent 
larger  earthquakes  to  the  northwest  near  El  Assam,  Algeria.  A  revised  crustal 
structure  model  for  this  source  area  was  also  obtained.  A  M.S.  thesis  (C.  Ng) 
describing  these  results  has  been  completed  and  a  paper  for  publication  has 
been  prepared.  These  results  are  Included  in  the  Semi-Annual  Technical  Report 
for  March  1,  1983  -  August  31,  1983  and  March  1,  1984  -  August  31,  1984. 

9.  Seismic  experiments  were  carried  out  to  study  the  vertical  and 
laterialheterogeneities  in  the  Earth's  interior  that  can  affect  body-wave 
arrivals  and  their  codas  (or  precursors  for  S  waves).  Various  signal  analysis 


techniques  were  employed  to  separate  near-receiver  effects  from  near-source 
effects,  and  near-receiver,  mode-convereted  (P-SV,  SV-P)  arrivals  were  interpreted 
in  terms  of  receiver  crust  and  upper  mantle  structure.  Significant  variations  in 
continental  crust  and  upper  mantle  structure  were  found  upon  comparing  results  for 
stations  RKON,  WH2YK,  and  LASA.  High  frequency  S  arrivals  from  P-SV  conversions 
from  the  400  km  discontinuity  at  RKON  from  Soviet  explosions  confirm  that  Q  is 
very  high  in  the  upper  mantle  beneath  that  part  of  the  Canadian  Shield.  A  study 
of  PcP-P  travel  time  residuals  show  that  significant  lateral  heterogeneity  exists 
near  the  core-mantle  boundary.  There  is  some  evidence  that  these  deep-seated 
heterogeneities  correlate  with  surface  tectonic  features  such  as  continental 
shields. 

10.  For  each  of  a  number  of  Soviet  explosions  two  different  spectral 
approaches  were  used  to  determine  the  average  Q  for  Lg  signals  arriving  in  a 
velocity  window  ranging  approximately  from  3.5  to  3.0  km/sec.  The  mean  spectral 
slope  based  on  all  the  events  analyzed  is  then  used  to  calculate  the  average 
attenuation  coefficient  (QU)~*.  Results  from  individual  stations  of  the 
Grafenberg  array  (6  to  13  stations)  were  averaged  for  each  event  and  compared 
to  the  single  station  (Al)  results.  Additionally,  transverse  component  data 
were  available  for  three  of  the  GRF  stations.  At  MAIO  only  vertical  component 
Lg  data  were  available.  For  the  propagation  path  from  Semipalatinsk  to  GRF, 
the  narrow-band  filtering  method  and  the  noise-corrected  power  spectra  of  both 
vertical  and  transverse  components  give  an  average  Q  for  Lg  of  approximately 
1000.  For  the  path  from  Azgir  to  GRF  the  average  Q  is  around  1100.  For  the 
path  from  Semipalatinsk  to  MAIO  the  average  Q  is  approximately  625.  There  appear 
to  be  no  significant  differences  in  the  results  among  individual  events  for  each 
path,  indicating  that  the  strong  tectonic  release  effects  in  evidence  at  long 
periods  for  some  of  the  Semipalatinsk  events  do  not  Influence  the  Lg  spectra  of 
those  events  appreciably  in  the  frequency  band  .3  to  1.2  Hz,  and  the  best  band  to 
use  to  establish  an  ®bL,g  vs  “b  relationship  is  centered  at  .6  to  .7  Hz. 

To  overcome  long-period  non-isotropic  source  effects  for  Shagan  River  events 
several  very  low  tectonic  release  events  are  used  to  obtain  an  empirical  "pure 
explosion”  Rayleigh  signature  (shape)  which  then  is  incrementally  scaled  and 
subtracted  from  the  total  Rayleigh  wave  signature  of  high  tectonic  release 
events  until  the  best  double?-couple  solution  that  also  fits  the  Love  wave 
signals  is  found  (i.e.  when  the  trace  of  the  moment  tensor  is  closest  to  zero); 
the  scaled  signal  is  then  the  appropriate  explosion  contribution.  The 
remaining  "pure  earthquake"  Rayleigh  and  Love  wave  signals  for  each  event  are 
interpreted  in  terms  of  a  best-fitting  equivalent  earthquake  model. 

To  determine  th<  approximate  threshold  for  obtaining  source  parameters  from 
long-period  surface  waves  in  an  intraplate  setting,  the  phase  matched  filtering 
method  is  used  to  analyze  the  1982  New  Brunswick  aftershock  sequence  with  the 
well-recorded  mainshock  seismograms  as  the  reference  signals.  The  source 
excitation  function  of  each  aftershock  is  extracted  from  the  observed  signal 
spectrum  by  correcting  for  the  path  transfer  function  and  the  main  shock  source 
excitation.  Preliminary  results  using  this  approach  on  recordings  at  GDSN 
stations  in  North  America  indicate  that  source  parameters  (including  depth)  can 
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If  tho  wavatraia  of  lnttrut  cannot  bo  completely 
taalaco4  la  cha  seismic  troco.  Spectral  leakage  will 
obscure  features  la  tba  amplitude  rolloff  obovo  tho 
coraor  frequency.  Moreover.  a  straightforward  FTT . 
ooao  after  tha  application  of  a  prudent  tapar  (a. g. . 
lUnaiac)  la  subject  to  largo  perlodogram  varlaaco  and 
alll  alaroproaaat  tha  frequency  content  of  a  pulse 
tnoan  to  bo  greatly  dlepereed.  We  will  daacrlbo 
miltiple-tapar  rlaa  aarloa  tochnlquos,  dovalopod  in 
tho  recant  eaglaearing  literature,  that  over cone 
aany  of  theae  41f f lcultloa.  We  have  teated  theae 
techniquaa  on  data  collected  (ran  the  Ante  Selamic 
TdMtind  Network  of  ten  digitally- recording  at  at  lone 
oanpled  at  250  tinea  per  aecond.  with  thla  high  quali¬ 
ty  data,  wa  have  enaainod  regional  wavetraina  froa 
earthquakes  and  one  NTS  blaat.  Several  dlatlnct 
•pactral  feature*  are  found,  both  above  end  below  the 
inferred  corner  frequency,  that  correlate  between 
horltontal  and  vertical  comments.  We  will  conpare 
tha  aUtltapar  apectra  with  ain«le  window  opactra. 
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TeleaelMlc  SV  wave  a  have  been  generally  Ignored  In 
wave  propagation  and  aource  etudlea  becauae  of  mown 
eoopilcationa  la  wave  propagation  for  etructure  near 
the  aource  eod  near  the  receiver.  The  validity  of 
conaoe  optic  ray  and  MRJ  aalanograa  aethoda  (or 
counting  SV  aynthetlc  aelanograaa  la  enenined  by 
confuting  aynthetlc  aelanogrena  uaing  theae  techniquaa 
and  coapariog  thee  to  SV  aynthetlca  produced  t roe  a 
vavenueoer  integration  technique,  both  ray  aathoda 
give  a  poor  epproataatlon  to  the  wave  propagation  for 
dlatancna  lean  than  60*.  Diffracted  bp  and  the  SPL 
wave  interfere  with  near  aource  phaaea.  auch  ea  S,  pS, 
and  aS  for  a  ahallow  aetaalc  aource,  producing 
anonaloualv  high  MpUiMai  and  coaflea  weveforw  in 
agreenent  with  oboervatlonal  enperience.  becauae  of 
the  Noho  Sp  phaae  and  diffracted  bp  phaae.  the 
*»ftlc»l  conpooeat  of  notion  above  greater  dlatortion, 
relative  to  tlw  ray  theory  raault,  than  doee  tne 
radial  couponeat  of  notion,  kay  theory  appeare  to  be 
appropriate  for  the  Initial  2u  aeconda  of  the  SV 
vavatrain  fron  a  ahallow  aource  for  rangea  greater 
than  bO*.  SV  wavea  iron  deep  aourcea  are  leaa 
affected  by  diffracted  bp  and  SPL  than  SV  iron  ahallow 
aourcea. 
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0'  IX  .101,1.0  0,  smllo. 

•■*rf>oua«fs  ‘  -*■  — 


JOHN  BOATWRIGHT  (USGS,  Hpnlo  Park,  CA  94025) 
5E0N6E  CHOY  (USCS,  Gold*!.  CO  90225) 


Thq  energy  f|g«  contains  in  thq  P-wave  orouos  ( P-pP 
•IP)  or  thf  S-wave  qrouos  (S*oS*sS)  raotatao  Oy  a  mal¬ 
low  cartnouaae  is  est  mated  from  me  assort  ion  tnat 
the  energy  f\u*  in  the  direct  ano  oeotn  onases  ados  in¬ 
coherently.  B,  denning  enera'.weo  radiation  patterns 

*ncn  incorporate  this  neutral  interference,  tne  wave 
J^ojos  hay  be  analyieo  as  tnougn  comprised  or  a  single 
ooase.  The  energy  density  soectra  or  oroad-oano  GOSn 
recordings  are  corrected  ror  oom  tne  oody-wave  attenu- 
itioh  and  tne  freauency  oand  of  tne  recoromg.  The 
corrected  energy  fiu«  measurements  are  then  used  to 
Mtimate  the  raoiateo  seismic  energy. 

This  analysis  is  aoolied  to  digital  recordings  of  tne 
teieseispic  wave-orouos  radiated  oy  me  Hay  2.  1983. 
^oannga.  California,  earthquake  ano  me  Octooer  38, 
*983.  Borah  Peak.  loano.  eartnauaae.  For  the  Coalinqa 
artnqua.e,  Es  ■  ( 1 . 51* . 40 ) xlO< 1  dvne-cm  was  calculated 
’top  si i  P-wave  orouos.  while  the  wave-group  rgcoro- 
«  «  station  COL  returned  E,  .  1. 39*10*1  dyn«-C".  cor 
**  Boran  Peak  eartnouake.  .  (3.07e.4«)«l<)21  dyne-cn 
calculated  from  seven  o-wave  orouos.  ’he  eneroies 
f  ein  estimates  of  tne  aooarent  stress  of  17  ano  9  oars 
o**  the  Coalmga  ano  me  Boran  Peak  earmegg.e.  respec- 
‘|veiy.  ▼nis  difference  aoocars  to  result  from  me 
•hi  *r*nt  tectonic  environment  of  me  events,  where  the 
•  rust  faultino  at  Coal i nog  i$  more  strongly  stressed 
l"*n  the  normal  faulting  at  Boran  Peak. 
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^Nvbrid  Space  and  Spatial  Transform  Domain  Nathod  for 

-a__Camputatioa  of  Sourca  Tima  functions  for  Ceaetal 

^aa-Dimwaa»oaal  Kuptura  Propagation 


^  WIN  (Lament -Doherty  Gaol.  Ob.. 
■••I-  Sci.  of  Coliabia  Uaiwaraity. 
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t^^TNrid  apoca  and  special  transform  method  for  tN 
°f  •0',rM  txm  Ewheciooa  baaed  on  an  iat*< 
•  i  aeuddoa  formulation  of  tbo  sued  boundary  valux 
i^*1"*  t9r  tewlt  diaplacamaata  and  tract iona  was  deve- 
iZT*  d0“iB  •*thod  *•  •"  •daptatton  ol 

of  in  which  tha  diaplacamaata  on  th. 

th*  traeaeiooa  off  tho  fault  are  computed  k) 
^•valving  tho  Groan' a  function  for  a  point  sourca  with 
.  tract iooa  from  all  previous  times  to  tha  ons«t  of 
ciog.  Tha  apacial  transform  domain  method  does  tha 


convolutions  by  Fourier  transforming  in  two  dimensions 
tbn  screened  and  Groan's  functions,  multiplying  them  in 
tha  transform  aomaia,  and  chan  mvarea  transforming 
back  to  tba  .pace  domain.  The  spatial  domain  is  nut h 
oora  accurst*  for  aiamiag  over  tha  tip  of  tha  Green's 
function  conn,  while  tha  tranaform  domain  method  io 
faatar  for  natng  over  tha  raat  of  tba  Croon's 
function  cone. 

Ao  a  chock  of  the  accuracy  of  thia  method  tba  results 
•ra  compared  against  tha  known  results  of  Kostrov  for  a 
circular  crocs  which  grows  without  stopping.  Good 

agreement  is  found . 
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Selamlcity  Simulations  With  a  Baca  and  State 

Dependant  Friction  Uw 

T.  CAO  (Deportment  of  Earth.  Atmospheric,  and 

Planetary  Sciences,  Naseachuaetts  Institute  of 

Technology.  Cambridge.  Naan.  02139) 

Tha  dynamic  nation  and  stability  of  a  single 
*Mr«a  ot  freedom  elastic  system  undergoing 
frictional  slip  is  examined.  Tha  svatom  coaalats 
of  a  eliding  block  connected  to  on  elastic  spring 
which  is  driven  at  a  constant  velocity.  Tho 
frictional  slip  is  governed  by  a  laboratory 
lafarrad  friction  law  proposed  by  Dlstsrich  and 
Kuiao.  Wo  further  extended  tho  solution  to  a 
one-d loans lonal  aaas-aprlag  aodal  with  the  some 
friction  law.  with  this  aodal.  we  can  simulate 
quite  different  behavior  exhibited  at  distinct 
sections  along  a  major  strike-slip  fault,  auch  as 
tho  San  Andrsaa.  Oo  same  sections  of  s  foult. 
there  is  little  oalsnlc  activity;  tho  fault 
•PPoere  to  be  locked.  On  other  sections  of  tha 
fault.  Mall  earthquakes  and  aaalsoic  creep 
relieve  at  loaot  a  fraction  of  the  aceiaaulaclng 
strain.  In  a  previous  seismicity  stoulatloa,  wa 
used  a  displacement  hardening-softening  friction 
law  propoaad  by  Stuart  and  found  that  when  we 
chose  a  critical  displacement  of  the  same  order 
of  nagnltude  aa  the  slip-weakening  critical 
displacement  estimated  by  Papegeorglou  end  4*1 
froa  strong  notion  data,  ws  got  s  normal 
seismicity  pattern  including  small  to  largo 
osrt hquakas .  In  the  present  approach,  wo  chess 
thoaa  nood loons lonal  constants  in  cha  friction 
law  according  to  experimental  rooults  and  found 
that  tha  charactartstlc  sliding  dlstanca.  as  it 
is  called  in  cha  friction  law,  la  consistent  with 
tho  slip-weakening  critical  displacement 
determined  by  Pspageorglou  and  4*1 . 
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**1  ■PH-*0<I<  Shear  Crsc.  Proososfno  Between  Q ay  1*1  Oh 
bid  S-wsve  •a'QC’fs  — - 

4"WEwS  (U.S.  Geoloeical  Survey,  Henlo  Park.  C* 


4  Oura  morp  JJ  ( i n_0 : any •  *«ear  crac*  cannot  O'OOa- 
o«»e  ioont»"»ouiiv  at  a  velocity  o»twe»«  tne  Sav'eion 
velocity  f  i”*  S-wave  velocity  /J  .  Out  a  3-0 
or  ?-D  mi seo-mode  s near  crack  can  orooaoate  in  this 
ranoe  of  velocities.,  n*in«  driven  Dv  tn«  moo#  III 
(*nti-olane‘  comoonont.  Thp^y  has  Dee*  snrq  disao roe- 
mqnt  in  the  literature  concermno  tne  solution  ror 
tne  mode  II  component  in  this  Cdse.  Tnq  analytic 
e*oression  that  is  tne  solution  for  cracs  velocity 
less  than  J  ,  aooiieo  to  velocity  creator  man  f  , 
nas  a  negative  stress  intensity  ractor,  whicn  would 
imoly  that  energy  is  generated  at  the  crack  tip. 
Burnaoe,  on  •»>*  other  hand,  round  that  a  difrerent 
analytic  expression,  having  a  stress  intensity  ractor 
o»  ier o  (a  singularity  of  oroer  smaller  tnan  1/2),  is 
•oolicaoie  for  crack  velocity  between  y  ino  & 
Yumencal  calculations  or  3-u  cracxs  have  oeen  a  one 
by  a  hounflarv  inteoral  method.  Thg  numerical  result 
for  a  mode  II  crack  prescribed  to  orooagate  at  con¬ 
stant  velocity  aorees  with  Bumooe  s  solution.  Soon- 
taneous  orooaoation  of  a  *i«eo-aiooe  rMoture  has  been 
cdlculateo  with  a  slio-wea«enma  friction  law,  in 
wnicn  the  siio  velocity  vector  i$  colli near  witn  tne 
total  traction  vector  and  tne  maomtufle  of  traction 
aeoenos  on  *ne  maomtude  of  siio.  The  soootaneous 
ruOture  velocity  can  n#  between  f  and  A  ,  in  wniCn 
case  enerov  is  absorbed  from  botn  components  of  tne 
•lastic  field,  but  tne  cohesive  eoeroy  is  taken 
primarily  from  the  mooe  II!  component. 
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Fault  Interaction  and  ehe  Frequency  of 
Aftershock  Occurrence  1 

T.  Y AHA SUITA 

t.  KJiOPOFF  (both  at  Ioetltute  of  Ceophyalca  aad 
Planetary  Physics,  Univerelty  of  California  Loa 
Angelas.  CA  90021) 


Tho  aftershock  aource  modal*  are  theoretically 
studied  la  detail.  The  effect  of  aource 
complexity  la  fully  introduced  la  these  aource 
•odela.  u  la  shown  that  tho  Omorl  formula  of 
aitarahock  activity  la  vary  w*U  axpiaload  by  our 
sourca  models. 


This  page  mav  be  freely  copied. 


TVo  kinds  of  sourca  models  era  considered.  Xn 
#n*  •°dal  wa  aaauma  tha  axlatooca  of  uobrokan 
patches  on  a  main-shock  fault  plaaa.  Thus* 
pacchoa  ara  aasumad  to  be  the  source  of 
aftarahocka.  la  another  sourca  modal  wa  coneldar 
tha  interaction  of  maio-ohock  fault  aad  its 
surround  log  small  faults.  Tha  catastrophic 
:oalaacanco  of  aalo-ehock  fault  and  each  small 
fault  can  ba  •**  aource  of  aftarahocka  la  thla 
oodal.  Afte-  «;k  act ivl ties  la  tha  interior  and 
exterior  of  the  main-shock  fault  plan*  era 
described  by  th*  former  and  latter  source  model*, 
respectively. 

Power  law  distributions  ar*  assumed  for  tho 
distances  ovor  which  catastrophic  and  quaalataclc 
rupture  ascensions  occur.  Th#  assumption  of  powmr 
law  distribution  is  justified  by  soma  axrarlmontal 
and  aolsmologlcal  observations.  Tho  atraoa 
corrosion  cracking  la  assumed  as  tha  physical 
oachaniam  of  quaalatatlc  ruptura  ascension. 
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a»d  Tectonic  Regime 


of  Nagnlcuds, 


Duko  University, 

Durham,  NC  27704 

Aftarahocka  ar*  associated  with  almost  all  aalsmlc 
•vonts.  Protarov  and  Dslawonskl  (1982)  hav*  datarmiaad 
a  productivity  or  influence  function  and  axamlnmd  its 
behaviour  as  a  function  of  magnitude  and  location. 

Sine#  theoretical  arguments  Indicate  that  aftershock 
production  is  s  reflection  of  tho  physical  processoa 
and  stress  distribution  at  tha  source,  it  Is  not 
unreasonebl*  to  find  substantial  variations  with  depth 
and  location  according  to  a  plat*  tectonic  framework. 

H°AA  •plc«nt#r  cataloguoa  for  tho  veers 
1970-1983  for  main  events  with  Ns  or  n^,  *  6.0  reveal 
that  plots  of  #  of  aftershocks  versus  Ha-*k  dislays  far 
mors  correlation  than  do  similar  plots  of  numoer  against 
Ns  or  mv,  Since  NaWmj,  indicates  ’slow'  or  non- impulsive 
•vent*,  it  is  evident  that  such  an  event  differs  in 
physical  processes  from  ocher  normal '  events.  However, 
the  feet  thee  theae  deviations  from  the  norm  lie  on  both 
sides  of  the  standard  line  would  argue  that  there  are 
two  cleeses  of  'slow'  events,  neither  of  which  is 


me  i . 

Examination  of  the  specifics  ot  these  events  end  ell 
other  well  situated  events  that  are  substantially  dif¬ 
ferent  froa  che  mean  events  of  their  class  ravaal  s 
tectonic  influence  though  conclusions  ere  necessarily 
restricted  by  the  amount  of  data. 


Proisrov  and  Diiowonskl  (1982)  A  Nethod  of  Studying 
Vsmtloni  in  th*  Clustering  Property  of  Eqs. 

J.  Geophys.  Res.  87  .  2829-2839. 
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►  fan iim  Focal  Depth  vs.  Megnituds  ( Moment )  For 
fcartnquaxas  in  eastern  North  America 

ALEXANDER,  a.S. .  (Department  of  Geosciences.  Th* 

Pennsylvania  Stats  University,  University  Park.  PA 

ib802)" 

sssed  on  che  entire  available  record  of  both 
historical  and  inscrumentally-record*  -arthqua*es, 
there  appears  to  be  no  documented  esse  ot  a  fault 
rupture  breaking  the  earth's  surtsc*  in  sasesm  North 
America.  This  implies  a  minimum  local  depth  as  a 
function  of  magnitude  (or  moment)  to  account  for  the 
increase  in  the  sire  of  che  fault  rupture  area  with 
Increasing  magnitude.  This  minimum  depth  depends  on: 
th*  tvpe  of  faulting  Involved  (generally  thrust  or 
strike-slip  in  sasesm  North  America);  eh*  stress 
drop;  che  distribution  of  deviecorlc  stress  with 
depth;  end  possibly  che  direction  of  propagation  of 
th*  rupture,  especially  for  earthquakes  involving 
dominantly  thrust  faulting.  To  establish  thia  depth 
relationship  quantitatively  and  eaplricslly,  accurate 
focal  depth  determinations  are  required;  nowever, 
reliable,  wall-conscralned  depths  are  not  routinely 
available  for  a  great  aany  of  cha  past  earthquakes  in 
eastern  North  America.  In  this  study  wa  have 
*tt**P**d  to  obtain  a  bound  on  th*  minimum  depth  vs 
magnitude  relationship  by  using  wall-constrslnad 
dapens  dactrmined  by  ochar  workers  and  by  further 
analysis  of  ochar  inatnasentally-racordad  earthquakes. 
Whan  available,  both  pP-P  and  Love  and  Raylaigh  wave 
•P*ctral  observation*  ara  used  to  constrain  th*  depth, 
because  seismograms  from  only  a  few  stations  ara 
required  whan  th*  epicenter  la  well-defined.  Although 
thaea  may  ba  geographical  variations,  our  p re 1 ini nary 
results  suggest  chat  m^  —  and-abov*  *  vanes  ara 
typically  deeper  chan  3  xa,  magnitude  3-and-abov* 
•vants  deeper  than  5  km  and  Mgnlcud*  6-and-abovs 
events  ara  deeper  than  10  km. 

Reproduced  from 
best  available  copy. 
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AFGL/DARPA  REVIEW  OF  NUCLEAR  TEST  MONITORING  BASIC  RESEARCH 
US  AIR  FORCE  ACADEMY ,  6-8  MAY  1985 


TITLE;  ESTIMATES  OF  SOURCE  AND  PATH  CHARACTERISTICS  IN  THE  USSR  AND  NORTH 

AMERICA  USING  SHORT  PERIOD  (Lg)  AND  LONG-PERIOD  SURFACE  WAVE  SPECTRA 

AUTHOR:  Shelton  S.  Alexander,  Department  of  Geosciences, 

The  Pennsylvania  State  University,  University  Park,  PA  16802 


SUMMARY:  For  each  of  a  number  of  Soviet  explosions  two  different  spectral 
approaches  were  used  to  determine  the  average  Q  for  Lg  signals  arriving  in  a 
velocity  window  ranging  approximately  from  3.5  to  3.0  km/sec.  The  mean  spectral 
slope  based  on  all  the  events  analyzed  is  then  used  to  calculate  the  average 
attenuation  coefficient  (QU)-*.  Results  from  individual  stations  of  the 
Grafenberg  array  (6  to  13  stations)  were  averaged  for  each  event  and  compared 
to  the  single  station  (Al)  results.  Additionally,  transverse  component  data 
were  available  for  three  of  the  GRF  stations.  At  MAIO  only  vertical  component 
Lg  data  were  available.  For  the  propagation  path  from  Semipalatinsk  to  GRF, 
the  narrow-band  filtering  method  and  the  noise-corrected  power  spectra  of  both 
vertical  and  transverse  components  give  an  average  Q  for  Lg  of  approximately 
1000.  For  the  path  from  Azgir  to  GRF  the  average  Q  is  around  1100.  For  the 
path  from  Semipalatinsk  to  MAIO  the  average  Q  is  approximately  625.  There  appear 
to  be  no  significant  differences  in  the  results  among  individual  events  for  each 
path,  indicating  that  the  strong  tectonic  release  effects  in  evidence  at  long 
periods  for  some  of  the  Semipalatinsk  events  do  not  Influence  the  Lg  spectra  of 
those  events  appreciably  in  the  frequency  band  .3  to  1.2  Hz,  and  the  best  band  to 
use  to  establish  an  m^L_  vs  sq,  relationship  is  centered  at  .6  to  .7  Hz. 

To  overcome  long-period  non-isotroplc  source  effects  for  Shagan  River  events 
several  very  low  tectonic  release  events  are  used  to  obtain  an  empirical  "pure 
explosion"  Rayleigh  signature  (shape)  which  then  is  incrementally  scaled  and 
subtracted  from  the  total  Rayleigh  wave  signature  of  high  tectonic  release 
events  until  the  best  double-couple  solution  that  also  fits  the  Love  wave 
signals  is  found  (i.e.  when  the  trace  of  the  moment  tensor  is  closest  to  zero); 
the  scaled  signal  is  then  the  appropriate  explosion  contribution.  The 
remaining  "pure  earthquake"  Rayleigh  and  Love  wave  signals  for  each  event  are 
interpreted  in  terms  of  a  best-fitting  equivalent  earthquake  model. 

To  determine  the  approximate  threshold  for  obtaining  source  parameters  from 
long-period  surface  waves  in  an  intraplate  setting,  the  phase  matched  filtering 
method  is  used  to  analyze  the  1982  New  Brunswick  aftershock  sequence  with  the 
well-recorded  mainshock  seismograms  as  the  reference  signals.  The  source 
excitation  function  of  each  aftershock  is  extracted  from  the  observed  signal 
spectrum  by  correcting  for  the  path  transfer  function  and  the  main  shock  source 
excitation.  Preliminary  results  using  this  approach  on  recordings  at  GDSN 
stations  in  North  America  indicate  that  source  parameters  (Including  depth)  can 
be  extracted  for  events  as  small  as  about  m^,  4  to  4.5. 

CONCLUSIONS  AND  RECOMMENDATIONS:  The  conclusions  are  included  above. 

The  principal  recommendations  are:  1)  carry  out  further  evaluations  of  Lg  as 
an  Independent  measure  of  yield,  including  especially  transverse  component  Lg 
observed  at  local  (e.g.  NORESS)  and  more  distributed  arrays  (e.g.  GRF,  NORSAR); 

2)  further  Investigate  source  parameter  estimation  methods  using  3-component 
broad-band  recordings  at  regional  distances;  and  3)  assess  relative  effects  of 
scattering  vs.  absorption  for  Eurasian  and  North  American  Lg  propagation  paths. 


III.  ESTIMATES  OF  SOURCE  AND  PATH  CHARACTERISTICS  IN 

THE  USSR  AND  NORTH  AMERICA  USING  SHORT-PERIOD  (Lg) 
AND  LONG-PERIOD  SURFACE  WAVE  SPECTRA 


ESTIMATES  OF  SOURCE  AND  PATH  CHARACTERISTICS 
IN  THE  USSR  AND  NORTH  AMERICA 

USING  SHORT  PERIOD  <Lg)  AND  LONG-PERIOD  SURFACE  WAVE  SPECTRA 


S.S.  Alexander 
Geosciences  Department 
The  Pennsylvania  State  University 
University  Park,  PA  16802 


Several  approaches  can  be  used  to  obtain  improved  absolute  estimates  of 
yield  from  Lg.  One  is  the  approach  of  von  Seggern  and  Alexander  (1983)  where 
path  effects  are  estimated  and  automatically  taken  into  account  to  obtain  the 
source  moment  (yield)  for  each  event.  Another  is  to  use  a  local  array  such  as 
Grafenberg  or  NORSAR  (NORESS)  to  obtain  the  mean,  noise-corrected  Lg  spectrum 
(with  station  structure  effects  averaged  out)  from  which  the  path  Q  can  be 
estimated,  provided  the  source  spectral  shape  for  Lg  excitation  can  be 
inferred.  Estimates  of  average  Q  of  Lg  signals  in  the  frequency  range  .3  to 
1.2  Hz  have  been  obtained  for  propagation  paths  from  the  Semlpalatlnsk  and 
Azglr  areas  of  the  USSR  to  the  Grafenberg,  West  Germany  (GRF)  array  and  from 
Semlpalatlnsk  to  Mashed,  Iran  (MAIO).  For  each  of  a  number  of  presumed 
explosions  two  different  spectral  approaches  were  used  to  determine  the 
average  Q  for  Lg  signals  arriving  in  a  velocity  window  ranging  approximately 
from  3.5  to  3.0  km/sec.  One  involves  using  a  sequence  of  narrow-band  filters 
at  different  center  frequencies  and  determining  the  spectral  slope  based  on 
both  the  maximum  amplitude  and  the  rms  amplitude  of  the  filtered  signal  in  the 
velocity  window.  The  other  involves  computing  the  spectral  slope  from 
noise-corrected  power  spectra  of  the  same  velocity  window.  The  mean  slope 
based  on  all  the  events  analyzed  is  then  used  to  calculate  the  average 
attenuation  coefficient  (QU)"*.  In  the  case  of  GRF  results  from  individual 
stations  of  the  Grafenberg  array  (6  to  13  stations)  were  averaged  for  each 
event  and  compared  to  the  single  station  (Al)  results.  Additionally, 
transverse  component  data  were  available  for  three  of  the  GRF  stations.  At 
the  MAIO  station  only  vertical  component  Lg  data  were  available.  For  the 
propagation  path  from  Semlpalatlnsk  at  GRF,  the  narrow-band  filtering  method 
gives  an  average  QU  of  2980  and  2940  for  vertical  and  transverse  components, 
respectively.  The  noise-corrected  power  spectra  also  give  a  comparable 
average  QU  of  3190.  These  results  Indicate  an  average  Q  for  Lg  propagation  to 
GRF  of  approximately  1000.  For  the  path  from  Semlpalatlnsk  to  MAIO  the 
average  QU  value  is  approximately  1480  corresponding  to  an  average  Q  of  around 
450.  For  the  path  from  Azglr  to  GRF  the  average  QU  is  approximately  2360 
corresponding  to  an  average  Q  of  around  680.  Examples  of  these  spectral 
results  for  different  events  are  shown  in  Figure  1. 

There  appear  to  be  no  significant  differences  in  the  results  among 
individual  events  for  each  path,  indicating  that  the  strong  tectonic  release 
effects  in  evidence  at  long  periods  for  some  of  the  Semlpalatlnsk  events  do 
not  influence  the  Lg  spectra  of  those  events  appreciably  in  the  frequency  band 
.3  to  1.2  Hz.  In  addition,  plots  of  Lg  spectral  amplitude  at  selected 
frequencies  vs  nq,  for  GRF  and  MAIO  paths  indicate  that  the  best  frequency  band 
to  use  to  establish  an  m^g  vs  m^  relationship  is  centered  at  .6  to  .7  Hz. 
Results  for  this  frequency  range  are  shown  in  Figure  2.  The  Lg  signal  level 
in  this  frequency  range,  corrected  back  to  the  source  using  the  estimated 
(QU)"*  value  appropriate  for  each  path,  as  discussed  above,  gives  an  estimate 
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of  source  strength  (yield).  Similar  analysis  of  NTS  and  Hoggar  events  of 
known  yield  is  in  progress  to  develop  the  appropriate  conversion  of  source 
terms  from  this  method  to  yield. 


To  overcome  the  biasing  perturbations  on  long-period  surface  waves  caused 
by  non-isotropic  source  effects,  a  formal  moment  tensor  decomposition 
(properly  constrained)  can  be  used  to  estimate  the  isotropic  (explosion) 
contribution  if  there  is  sufficient  azimuthal  coverage.  However,  in  many 
instances  surface  wave  data  from  relatively  few  stations  are  available,  so 
other  techniques  are  needed,  especially  for  thrust-type  tectonic  release 
associated  with  Shagan  River  events.  Our  alternative  approach  for  the  Shagan 
River  area  is  to  use  several  very  low  tectonic  release  events  to  obtain  an 
empirical  “pure  explosion"  Rayleigh  signature  (shape)  which  then  i3 
incrementally  scaled  and  subtracted  from  the  total  Rayleigh  wave  signature  of 
high  tectonic  release  events  until  the  best  double-couple  solution  that  also 
fits  the  Love  wave  signals  is  found  (i.e.  when  the  trace  of  the  moment  tensor 
is  closest  to  zero).  This  scale  factor  times  the  "pure  explosion”  signature 
then  gives  approximately  the  equivalent  Rayleigh  wave  for  the  explosion 
without  tectonic  release.  For  each  event,  these  derived  “pure  explosion" 
Rayleigh  waves  for  the  network  stations  are  then  used  to  estimate  Mg  (yield). 
The  assumptions  in  this  approach  are  (a)  that  “pure  explosion"  Rayleigh 
waveform  shapes  are  the  same  for  events  of  comparable  yield  in  the  same  source 
region  and  (b)  the  long-period  Rayleigh  and  Love  wave  signatures  generated  by 
non-isotropic  source  effects  can  be  represented  by  a  double-couple  earthquake 
source.  The  remaining  "pure  earthquake"  Rayleigh  and  Love  wave  signals  for 
each  event.  Interpreted  in  terms  of  a  best-fitting  equivalent  earthquake  model 
(strike,  dip,  rake,  depth,  moment),  characterize  the  corresponding  tectonic 
stress  field  at  the  source.  To  the  extent  that  similar  mechanisms  are  found 
for  events  located  in  sub-areas  of  the  test  site,  a  further  iteration  can  be 
made  to  correct  for  the  tectonic  release  contribution  to  the  Rayleigh  waves 
observed  for  the  low-tectonic  release  calibration  events. 

In  order  to  determine  the  approximate  threshold  for  determining  source 
parameters  from  long-period  surface  waves  in  an  intraplate  setting,  the  phase 
matched  filtering  method  is  used  to  analyze  the  1982  New  Brunswick  aftershock 
sequence  by  using  the  well-recorded  mainshock  seismograms  as  the  reference 
signals.  This  method  not  only  significantly  increases  S/N  for  the  aftershock 
records,  but  also  separates  primary  Rayleigh  or  Love  wave  trains  from 
multipaths  and  other  contaminations.  The  well-documented  mainshock  source 
parameters  have  been  used  to  synthesize  the  source  excitation  function  of  the 
reference  filter.  Hence,  we  obtain  the  path  transfer  function  by  division  of 
the  observed  mainshock  signal  spectrum  at  each  station  by  the  complex-valued, 
mainshock  source  excitation  function  for  that  source-station  azimuth.  The 
source  excitation  function  of  each  aftershock  can  then  be  extracted  from  the 
observed  signal  spectrum  through  division  by  the  path  transfer  function. 
Finally,  the  resulting  Love  and  Rayleigh  source  spectral  excitation  vs 
frequency  and  azimuth  Is  used  in  an  Inverse  calculation  to  derive  the  strike, 
dip,  slip,  depth,  and  moment  of  each  aftershock.  Preliminary  results  using 
this  approach  on  recordings  at  GDSN  stations  in  North  America  indicates  that 
source  parameters  (Including  depth)  can  be  extracted  for  New  Brunswick 
aftershocks  as  small  as  about  mj,  4  to  4.5.  Figure  3  and  Table  l  give  examples 
of  results  for  four  of  the  New  Brunswick  events. 
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Figure  >2b  Observed  Lg  amplitude  vs  m  for  the  center  frequency  Indicated 
for  Semlpalatinsk  explosions.  Left-hand  figure  is  the  maximum 
amplitude  and  the  right-hand  figure  is  the  RMS  in  the  Lg  velocity 


Figure  2c  Observed  Lg  amplitude  vs  m^  for  the  center  frequency  indicated 
for  Semipalatinsk  explosions.  Left-hand  figure  is  the  maximum 
amplitude  and  the  right-hand  figure  is  the  RMS  in  the  Lg  velocity 
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ABSTRACT 

The  1982  New  Brunswick  earthquake  series,  including  the  m^,  5.7 
mainshock  and  several  high  magnitude  aftershocks,  posed  interesting 
problems  in  determining  source  mechanism  and  focal  depth.  The  high 
quality,  long  period  and  broad-band  DWWSSN  digital  recordings  from  SCP 
Inspired  a  single  station  surface  wave  study  to  address  these  questions. 
Band-pass  filters  applied  over  the  entire  frequency  range  provided  group 
velocity  and  spectral  amplitude  data  for  Rayleigh  and  Love  inodes.  Syn¬ 
thetic  amplitude  spectra  were  generated  assuming  structure  parameters 
common  to  both  the  source  area  and  the  propagation  path.  Key  aspects  of 
the  source  mechanism  could  be  Inferred  from  the  fundamental  mode  observa¬ 
tions,  particularly  for  the  mainshock  on  January  9,  1982,  and  the  largest 
magnitude  aftershock  which  occurred  on  January  11.  Although  there  are 
significant  differences  In  modal  excitation  at  the  higher  frequencies 
from  event  to  event,  only  qualitative  support  for  source  mechanism 
could  be  Inferred  because  higher  order  modes  could  not  be  identified. 
Additionally,  published  NEIS  teleseismic  m{j,  network  n^g,  and  Ms  values 
were  compared  with  m^g  and  Ms  values  computed  with  the  SCP  data.  These 
results  Indicated  that  the  single-station  m^g  values  for  the  four  New 
Brunswick  events  studied  agreed  remarkably  well  with  the  network  average 
magnitudes,  whereas  Ms  values  commonly  disagreed.  In  spite  of  the 
observed  variable  excitation  of  the  individual  higher  order  modes, 
consistent  with  the  theoretical  excitation  appropriate  for  the  different 
mechanisms  for  these  events,  the  ensemble  of  modes  that  comprises  the 
Lg  signal  gives  every  indication  of  providing  a  stable  measure  of 


source  strength  when  m^Lg  measured  according  to  Nuttli's  magnitude 
relationship. 


vi 


The  fundamental  mode  Love-to-Rayleigh  surface  wave  spectral  ratios 
for  periods  greater  than  ten  seconds  provided  useful  qualitative  evidence 
for  source  mechanism  and  focal  depth.  The  mechanism  proposed  by  Choy 
et.al.  (1983)  was  supported  for  the  mainshock,  with  fault  plane  striking 
at  195  degrees  and  dipping  65  degrees  to  the  northwest  with  a  rake  of 
70  degrees,  but  with  preferred  focal  depth  of  7  km  rather  than  the  9 
km  given  by  Choy  et.al.  (1983).  The  results  for  the  aftershock  on 
January  11  are  wel 1-expl ained  by  the  mechanism  proposed  by  Wetmiller 
et.al.  (1983),  indicating  thrusting  on  a  plane  conjugate  to  the  fault 
plane  of  the  mainshock,  with  strike  of  332  degrees,  dipping  48  degrees 
to  the  southeast  with  a  rake  of  59  degrees.  The  focal  depth  for  this 
event  was  found  to  be  6  km  from  Love-to-Rayleigh  spectral  ratios  at 
SCP.  While  results  for  the  aftershock  on  January  9  were  less  con¬ 
clusive,  they  suggest  a  shallower  depth  than  the  mainshock.  Although 
the  thrusting  mechanisms  are  consistent  with  eastern  North  American 
crustal  stress  patterns,  the  inferred  dip  and  rake  of  the  mainshock 
are  not  consistent  with  faulting  under  the  action  of  the  ENE  to  NE 
horizontal  compressive  stress  prevailing  over  much  of  the  northeastern 
United  States.  The  inferred  source  parameters  for  the  January  11 
aftershock  are  consistent  with  this  stress  orientation. 
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CHAPTER  1 
INTRODUCTION 

The  recent  series  of  New  Brunswick  earthquakes,  beginning  with  the 
January  9,  1982,  mainshock  (mb  «  5.7),  has  provided  high  quality  data 
from  which  Inferences  on  source  mechanism  and  crustal  structure  may  be 
drawn.  Interest  in  these  events  has  been  particularly  great  due  to  the 
unusually  large  mainshock;  most  previous  studies  in  eastern  North 
America  involve  earthquakes  of  magnitudes  less  than  5  in  ascertaining 
seismicity  trends,  stress  patterns,  and  implications  for  regional  tec¬ 
tonics,  all  Important  considerations  in  assessing  seismic  hazard.  As 
of  its  occurrence,  the  mainshock  was  the  largest  event  in  eastern  North 
America  since  the  Institution  of  the  Global  Digital  Seismograph  Network 
(GDSN),  and,  as  such,  it  offers  an  unusual  opportunity  for  analysis  at 
both  regional  and  teleseismic  ranges.  Resolution  of  source  parameters 
for  all  the  events  from  the  available  body  wave  data  was  complicated 
by  ambiguities  In  depth  determination  resulting  from  depth  phase  identi¬ 
fication  and  noisy  teleseismic  recordings  of  the  aftershocks,  conflict¬ 
ing  evidence  from  locally  monitored  aftershock  activity,  and  the  sugges¬ 
tion  of  conjugate  rupturing  along  Intersecting  fault  limbs.  Signifi¬ 
cantly  different  source  parameters  have  been  Inferred  for  the  mainshock 
by  different  Investigators.  There  Is  even  more  uncertainty  for  the 
larger  aftershocks,  particularly  with  regard  to  focal  depth.  Long 
period  and  broad-band  surface  waves  recorded  at  the  State  College, 
Pennsylvania  (SCP)  station  of  the  Digital  World-Wide  Seismographic 
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System  Network  (DWWSSN)  had  very  high  signal-to-noise  ratios.  The  fun¬ 
damental  mode  Rayleigh  and  Love  wave  signals  exhibited  such  classic  dis¬ 
persion  features,  along  with  a  striking  two-to-one  Love-to-Rayleigh 
excitation  that  surface  wave  analysis  seemed  a  viable  alternate  approach 
to  acquire  source  information  including  depth  of  focus.  The  Lg  contri¬ 
bution  was  also  well-recorded  for  the  mainshock  and  the  larger  after¬ 
shocks,  so  that  source  differences,  reflected  in  the  excitation  of 
these  higher  order  modes,  could  be  investigated.  This  is  of  interest 
because  mbl_g  at  SCP  agrees  closely  with  the  magnitudes  obtained 
from  the  total  set  of  regional  and  teleseismic  observations  for  these 
events. 

Directions  of  the  Present  Study 

Determination  of  focal  depth  using  surface  wave  analysis,  in 
particular,  surface  wave  spectra  and  Love-to-Rayleigh  spectral  ratios, 
has  been  addressed  by  a  number  of  investigators  (Harkrider,  1970;  Tsai 
and  Aki,  1970a,  1970b;  Canitez  and  Toksoz,  1971;  Masse  et.al.,1973; 
Turnbull,  1976;  Mendiguren,  1977;  Mitchell  and  Herrmann,  1979; 
Cybriwsky,  1979).  Multi-station  studies  reveal  that  azimuthal  varia¬ 
tion  In  spectral  amplitudes  permits  determination  of  source  mechanism 
as  well  as  depth.  Single-station  studies  report  precise  results  on 
focal  depths  when  fault-plane  solutions  are  well  known.  Although 
sensitive  to  the  accuracy  of  all  source  parameters  and  models  chosen 
for  earth  structure,  precise  focal  depth  determination  via  surface  wave 
analysis  Is  viable  because  the  spectra  are  so  strongly  affected  by 
depth  of  source.  Attenuation,  crustal  thickness,  spatial  and  temporal 
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behavior  of  the  source  are  complicating  factors  which  affect  absolute 
amplitude  spectra,  but  which  do  so  to  a  far  smaller  extent  than  does 
depth  (Tsai  and  Aki,  1970a)  and  which  can  be  cancelled  out  when  Love-to- 
Rayleigh  amplitude  ratios  are  used  (Canitez  and  Toksoz,  1971).  Given 
the  shallow  nature  of  the  New  Brunswick  events,  significant  surface 
wave  excitation  was  expected.  The  digital  data  at  SCP  -  especially 
for  the  mainshock  -  quite  obviously  contained  meaningful  surface  wave 
Information,  the  usefulness  of  which  compounded  substantially  with  the 
number  of  modes  that  could  be  identified.  It  was  hoped  that  a  single¬ 
station,  multi-mode  surface  wave  analysis  of  the  SCP  data,  incorporating 
source  parameters  available  from  existing  body  wave  studies,  would 
provide  additional  constraints  on  source  mechanism,  especially  focal 
depth,  for  these  New  Brunswick  events. 

Spectral  amplitude  of  surface  wave  modes  observed  at  a  recording 
station  depend  not  only  on  the  depth  and  focal  mechanism  of  the 
earthquake  source,  but  also  on  the  elastic  and  anelastic  properties  of 
the  crust  and  upper  mantle  between  the  source  and  the  receiver  and  on 
the  instrument  response  of  the  system  recording  the  signal;  these 
propagation  and  instrument  effects  must  be  removed  from  the  observed 
signals  before  conclusions  can  be  drawn  about  any  single  source  para¬ 
meter.  This  Is  far  from  a  trivial  operation,  since  different  combina¬ 
tions  of  source  and  propagation  effects  can  in  some  Instances  produce 
very  different  signals.  The  special  advantage  associated  with  surface 
waves  is  that,  for  a  given  path  structure  and  receiving  Instrument, 
excitation  of  each  of  the  various  fundamental  and  higher  modes  of 
Love  and  Rayleigh  waves  is  extremely  sensitive  to  source  parameters. 
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Including  strike,  dip,  rake  and  focal  depth.  In  principle,  each 
mode  that  can  be  resolved  from  the  observed  surface  wave  seismogram 
and  compared  by  spectral  amplitude,  frequency  by  frequency,  provides  an 
independent  piece  of  evidence  with  which  to  evaluate  the  source.  In 
practice,  this  can  be  a  challenging  undertaking;  however,  relative  modal 
excitation  among  events  can  be  used  more  easily  to  characterize 
similarities  and  differences  in  the  sources. 

The  New  Brunswick  earthquakes  provided  an  unusual  opportunity  for 
surface  wave  analysis  of  source  mechanism  and  depth,  since  the 
epicenters  of  the  ten  larger  events  recorded  at  SCP  were  all  within  a 
few  tenths  of  a  degree  in  latitude  and  longitude  from  one  another.  This 
ensured  that  the  paths  traveled  from  source  to  receiver  were  essen¬ 
tially  identical  and,  thus,  dispersion  and  attenuation  effects  were 
assured  to  be  the  same  for  each  surface  wave  mode.  Short,  intermediate, 
and  long  period  digital  data  at  SCP  were  all  used  in  the  analysis. 

Since  transfer  functions  for  each  instrument  had  been  successfully 
developed  to  deconvolve  amplitude  and  phase  effects  introduced  by  the 
recording  instruments,  the  full  spectral  range  recorded  could  be  used  to 
compare  source  excitation  among  the  events. 

The  analysis  and  Interpretation  approach  that  was  used  is 
summarized  as  follows: 

1.  By  applying  multiple  band-pass  filters  to  the  instrument- 
corrected  data,  group  arrivals  associated  with  the  center  frequency  of 
each  band-pass  were  Identified  and  associated  with  individual  modes,  and 
group  velocities  were  calculated  from  the  arrival  times  of  envelope 
peaks.  Theoretical  dispersion  curves  for  crustal  structures  represent- 


5 


1 ng  reasonable  average  paths  between  New  Brunswick  and  SCP  were  compared 
with  the  band-pass  dispersion  results,  and  the  best-fitting  model  was 
selected  to  represent  the  structure  in  subsequent  calculations  of 
theoretical  seismograms.  A  source  structure  appropriate  for  the 
New  Brunswick  epicentral  area  was  used  to  calculate  the  modal  excita¬ 
tion  for  various  source  mechanisms  and  focal  depths.  In  this  case, 
the  best-fitting  path  structure  was  also  appropriate  for  the  source 
structure,  so  the  two  were  chosen  to  be  the  same. 

2.  Programs  provided  by  D.  Harkrider  were  adapted  for  use  on  the 
IBM  370/3081  at  PSU  (under  the  H-extended  compiler  option)  to  compute 
multi-mode  Love  and  Rayleigh  wave  excitation  and  synthetic  surface  wave 
seismograms.  After  inputting  eigenvalues  of  the  preferred  multilayer 
structure  and  assuming  various  source  mechanisms  obtained  in  previous 
studies  of  the  New  Brunswick  events,  amplitude  spectra  were  generated 
for  fundamental  and  higher  order  Love  and  Rayleigh  modes  for  varying 
source  depths.  Comparison  of  the  predicted  excitations  with  the  rela¬ 
tive  amplitudes  of  observed  envelope  peaks  helped  identify  specific 
modes  in  the  data.  Since  the  mechanism  for  the  January  11,  1982, 
aftershock  was  postulated  to  be  different  from  the  mainshock,  synthetic 
spectra  were  generated  using  source  parameters  proposed  by  Wetmiller 
et.al.  (1983)  In  an  attempt  to  substantiate  this  mechanism  and  to  use  the 
predicted  amplitude  excitation  to  resolve  the  depth.  With  the  evidence 
accumulated  from  analysis  of  the  January  sequence,  application  of  the 
same  technique  to  the  data  from  the  June  16,  1982,  earthquake  permitted 
Inferences  of  Its  depth  and  source  mechanism.  Since  it  is  located  about 
30  km  west  of  the  other  events,  and  no  earthquake  activity  has  been 
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detected  in  the  Intervening  region  (Wetmiller  et.al.,  1982,  1983),  its 
relationship  to  the  others  is  of  interest,  particularly  with  regard 
to  Implications  for  further  seismic  activity  in  New  Brunswick. 

The  objectives  of  this  study  have  been:  (a)  to  use  multi-mode 
surface  wave  analysis  to  develop  additional  evidence  for  constraining 
focal  depth  and  source  mechanism  for  the  New  Brunswick  events;  (b)  to 
delineate  arrival  times  versus  frequency  for  higher  order  modes  from 
band-pass  filter  results,  in  hopes  of  providing  more  complete  group 
velocity  dispersion  curves;  (c)  to  determine  the  extent  to  which 
conventional  single-station  Lg  magnitudes  (m^g)  give  good  estimates 
of  source  size  when  individual  higher  mode  excitation  varies  signifi¬ 
cantly  from  event  to  event  in  the  New  Brunswick  sequence;  and  (d)  to 
characterize  the  substantial  differences  in  path  effects  on  Lg  in 
eastern  United  States  compared  to  the  western  United  States. 
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CHAPTER  2 

GEOLOGICAL  BACKGROUND  AND  REGIONAL  SEISMICITY 


Introduction 


The  present  geological  activity  of  the  eastern  United  States 
does  not  offer  the  excitement  of  a  lithospheric  plate  plunging 
under  another  along  a  Benioff  plane  or  the  rapid  grinding  of 
one  plate  against  another  along  a  system  of  transcurrent 
faults;  nor  is  there  active  volcanism  or  continental  rifting. 
No  doubt,  most  of  these  processes  were  vigorously 
operative  in  the  past.  Today,  we  are  left  with  a  battered, 
contorted,  and  creaking  remnant  of  the  continental  lithosphere 
that  evolved  over  a  long  period  of  time  and  that  somehow 
reflects  these  past  processes.  The  challenge  is  to  unravel 
the  events  and  the  influence  of  older  fabric  upon  younger 
structure. 


-  Diment  et.al.,  1966 


Epicentral  data  for  the  New  Brunswick  events  chosen  for  this  study, 
as  reported  in  the  Preliminary  Determination  of  Epicenters,  are  shown  in 
Table  1.  While  this  region  is  clearly  not  aseismic,  it  is  unusual  for 
events  greater  than  magnitude  5  to  occur  there.  Figure  1(a)  presents  a 
record  of  historical  seismicity  for  eastern  North  America  which  can  be 
expected  to  be  subject  to  the  usual  population  bias.  Yang  and  Aggarwal 
(1981)  have  attempted  to  revise  this  record  in  their  plots  of  epicentral 
locations  of  events  of  mf,  greater  than  or  equal  to  two  for  the  period 
1970  -  1979,  as  instrunental ly  recorded  by  the  Northeastern  Seismic 
Network  [Figure  1(b)].  A  more  detailed  picture  of  activity  In  the 
epicentral  region  has  been  presented  by  Wetmiller  et.al.  (1983),  who 
Illustrate  the  regional  seismicity,  distinguishing  events  of  magnitudes 
less  than  3  from  four  documented  earthquakes  of  assigned  magnitudes 


between  4  and  5.5  (Figure  2).  They  report  that,  while  it  is  likely  the 
largest  earthquake  to  occur  in  New  Brunswick  within  historic  times,  the 
January  9  mainshock  is  not  inconsistent  with  existing  probabilistic 
ground  motion  maps  of  Canada  derived  for  seismic  risk  purposes 
(referenced  to  Basham  et.al.,  1982).  From  these  observations,  it  is 
evident  that  the  remote  New  Brunswick  area  is  considerably  more  seis- 
mically  active  than  the  historical  data  indicate;  however,  magnitudes 
on  the  order  of  m^  *  5  still  are  not  common  there. 


Table  1 

Earthquake  Parameters  of  New  Brunswick  Mainshock 
and  Four  Principal  Aftershocks 

(As  Reported  in  the  PDE) 


Date 

Origin  Time 
(h  m  s) 

Latitude 
(Deg.  N) 

Longitude 
(Deg.  W) 

mb 

January  9,  1982 

12:53:52 

46.98 

66.66 

5.7 

January  9,  1982 

16:36:43 

47.02 

66.65 

5.1 

January  11,  1982 

21:41:08 

46.98 

66.66 

5.4 

March  31,  1982* 

21:02:20 

47.00 

66.60 

5.0 

June  16,  1982 

11:42:30 

46.97 

66.99 

4.7 

*  Not  analyzed  in  the  present  study 
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Understanding  the  causes  of  large  eastern  North  American 
earthquakes  has  been  clouded  by  the  fact  that  geologic  and  tectonic 
correlation  with  their  hypocentral  locations  and  focal  mechanisms  have 
In  general  been  lacking  (Herrmann,  1979;  Fletcher  et.al.,  1978). 
Current  hypotheses  tend  toward  one  of  two  extremes:  a)  fractures  in 
previously  unbroken  rock,  or  b)  slippage  on  major  pre-existing  faults 
(reported  in  Choy  et.al.,  1983).  The. first  of  these  explanations  has 
been  linked  tentatively  with  high  stress  drops  that  are  important  for 
consideration  of  strong  ground  motion;  however,  predicting  future 
high  stress-drop  earthquake  locations  via  traditional  geologic  methods 
would  appear  to  be  very  difficult.  The  second  should  be  intimately 
related  to  neotectonic  fault  patterns  and  contemporary  crustal  move¬ 
ments.  Now  that  instrumental  coverage  of  this  area  is  so  greatly  im¬ 
proved,  analysis  of  large-magnitude  events  such  as  the  New  Brunswick 
series,  coupled  with  geological  studies  of  their  epicentral  locations, 
including  high-resolution  seismic  profiling,  may  make  it  possible  to 
determine  conclusively  whether  earthquake  sources  are  associated  with 
new  or  with  pre-existing  fault  structures. 

Geologic  and  Tectonic  Setting 

The  New  Brunswick  earthquakes  occurred  within  a  granitic  pluton  of 
the  Miramichi  Anticl inorium  (Figure  3)  which  is  located  within  the 
Gander  Zone,  one  of  five  major  geologic  provinces  defined  by  sharp 
contrasts  In  the  features  of  the  late-Precambrian  and  early  Paleozoic 
rocks  across  the  Canadian  Appalachians:  in  northwest  to  southeast  pro¬ 
gression,  these  are  Humber,  Dunnage,  Gander,  Avalon  and  Meguma  zones 
(Williams,  1982;  Bird  and  Dewey,  1970).  In  particular,  the  Dunnage, 
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Gander,  and  Avalon  zones  are  likely  extensions  of  the  Connecticut-Valley 
Gaspe  Syncl inorium,  the  Bronson-Hill  Antic! inorium,  and  the  Merrimac 
Syncl inorium.  New  England  counterparts  with  somewhat  different  features 
(King,  1977;  Taylor  and  Toksoz,  1979;  Yang  and  Aggarwal,  1981).  While 
their  tectonic  evolution  is  similar,  these  zones  do  not  share  the  meta- 
morphic  and  igneous  character  of  the  New  England  province  nor  the 
distinctive  fold-fault  features  and  longitudinal  two-fold  division  so 
apparent  in  the  trends  of  the  U.S.  Appalachians  (Spencer,  1969). 

In  accord  with  plate  tectonic  theory,  the  current  thinking  on  the 
evolution  of  the  Northern  Appalachians  centers  on  the  rifting  and  sub¬ 
sequent  closing  of  a  proto-Atlantic,  Iapetus  ocean  during  late- 
Precambrian  through  Paleozoic  time,  with  the  Taconic,  Acadian,  and 
Alleghenian  orogenies  all  associated  with  various  episodes  in  the 
closing  process  (Bird  and  Dewey,  1970;  King,  1977;  Rast  and  Stringer, 
1974;  Taylor  and  Toksoz,  1979,  1982;  Williams,  1979,  1982).  The  Humber 
zone,  with  its  Grenvillian  basement,  constituted  the  stable  continental 
margin  that  was  left  behind  as  the  Iapetus  Ocean  opened  while  the  Dunnage 
and  Gander  zones  comprised  the  oceanic  region.  The  Avalon  and  Meguma 
zones  may  have  been  the  eastward  drifting  continental  mass  (King,  1977), 
although  the  Precambrian  and  Paleozoic  evolution  of  these  zones  is 
unclear  and  allows  a  wide  range  of  explanations  (Williams,  1982).  As 
oceanic  closing  proceeded  during  early  to  mid-Ordovician  through  Per¬ 
mian  time.  Dunnage  and  Gander  zones  were  juxtaposed  with  the  Humber  zone 
and  with  one  another,  producing  at  their  boundaries  the  most  dramatic 
Impacts  of  the  Taconic  orogeny.  Of  particular  Interest  were  the  steep 
faults  with  locally  wide  zones  of  ductile  deformation  produced  in  the 
Gander  zone.  In  New  England,  the  Bronson-Hill  Antic! inorium  experienced 
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considerable  volcanic  activity  during  this  period,  consistent  with 
island-arc  volcanism  with  subduction  strongly  suggested  on  an  eastward¬ 
dipping  Benioff  zone.  Similar  trends  in  New  Brunswick  implied  south¬ 
easterly  subduction  under  the  Gander  zone  (Bird  and  Dewey,  1970;  Taylor 
and  Toksoz,  1979,  1982;  Williams,  1979,  1982). 

After  a  period  of  relative  quiescence  during  Silurian  to  early 
Devonian  time,  tectonic  activity  resumed  with  the  Acadian  orogeny; 
renewed  volcanism,  high-grade  metamorphism  in  New  England,  intrusion  of 
granitic  plutons  with  accompanying  large-scale  recumbant  folding  in 
many  cases,  and  deformation  and  development  of  large  thrust  belts  were 
all  featured,  reaching  peak  activity  in  mid-Devonian  time.  The  con¬ 
trolling  features  of  this  episode  of  major  deformation  remain  enigmatic, 
but  their  impact  affected  the  entire  Canadian  Appalachian  sequence, 
including  the  western  zones  previously  deformed  in  the  Taconic  orogeny 
(Williams,  1982).  King  (1977)  asserts  that  almost  all  major  granitic 
plutons  were  emplaced  during  this  period  of  continental  collision 
between  North  America  and  the  Avalon  block  and  that  their  features  bear 
little  similarity  to  modern  island-arc  counterparts,  thus  disfavoring  a 
subduction  mechanism  for  their  origin  (Williams,  1982).  The  subsequent 
Alleghenian  orogeny  was  thought  to  be  associated  with  the  final  continent 
continent  collision,  but,  since  its  Impact  was  restricted  to  the  Avalon- 
Meguma  boundary,  far  eastward  of  the  region  of  interest,  it  is  not  of 
Immediate  concern  to  this  study. 
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Recent  Seismic  Activity  and  Associated  Geologic  Structures 
in  the  Study  Area 

To  depict  the  epicentre!  region,  Wetmiller  et.al.  (1983)  show  the 
aftershock  activity  from  the  January  and  March  events  superimposed  upon 
a  geologic  map  of  the  area  (Figure  4).  Although  not  included  in  this 
picture,  location  of  the  June  event  (Figure  3)  supports  the  observation 
that  all  of  the  earthquakes  occurred  within  a  single  massive  granitic 
pluton  of  Devonian  age,  which  has  yet  to  reveal  any  indication  of 
previous  deformation.  Preliminary  geological  mapping  and  gravity 
investigations  have  been  carried  out  by  Burke  et.al.  (1982),  yielding 
some  initial  estimates  on  the  size  and  extent  of  this  pluton.  They 
report  that  the  area  enclosing  the  epicenters  of  the  mainshock  and 
aftershocks  is  mainly  underlain  by  granitic  rock  with  small  dioritic 
inliers,  all  part  of  the  extensive  800  km?  "North  Pole  pluton”  which 
Intruded  deformed,  Cambro-Ordovician  metamorphic  rocks.  A  large  gravity 
low  coincides  with  the  position  of  the  pluton,  indicating  a  100  -  200 
kg/m3  density  contrast  between  the  plutonic  and  surrounding  metamorphic 
rocks.  The  low  encompasses  some  smaller  positive  anomalies  suggesting 
that  the  higher  density  diorite  is  more  extensive  at  depth  than  is 
indicated  In  the  surface  exposures.  These  gravity  anomalies,  combined 
with  the  best  available  seismic  evidence  for  depth  of  the  earthquakes. 
Indicate  that  the  events  were  confined  to  the  plutonic  body,  which  is 
thus  Inferred  to  be  more  brittle  and  weaker  than  the  surrounding 
metamorphic  rocks  when  exposed  to  the  same  stress  regime  (Burke  et.al., 
1982). 

The  New  Brunswick  events  occurred  within  an  extensive  granitic 
pluton  and  produced  "the  first  observed  earthquake- Induced  bedrock  break 
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In  eastern  Canada,  but  the  displacement  [at  the  surface]  was  small  (25 
mn)  and  not  obviously  a  primary  movement."  (p.2,  Wetmlller  et.al., 
1983).  Located  on  a  pre-existing  joint  oriented  parallel  to  the 
Inferred  rupture  of  the  mainshock,  the  surface  break  occurred  barely 
within  a  less-prominent  joint  system  of  steeply  dipping  joints  (Burke 
et.al.,  1982).  If  it  could  be  established  that  there  is  a  fundamental 
relationship  among  pre-existing  fault  structures,  the  prevailing  joint 
system,  and  the  earthquake  rupture  geometry  (and  associated  stress 
orientation),  it  would  certainly  go  a  long  way  toward  development  of  a 
coherent  tectonic  picture,  but,  so  far,  the  available  evidence  has 
resisted  such  generalization.  Additional  field  studies,  including 
trenching  in  the  epicentral  area,  are  still  in  progress,  so  some 
definite  conclusions  regarding  the  influence  of  earlier  episodes  of 
deformation  may  become  available  when  these  investigations  are 
completed. 
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CHAPTER  3 

THEORY 


The  historical  development  of  the  theory  for  calculating  surface  wave 
dispersion  and  far-field  time  and  frequency  domain  values  for  a  variety 
of  earthquake  source  types  is  traced  very  completely  in  Wang  (1981)  and 
Masse  et.al.  (1973).  One  of  the  central  figures  in  this  development  is 
David  G.Harkrider  who  developed  the  programs  for  modal  excitation  for  earth 
quake  and  explosion  sources  in  multilayered  earth  models  and  for  generat¬ 
ing  the  corresponding  synthetic  surface  wave  seismograms  that  were  used 
in  this  study. 

Following  his  conventional  notation,  as  revised  to  include  attenua- 
/ 

tion  factors  (Masse  et.al.,  1973),  the  far-field  spectral  amplitudes  for 
Rayleigh  (Ur)  and  Love  (U^)  waves  may  be  written  as: 


.,/  x  c/  m"  >  E.(tu)  e'‘v  -V-" 

Ub(uj)  -  S(uu)  e  A  t(0,h,ai)  -> - - e 


n(»)  -  s< »  ji,;  .  . 'u;  <  e.  h.* )  i^vp' v- 


where  S(uj)  is  the  spectral  source  function,  that  Is,  the  Fourier-trans¬ 
formed  dislocation-time  history  (the  same  for  Rayleigh  and  Love  waves  for 
non-propagating  point  sources);  m«0  for  a  point  force,  m«l  for  a  couple 
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or  double-couple  source;  k  (as  iU/C,  angular  frequency  divided  by  phase 
velocity)  is  the  Love  or  Rayleigh  wave  number;  h  is  the  source  depth; 
f  is  azimuth  from  the  strike  of  the  fault,  as  specified  in  Ben- 
Menahem  and  Harkrlder  (1964);  r  is  the  epicentral  distance  and 
Y  r(u»)  and  Yl(uj)  are  the  coefficients  of  anelastic  attenua¬ 
tion  which  are  frequency  dependent.  The  factors  Er  and  El 
are  defined  separately  for  Love  and  Rayleigh  waves  as  their  amplitude 
response  divided  by  the  square  root  of  their  wave  number,  multiplied  by 
ellipticity  at  the  surface  in  the  case  of  Rayleigh  waves. 


E.=  £oA.k. 


-V. 


\K 


The  amplitude  layer  responses,  Ar  and  Al,  are  dependent  on  wave  type, 
frequency,  mode  order,  and  the  elastic  properties  of  the  multilayered 
medium  through  which  the  waves  pass,  and  are  Independent  of  source  type 
and  depth.  Harkrider's  definitions  for  Love  and  Rayleigh  amplitude 
response  are  based  on  Thomson-Haskell  displacement-stress  vector  elements 
(Haskell,  1953),  terms  related  to  Thomson-Haskell  displacement-stress 
matrices  and  their  products,  and  the  rate  of  change  of  the  period  equa¬ 
tions,  Fr  and  Fl,  with  respect  to  wave  number  as: 


A 


#  * 

G  N  -  L  H 
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Definition  of  the  amplitude  response  functions  in  terms  of  depth 
integrals,  which  are  proportional  to  the  kinetic  energy  densities  of 
their  respective  surface  waves,  provides  an  alternative  method  of  cal¬ 
culation  and  a  useful  check  on  accuracy  of  the  values  determined  by  ma¬ 
trix  methods.  Both  techniques  are  used  in  Harkrider's  program.  Numeri¬ 
cal  instabilities  commonly  appear  in  the  Thomson-Haskel 1  matrix  technique 
as  frequency  or  mode  order  is  increased  to  the  point  that  differences 
calculated  between  rapidly  increasing  exponential  terms  lose  numerical 
significance  as  they  become  too  small  for  the  precision  of  the  machine; 
these  difficulties  were  overcome  in  the  Harkrider  program  by  incor¬ 
porating  the  compound  matrix  formulations  developed  by  Dunkin  (1965) 
and  Thrower ( 1965 ) . 

The  complex,  frequency-dependent  radiation  pattern,  .X*(uj),  is 
dependent  on  the  source  parameters:  dip,  4,  rake.  A,  azimuth,  i, 
and  depth  (h)  (see  fault-plane  geometry  illustrated  in  Appendix  1) 
according  to  the  relationships  given  in  Table  2.  A,  B,  C,  S,  V,  and  W, 
all  frequency-dependent  functions  of  source  depth,  are  defined  in 
terms  of  the  displacement-stress  vector  elements  (the  quantities  in 
brackets).  These  elements  are  determined  by  specifying  the  elastic 
parameters  (compresslonal  and  shear  wave  velocities,  densities,  and 
layer  thicknesses)  for  each  layer  of  an  assumed  multilayered  structure. 
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and  demanding  that  the  wave  equation  be  satisfied  under  boundary  condi¬ 
tions  of  zero  stress  at  the  free  surface,  zero  displacement  at  (effec¬ 
tively)  infinite  depth,  and  continuity  of  stress  and  displacement  at 
each  layer  interface  within  the  medium. 

The  coefficients  of  anelastic  attenuation,  y^Cm)  and  y^ui), 
are  related  to  the  specific  quality  factor,  Q,  as: 

'  2QbU(uj) 

where  Q.  is  mode-  and  frequency-dependent  and  U(u>)  is  the  group  velocity 
for  the  specified  mode  at  the  frequency  u>.  Although  there  is  evidence 
suggesting  that  Q ^  is  not  the  same  for  Love  and  Rayleigh  waves  (Mitchell 
and  Herrmann,  1979),  there  is  some  support  for  assuming  that  these  differ¬ 
ences  can  be  ignored  (Tsai  and  Aki,  1970a),  at  least  for  low  frequencies. 
In  this  study,  where  the  propagation  path  is  relatively  short  (1130  km) 
and  identical  from  event  to  event,  it  has  been  assumed  that  7r  and 
yL  are  not  significantly  different  throughout  the  period  range  over 
which  the  spectral  ratio  is  calculated. 

Since  Up  and  U|_  are  dependent  on,  and  particularly  sensitive 
to,  source  finiteness  effects,  errors  in  assumed  time  and  space 
histories  for  the  source  mechanism  will  affect  predicted  results  for 
these  absolute  spectral  amplitudes  adversely.  However,  assuming 
the  source  factors  to  be  approximately  the  same  for  both  Love 
and  Rayleigh  waves  (Masse,  Lambert,  and  Harkrider,  19  ),  forming 

the  amplitude  ratio  cancels  their  Influence  completely.  Taken 
In  combination  with  the  additional  simplifications  Introduced 
through  cancellation  of  common  distance  and  attenuation  terms,  this 
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advantage  strongly  suggests  the  merit  of  using  Love-to-Rayleigh  spectral 

ratios  in  the  evaluation  of  source  parameters,  particularly  depth.  The 

L/R  spectral  ratio  is  very  sensitive  to  focal  depth  (Ben-Menahem  and 

Harkrider,  1964;  Canitez  and  Toksoz,  1971;  Cybriwsky,  1979;  Hark- 
/ 

rider,  1970;  Masse  et.al.,  1973;  Tsai  and  Aki,  1970a;  Turnbull,  1976; 
and  many  other  recent  studies  of  specific  earthquakes). 

The  flow  of  the  numerical  procedures  used  to  calculate  theoretical 
spectral  amplitudes  is  detailed  in  Appendix  1;  a  summary  of  the  steps 
is,  however,  appropriate  here.  A  routine  (RLSCFRQ),  which  scans  fre¬ 
quency/phase-velocity  space  for  appropriate  roots  in  frequency  and 
phase  velocity  by  mode,  greatly  facilitates  the  more  accurate  root¬ 
finding  process  in  the  Love  and  Rayleigh  eigenvalue-generating 
programs.  These  programs  (LOVIT.UP)  and  (RAYLIT.UPEC)  then  calculate, 
for  each  mode  specified,  the  disperion  and  depth- independent  quantities 
for  given  input  frequencies  as  well  as  displacement-stress  values  for 
the  midpoint  depth  of  each  layer.  In  finding  the  roots  of  the  period 
equation,  F(tu,  C,  layer  parameters)  ■  0,  the  Love  wave  program  uses  the 
Thomson-Haskel 1  formulation  of  F|_  while  the  Rayleigh  wave  program 
uses  the  compound  layer  matrix  formulation  of  Fr.  Once  bracketed  to  the 
desired  accuracy  for  the  given  frequency,  successive  roots  (resulting  from 
perturbing  phase  velocity  slightly)  are  used  to  calculate  the  group  velocity, 
U,  the  kinetic  energy  density  and  the  amplitude  response  of  the  medium,  as 
well  as  the  surface  ellipticity  in  the  case  of  the  Rayleigh  program.  The 
group  velocities  and  amplitude  responses  are  computed  by  the  two  different 
techniques  referred  to  earlier,  using  both  partial  derivative  and  energy  in¬ 
tegral  formulations  (Harkrider,  1970).  Output  of  these  programs,  after  minor 
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formatting  adjustments,  are  then  used  as  input  to  the  source  depth 
program  (EQTR.ZBIG)  where  source  type  and  associated  source  parameters 
are  selected,  and  spectral  amplitudes  or  the  seismograms  themselves  can 
be  calculated  (and  plotted)  for  each  mode,  for  a  given  receiver  a2imuth 
and  depth,  for  any  number  of  desired  source  depths.  It  is  important  to 
point  out  that,  at  present,  the  program  assumes  that  the  structures  are 
the  same  at  the  source  and  over  the  propagation  path  to  the  receiver. 
Modification  of  the  program  to  accommodate  excitation  calculations  at 
the  source  using  one  structure  and  propagation  effects  to  the  receiver 
using  another  structure  is  relatively  straightforward.  In  this  study, 
however,  the  chosen  source  structure  fortuitously  also  adequately  matched 
the  observed  path  dispersion  effects  based  on  empirical  band-pass  group 
velocity  determinations.  Therefore,  the  synthetics  generated  here  assume 
an  average  structure  from  New  Brunswick  to  SCP  that  is  the  same  as  the 
source  structure.  Using  various  proposed  source  parameters  obtained  from 
the  fault-plane  solutions  and  moment  tensor  inversions  of  recent  studies 
of  the  New  Brunswick  events,  the  main  program  was  run  to  generate 
theoretical  Rayleigh  and  Love  mode  spectra  and  spectral  ratios  for 
source  depths  ranging  from  1  to  10  km.  These  were  compared  with  the 
band-pass  amplitude  results  In  an  attempt  to  identify  a  preferred  depth 
for  each  of  the  events. 
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CHAPTER  4 

DATA  PROCESSING  AND  ANALYSIS 


Acquisition  of  the  Data 

The  original  DWWSSN  short,  intermediate  and  long  period  seismograms 
for  the  four  events  described  in  Table  1  were  subset  from  the  SCP  digital 
tapes;  they  are  shown  in  Appendix  2.  The  long  period  data  were  continu¬ 
ously  recorded  at  one  sample  per  second.  Short  and  intermediate  period 
signals  were  triggered  and  sampled  at  twenty  and  ten  samples  per  second, 
respectively.  Only  the  vertical  component  is  available  for  the  short 
period;  vertical  and  two  horizontal  components  are  available  for  inter¬ 
mediate  and  long  period  data,  which,  along  with  epicentral  coordinate 
Information,  allow  rotation  to  obtain  radial  and  transverse  components. 

With  reference  to  the  seismograms  in  Appendix  2,  it  is  evident  that, 
at  an  epicentral  distance  of  1130  km  (  a  «  10.2*),  SCP  recorded  a  short 
period  signal  for  the  mainshock  that  clipped  for  about  fifty  seconds 
while  the  Intermediate  period  signal  did  not.  Ratios  of  peak  amplitudes 
measured  on  transverse  and  vertical  components  for  both  Intermediate  and 
long  period  records  showed  a  high  love-to-Rayleigh  excitation;  the  tove- 
to-Rayleigh  ratio  was  nearly  three  on  the  intermediate  period  and 
approached  two  on  the  long  period  recordings.  Classic  dispersion 
features  characterize  the  long  period  signals,  suggestive  of  pure  funda¬ 
mental  mode  arrivals.  Although  substantially  noisier,  seismograms  for 
the  January  9  aftershock  exhibit  similar  features,  particularly  in 
Love-to-Rayleigh  excitation  on  Intermediate  and  long  period  signals. 
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The  January  11  aftershock  displays  considerably  different  waveform  fea¬ 
tures  on  the  long  period  components  with  a  lower  Love-to-Rayleigh  ratio, 
suggestive  of  a  different  source  mechanism  than  the  mainshock.  The 
earthquake  which  occurred  on  June  16,  of  substantially  lower  magnitude 
than  the  other  events,  has  Love-to-Rayleigh  ratios  which  agree  more 
closely  with  the  January  11  aftershock  than  with  the  other  events. 

Instrument  Correction 


In  preparation  for  applying  band-pass  filters,  the  digital  seismo¬ 
grams  were  transformed  using  the  Cooley-Tukey  (1965)  algorithm  for  the 
Fast  Fourier  Transform  (FFT),  yielding  the  complex-valued  frequency 
domain  representation  of  the  time  domain  seismograms.  Choosing  a  power 
of  two  at  least  one  greater  than  that  associated  with  the  number  of  data 
points  sampled  reduced  cyclic  convolution  effects  to  a  minimum  when  the 
filters  were  applied.  Instrument  correction  was  accomplished  In  the 
frequency  domain  by  using  the  complex-valued  transfer  function, 
developed  using  poles  and  zeros  given  in  the  Directory  of  World  Digital 
Seismic  Stations  for  each  of  the  SCP  digital  Instruments.  The  transfer 
function  was  defined  as: 

A*  •  OS  •  (S  -Zj)  •  (S  -  Z2) . (S  -  Z,n) 

T(S)  ■  _ _ _ _ _ _ 

(S  -  Pi)  •  (S  -  ?i) . (S  -  Pn) 

where  S  ■  1  uj 

u j  ■  angular  frequency 

A$  *  a  scalar  quantity  -  different  for  each  Instrument 


DS  *  digital  sensitivity  expressed  in  counts  per  micron 
m  *  the  number  of  complex  zeros 
n  *  the  number  of  complex  poles 
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Table  3  gives  the  relevant  constants  for  each  instrument. 

Deconvolution  of  the  instrument  response  using  this  transfer  function 
is,  then,  a  straightforward  complex  division  in  the  frequency  domain.  Depen¬ 
dent  as  it  is  on  frequency,  the  function  can  accommodate  any  sampling  in¬ 
terval  required  by  the  time  window  of  the  chosen  seismogram.  These  correc¬ 
tions  were  applied  to  the  SCP  data  in  this  study  and  can  now  be  routinely 
applied  to  any  acquired  short,  intermediate,  and  long  period  digital  data 
as  an  option  within  the  bandpass  program. 

Band-Pass  Filtering 

The  band-pass  method  is  a  filtering  procedure  pioneered  in  its 
application  to  multi-mode  surface  wave  analysis  by  Alexander  (1963), 
further  developed  by  Dziewonski  et.al.  (1969)  and  tested  with  regard  to 
group  velocity  and  amplitude  determination  by  Herrmann  (1973)  and 
Mitchell  (1973),  among  others.  The  method  exploits  the  simple  fact 
that  energy  associated  with  different  surface  wave  modes  comprising  a 
dispersed  wave  train  generally  arrive  at  different  times  for  each  fre¬ 
quency;  if  only  a  single  station  (rather  than  an  array)  is  available, 
the  method  fails  when  more  than  one  mode  arrives  at  the  same  time  with 
the  same  frequency.  In  principle,  by  nulling  out  all  but  a  very 
limited  range  of  frequencies  around  a  chosen,  central  value,  the  narrow 
band-pass  filter  gives  the  arrival  time  and  spectral  amplitude  of  each 


24 


Table  3 


Complex  Poles  and  Zeros  of  Transfer  Functions 
for  SCP  DWWSSN  Instrunents 


Short  Period 

Intermediate  Period 

Long  Period 

A<*> 

1.11  x  10? 

354.8 

.02755 

DS 


(counts/ 

micron) 

1000  125 

500 

zeros 

Z1  *  z2  *  z3  *  Z4  *  Z5  »  0 

for  each  instrument 

poles 

Pi 

(-16.63,  0) 

(-.377,  .1827) 

(-.377,  .1827) 

P2 

(-4.241,  6.432) 

(-.377, -.1827) 

(-.377, -.1827) 

P3 

(-4.241,-6.432) 

(-2.407,  5.804) 

(-.654,  0) 

P4 

(-8.415,  .3575) 

(-2.407,-5.804) 

(-.2318,  0) 

p5 

(-8.415,-. 3575) 

(-5.848,  2.366) 

-.2318,  0) 

Pe 

(-111.1,  111.1) 

(-5.848,-2.366) 

(-.2318,  0) 

P7 

(-111.1,-111.1) 

(-.02108,  0) 

(-.3276,  0) 

Pa 

(-.04442,  .04444) 

(-.02108,  0) 

(-.3276,  0) 

P9 

(-.04442, -.04444) 

(-.3276,  0) 

l10 

pn 

(-.0214,  0) 
(-.0214,  0) 
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Individual  signal  comprising  the  seismogram  at  that  frequency.  It  is 
especially  useful  for  separating  out  the  contributions  of  various  sur¬ 
face  wave  modes,  because  modes  of  different  frequency  typically  arrive 
simultaneously  giving  rise  to  the  complex  wave-forms  associated  with  the 
Lg  portion  of  most  seismograms.  Each  envelope  peak  of  the  band-pass 
filtered  signal  corresponds  to  the  group  arrival  time  of  a  mode  (the 
arrival  time  of  the  energy  of  each  wave  packet  at  that  frequency)  and 
the  amplitude  at  each  peak  is  proportional  to  the  Fourier  spectral  ampli¬ 
tude  of  that  arrival. 

The  uncertainty  principle  in  quantum  theory  relates  particle  momentum 
and  position  through  the  Fourier  Transform;  resolution  of  one  quantity  is 
inversely  related  to  resolution  of  the  other.  Similarly,  there  are  time 
and  frequency  domain  tradeoffs  on  resolution  which  must  be  considered  in 
filter  shape  and  design.  Optimally,  only  a  single  frequency  should  be 
passed;  however,  inverse  transformation  of  a  spectrum  windowed  by  such  a 
narrow  filter  would  produce  a  broad  sinusoid  in  time,  whose  envelope  peak 
would  be  difficult  to  pick  with  accuracy  and  whose  interferences  with  the 
carrier  of  nearby  arrivals  could  be  easily  mistaken  for  additional  group 
arrivals.  According  to  the  frequency  scaling  property,  broadening  the 
frequency  domain  function  narrows  the  corresponding  time  domain  function 
in  inverse  proportion.  By  using  a  somewhat  larger  filter  bandwidth  and 
introducing  a  cosine  taper  to  the  filter  in  the  frequency  domain  to 
minimize  sldelobes,  these  spurious  effects  are  reduced.  Group  arrivals 
are  nearly  Invariant  under  changes  in  filter  bandwidth.  The  filter  that 


FIC 


FHC 


was  used  is  defined  as: 

in  shape: 
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=  -< 


^  1 
cos 


from  FLC  to  FHC 
I  Tt  (  F  -  FIC)  { 

I  2  (FLC-  FL)  j  , 


the  bandwidth 

FL  <  F  <  FLC 


V.  cos 


FHC  <  F  <  FH 


Any  desired  filter  shape  could  have  been  used.  For  example,  a 
Gaussian  filter  has  the  advantage  of  not  having  sidelobes,  but  it  is 
susceptible  to  spectral  leakage  when  the  spectrin  being  filtered  is 
highly  colored  as  is  typical  of  individual  modal  excitation  spectra. 

The  shape  chosen  with  a  cosine  taper  avoids  this  problem  at  the  expense 
of  small  sidelobes  and  preserves  the  true  signal  energy  between  the 
taper  points  (FLC  and  FHC).  It  also  allows  assymmetrical  filter  shapes 
to  be  used  (FL  -  FLC  i  FHC  -  FH). 

Application  of  the  filter  is  accomplished  conveniently  by  taking 
Its  product  In  the  frequency  domain  with  the  real  and  Imaginary  spectral 
components  of  the  transformed  seismogram  after  Instrunent  correction;  it 
Is  a  phaseless  filter.  A  suite  of  these  cosine-tapered,  narrow  band-pass 
filters  was  applied  to  the  short,  intermediate,  and  long  period  data, 
limiting  the  range  of  center  frequencies  to  well  within  Nyquist[l/(2 At)] 
where  At  is  the  sampling  Interval  in  time  In  each  case.  Using  a  com¬ 
putational  trick  consisting  of  nulling  the  negative  frequency  contribu¬ 
tions  and  halving  the  spectral  amplitude  at  DC  and  Nyquist  frequencies, 
the  Inverse  transform  will  be  complex-valued  (the  analytic  signal)  where 
the  real  part  Is  the  filtered  seismogram  (times  0.5)  and  the  Imaginary 
part  is  Its  Hilbert  transform  (times  0.5).  The  envelope  was  calculated 
as  the  square  root  of  the  sum  of  the  squares  of  the  real  and  Imaginary 
parts  of  the  returned  time  series.  Envelope  peaks  were  identified,  and 
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group  velocities  calculated  very  straightforwardly  as  the  distance  tra¬ 
veled  divided  by  the  travel  time  to  each  envelope  peak. 

For  filters  centered  around  lower  frequencies,  approaching 
zero  Hz,  narrowing  the  bandwidth  of  the  filter  becomes  a  necessity. 

Since  the  envelope  amplitudes  of  the  filtered  signals  are  dependent  on 
the  bandwidth,  indeed  the  filter  shape  in  general,  this  effect  had  to  be 
considered.  Ideally,  of  course,  it  is  desirable  to  always  use  the  same 
filter  bandwidth  and  shape.  The  narrow  bandwidth  required  for  low  fre¬ 
quencies  was  Inadequate  for  resolution  of  the  various  modal  contributions 
at  the  higher  frequencies.  The  compromise  chosen  was  to  use  two  filters 
with  identical  cosine  tapers  but  with  a  bandwidth  of  .05  Hz  for  central 
frequencies  greater  than  .1  Hz  and  a  bandwidth  of  .005  Hz  for  those  less 
than  .1  Hz.  Clearly,  this  bandwidth  choice  has  a  significant  effect  on 
the  amplitude  of  the  envelope,  so  comparison  of  absolute  amplitudes  was 
restricted  to  frequency  ranges  over  which  the  filter  shape  was  constant. 
Of  course,  effects  of  the  filter  shape  can  be  cancelled  completely  by 
examining  spectral  ratios  instead.  Figure  5  illustrates  probably  the 
most  extreme  Influence  this  effect  can  be  expected  to  have  on  group 
velocity  and  amplitude  values.  Even  with  the  dispersion  apparent  in 
the  filtered  signal  of  .05  Hz  bandwidth,  the  group  velocity  is 
3.37  km/sec,  as  compared  with  3.27  km/sec  for  the  filter  of  .005  Hz 
bandwidth.  This  suggests  that  the  group  velocity  can  be  expected  to  be 
determined  to  within  .1  km/sec,  even  with  a  factor  of  ten  change  in 
bandwidth.  Amplitude  scaling  as  illustrated  by  the  peak  amplitudes  of 
the  filtered  signals  Is  clearly  a  problem,  however. 

Another  problem  which  affects  both  amplitude  and  group  velocity 
results  is  spectral  leakage  due  to  radically  varying  spectral  amplitudes 
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in  the  region  of  the  central  frequency.  This  can  have  an  obvious  effect 
on  the  amplitudes  and  will  shift  the  peaks  in  the  direction  of  high 
amplitude  contributions,  thereby  affecting  both  the  group  velocity  and 
amplitude  calculations.  The  test  for  this  effect  is  the  instantaneous 
frequency  of  the  band-pass  filtered  signal;  It  should  be  very  close 
to  the  center  frequency  of  the  bandpass.  Pre-whitening  before  filtering 
is  expected  to  alleviate  these  leakage  problems. 

The  results  of  applying  the  band-pass  filters  to  the  short,  inter¬ 
mediate,  and  long  period  seismograms  of  Appendix  2  are  compiled  in 
Appendix  3.  Figure  6  is  an  illustration  of  how  dispersive  effects  can 
be  observed  as  arrival  times  of  the  envelope  peaks  change  for  a  given 
modal  arrival.  These  results  for  the  January  11  aftershock  are  remark¬ 
able  in  that  a  single  mode  seems  to  predominate  throughout  the  suite  of 
filters  centering  on  frequencies  ranging  from  5  to  8.5  Hz.  With  Nyquist 
frequency  at  10  Hz,  it  is  surprising  to  see  such  clean  evidence  of  energy 
from  a  single  mode  arriving  at  such  high  frequencies.  Figure  7  is  a 
sobering  example  of  some  of  the  difficulties  to  be  faced  in  Identifying 
different  modes.  Spurious  peaks  can  be  the  result  of  truncation  effects, 
interfering  waveforms,  or  multipathed  arrivals.  Only  the  filtered  signal 
centered  at  5.5  Hz  bears  any  resemblance  to  the  clean  modal  arrivals  of 
Figure  6.  The  result  of  applying  a  suite  of  filters  to  the  vertical 
component,  intermediate  period  signal  of  the  January  9  mainshock  is 
shown  in  Figure  8.  A  combination  of  simple  dispersed  fundamental  mode 
arrivals  at  the  lower  frequencies  and  interfering  waveforms  of  uncer¬ 
tain  origin  at  higher  frequencies  Illustrates  the  typical  appearance  of 
the  filtered  results. 

Since  the  signals  for  all  events  traveled  identical  paths,  disper- 
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Sion  is  assured  to  be  the  same  among  events;  amplitudes  of  the  arrivals 
were  expected  to  vary  with  modal  excitation  at  the  source  and  resultant 
variations  in  Interfering  waveforms.  In  an  attempt  to  present  the  band¬ 
pass  results  in  a  way  that  might  help  isolate  the  modes,  the  amplitudes 
of  each  envelope  peak  on  a  single  band-pass  were  normalized  to  the  peak 
amplitude  for  that  band-pass.  A  crude  scale  was  then  set  to  distinguish 
among  these  normalized  amplitudes,  symbols  were  assigned  to  correspond 
with  these  values,  and  composite  plots  of  group  velocity  vs.  frequency 
for  all  events,  with  relative  amplitudes  roughly  Indicated  by  the 
symbols,  were  drawn  up  separately  for  Love  and  Rayleigh  waves.  (See 
Figure  10  and  its  caption  for  more  detail.) 

The  next  step  in  sorting  through  this  profusion  of  data  consisted 
of  assigning  an  average  path  structure  and  calculating  theoretical 
dispersion  curves  for  comparison  with  band-pass  results. 

Theoretical  Path  Dispersion  -  Choice  of  Average  Path  Structure 

The  search  for  a  model  to  generate  an  appropriate  theoretical 
dispersion  curve  representative  of  an  average  structure  along  the  path 
from  New  Brunswick  to  SCP  was  confounded  by  the  fact  that  most  of  the 
existing  studies  focused  on  geological  provinces  bordering  on  the 
desolate  New  Brunswick  region,  but  not  encompassing  it.  A  comprehensive 
summary  of  eastern  U.S.  structures  Is  presented  by  O'Neill  (1983). 

Brune  and  Dorman  (1963)  used  least-squares  inversion  of  Love  and  Ray¬ 
leigh  phase  velocities  to  determine  their  three-layered  crust  for  the 
Canadian  Shield.  Dorman  and  Ewing  (1962)  refined  the  previous  work  of 
Oliver,  Kovach,  and  Dorman  (1960),  using  Rayleigh  wave  phase  velocity 
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dispersion  to  deduce  structure  In  the  New  York-Pennsyl vania  region  that 
was  consistent  with  the  seismic  refraction  study  performed  by  Katz 
(1955).  Mitchell  and  Herrmann  (1979)  used  single-station  group  velocities 
for  earthquakes  originating  in  the  central  U.S.  with  known  fault  plane 
solutions  to  Invert  for  shear  velocity  structure  in  the  Eastern  United 
States.  Recent  work  in  New  England  by  Curtin  et.al.  (1983)  refined 
previous  compressional  and  shear  wave  velocity  structures  by  inverting 
short  period  data  from  quarry  blasts  and  local  earthquakes.  Taylor  and 
Toksoz  of  MIT  have  recently  completed  several  studies  of  northeast  United 
States  structure.  Their  1979  Inversion  of  teleseismic  P-wave  travel 
time  residuals  recorded  by  the  NE  United  States  Seismic  Network  explored 
features  of  regional  structures  whose  crustal  portion  was  single-layered. 
In  their  1982  paper,  Taylor  and  Toksoz  (1982a)  summarized  previous 
studies,  keying  velocity  structures  to  their  corresponding  regions  on  a 
geologic  map  of  the  northeastern  stretch  of  the  Appalachian  orogen.  The 
most  fruitful  report  for  the  present  study  proved  to  be  their  1982  inter¬ 
station  study  of  fundamental  mode  Rayleigh  phase  and  group  velocities  of 
teleseismic  events  measured  along  great  circle  paths  which  sampled 
different  regions  In  the  Appalachians  and  the  Grenville  Province  (Taylor 
and  Toksoz,  1982b).  Figure  9  Illustrates  the  various  paths  studied. 

Their  path  A1  sampled  structure  parallel  to  the  tectonic  fabric  of  the 
Appalachians  utilizing  vertical  component,  long-period  seismograph 
recordings  from  St.  Johns,  Newfoundland;  Weston,  Massachusetts;  and 
Ogdensburg,  New  Jersey,  of  earthquakes  originating  In  Sicily  and  Crete. 
The  central  New  England  Appalachians  (Path  A2)  were  sampled  with  events 
from  the  Kuriles  and  the  mid-Atlantic  Ridge  as  recorded  at  Montreal - 
Weston  and  Ottawa-Weston  station  pairs.  Path  G  sampled  the  eastern 
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Grenville  Province  using  events  from  South  America  and  the  Caribbean  as 
recorded  by  the  Ogdenburg-Ottawa,  and  Ogdenburg-Montreal  station  pairs. 
On  the  basis  of  the  regions  sampled,  Paths  A1  and  Path  G  were  selected 
along  with  Curtin's  model  as  the  starting  structures.  The  layer  para¬ 
meters  obtained  for  each  proposed  structure  are  given  in  Table  4. 

Using  Harkrider's  scanning  routine  (RLSCFRQ)  and  the  Rayleigh  and 
Love  eigenvalue  generating  programs  (RAYLIT.UPEC  and  LOVIT.UP),  theore¬ 
tical  fundamental  mode  dispersion  curves  were  generated  for  the  three 
trial  structures  and  superimposed  on  the  plots  of  group  velocity  vs. 
frequency  which  had  been  assembled  for  the  January  9,  1982,  mainshock. 
Agreement  between  the  theoretical  dispersion  curve  for  Path  G  and  the 
high  amplitude  group  arrivals  was  remarkable,  especially  for  the  funda¬ 
mental  mode  Love  waves  [Figures  10  (a)  and  (b)].  While  the  fit  was  not 
quite  as  good  for  the  Rayleigh  dispersion  curve,  it  seemed  to  provide 
the  best  overall  match  between  the  band-pass  data  and  the  theoretical 
curves. 

In  an  attempt  to  Improve  the  fit  of  the  theoretical  dispersion 
curves  with  the  result  of  the  band-pass  filtering,  the  Path  G  layer 
parameters  were  modified  In  keeping  with  the  expected  effects  of 
changing  crustal  thickness  and  velocity  contrast  at  the  base  of  the 
crustal  layer.  Since  separation  In  frequency  between  successive  Airy 
phases  Is  controlled  by  the  velocity  contrast  at  the  base  of  the 
crustal  layer  (Alexander,  1963),  the  two  effects  could  accommodated 
by  simply  reducing  the  shear  wave  velocity  In  layers  three  and  four 
to  3.5  km/sec,  thereby  thickening  the  uniform  layer  to  35  km  and 
Increasing  the  velocity  contrast  from  3.8  to  4.5  km/sec  to  3.5  to 
4.5  km/sec  at  the  base  of  the  layer.  Dispersion  curves  were  gener- 
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Table  4 


Taylor  and  Toksoz  Proposed  Structures 


PATH  A1 


Layer 

Thickness, 

Depth  (km) 

VP 

vs 

Density  (g/cc) 

1 

10. 

0. 

6.0 

3.4 

2.7 

2 

10. 

10. 

6.2 

3.5 

2.75 

3 

10. 

20. 

7.0 

4.0 

3.0 

4 

10. 

30. 

7.0 

4.0 

3.0 

5 

20. 

40. 

8.1 

4.6 

3.3 

6 

20. 

60. 

8.1 

4.5 

3.3 

7 

20. 

80. 

8.2 

4.5 

3.3 

8 

100. 

8.2 

4.6 

3.3 

PATH  G 

1 

5. 

0. 

6.0 

3.5 

2.7 

2 

10. 

5. 

6.0 

3.5 

2.7 

3 

10. 

15. 

6.6 

3.8 

2.8 

4 

10. 

25. 

6.8 

3.8 

2.85 

5 

20. 

35. 

8.1 

4.5 

3.3 

6 

20. 

55. 

8.1 

4.5 

3.3 

7 

20. 

75. 

8.2 

4.7 

3.3 

8 

95. 

8.2 

4.7 

3.3 

Curtin's 

Structure 

Thickness  (km) 

VP 

Vs 

Thickness  (km) 

8.8 

6.00 

3.51 

11.3 

11.3 

6.35 

3.86 

7.2 

19.8 

6.95 

4.04 

19.9 

8.18 

4.67 
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ated  for  this  modified  structure  and  plotted  in  Figures  11  (a)  and  (b) 
as  dotted  lines  in  contrast  with  the  solid  line  corresponding  to  the 
original  model.  Clearly,  this  modification  overcompensated  for  the 
previously  noted  discrepancies  in  the  Rayleigh  dispersion  fit;  more¬ 
over,  it  worsened  the  fit  for  the  Love  wave  dispersion.  As  a  result, 
the  original  Path  G  model  was  retained  as  the  preferred  structure  to 
represent  the  average  path  from  New  Brunswick  to  SCP  for  the  remainder 
of  the  study. 

While  the  fits  are  good  in  the  long  period  range  of  the  plots,  it 
was  not  clear  what  was  happening  to  the  fundamental  mode  in  the  shorter 
periods,  below  ten  seconds.  Dispersion  of  short  period  fundamental  mode 
surface  waves  is  governed  primarily  by  shallow  structure  and  a  low  veloc¬ 
ity  sedimentary  layer  is  assured  to  exist  along  at  least  some 
substantial  portion  of  the  path.  Thus,  it  is  likely  that  velocities  at 
periods  below  10  seconds  are  lower  than  the  theoretical  curves  predict  in 
Figures  10  and  11.  Evidence  of  later,  high  amplitude  arrivals  (indicated 
by  stars  in  the  band-pass  results)  support  this  interpretation. 

Higher  order  modes  were  expected  in  the  period  range  below  10  seconds. 
In  hopes  of  Identifying  these  arrivals,  theoretical  dispersion  curves  were 
generated  for  the  first  four  higher  order  modes.  Figures  12  and  13  illus¬ 
trate  the  predicted  dispersion  curves  for  the  Path  G  structure,  carried 
to  the  highest  frequency  for  which  stable  calculations  were  possible. 

Group  velocities  from  the  bandpass  results  for  the  mainshock  were  plotted 
for  comparison.  Appendix  IV  shows  these  same  dispersion  curves  super¬ 
imposed  on  the  band-pass  results  for  each  individual  event  as  well  as  on  the 
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composite  plots  of  all  the  data.  Although  individual  mode  arrivals  can't 
be  identified,  the  overall  theoretical  patterns  are  consistent  with  the 
band-pass  results. 

Theoretical  Amplitude  Spectra 

Ip  order  to  generate  theoretical  spectral  amplitudes  for  a  number 
of  modes  at  various  depths,  a  structure  at  the  source  had  to  be  specified 
and  a  source  mechanism  chosen.  Since  Harkrider's  program  to  generate 
synthetic  seismograms,  EQTR.ZBIG,  was  set  up  to  utilize  the  same  struc¬ 
ture  at  the  source  as  was  assumed  for  the  propagation  path,  the  Path  G 
structure  was  taken  to  also  represent  the  structure  at  the  source.  This 
was  not  unreasonable  in  light  of  the  P-wave  velocity  structures  present 
for  nearby  regions  in  Taylor  and  Toksoz  (1982a)  and  density  values  for 
typical  crustal  rocks  In  the  epicentral  region  provided  by  Dainty  et.al. 
(1966).  Wetmiller  et.al.  (1983)  used  a  constant  velocity  half-space  in 
their  hypocentral  location  calculations  with  P  and  S  velocities  of  6.20 
and  3.57  km/sec.  Some  limited  refraction  data  in  the  area  suggest  a 
significantly  lower  P  velocity  near  the  surface  (reported  by  Wetmiller 
et.al.  (1983)  to  G.  LeBlanc,  personal  communication),  which  places  both 
values  very  much  In  line  with  the  P  and  S  velocities  of  the  chosen 
Path  G  structure. 

A  number  of  investigators  have  proposed  source  mechanisms  for  the 
New  Brunswick  events  based  on  a  variety  of  methods  of  analysis.  It  is 
appropriate  to  sketch  here  the  major  conclusions  of  these  studies  since 
there  Is  considerable  variability  among  the  mechanisms  for  these  events. 
Including  the  malnshock. 

Well-defined  P-wave  pulse  shapes  were  analyzed  by  Choy  et.al. 
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(1983)  to  obtain  source  characteristics  for  the  January  9,  1982,  earth¬ 
quake  from  broadband  GDSN  data.  Moment  tensor  inversions  of  the  long 
period  P  and  SH  waves  provided  largely  corroborative  evidence  for  their 
choice  of  focal  parameters;  some  ambiguity  was  noted  and  attributed  to 
differences  in  the  data  sets  used  in  the  two  analyses  and  to  the  evident 
violation  of  a  simple  point  source  assumption  for  this  event.  Although 
acknowledging  as  significant  the  discrepancies  in  depth  between  their 
9  km  hypocenter  location  and  the  shallower  (0-7  km)  locally  recorded 
aftershocks  (Wetmiller  et.al.,  1982),  Choy  et.al.  (1983)  maintained  that 
their  interpretation  favored  the  greater  focal  depth.  Nabelek  et.al. 
(1982)  inverted  short  and  long  period  body  waveforms  (P  and  SH)  and 
amplitude  spectra  of  vertical  component  Rayleigh  waves  to  determine 
focal  depth,  source  orientation  and  scalar  seismic  moment.  Although 
consistent  in  type  and  orientation,  their  earthquake  parameters  fell 
outside  the  range  of  uncertainty  of  Choy’s  results  and,  in  particular, 
indicated  a  hypocentral  depth  of  7  km.  In  a  later  report  of  this  work, 
Nabelek  and  Toksoz  (1983)  reported  revised  source  parameters  more  in  line 
with  Choy’s,  but  maintained  the  shallower  depth  of  7  km. 

It  was  clear  that  the  teleseismically  recorded  aftershocks  of 
January  9  and  11  had  hypocenters  which  were  significantly  different  from 
one  another  and  from  the  mainshock  (Choy  et.al.,  1983).  As  tempting  as  it 
was  for  simplicity  to  assume  that  these  events  occurred  on  a  single  plane 
with  a  similar  mechanism,  it  was  clear  both  from  teleseismic  data  and  the 
locally  recorded  aftershocks  that  this  was  not  the  case.  If  it  could 
be  Inferred  that  the  mainshock  and  the  January  9  and  11  aftershocks 
occurred  on  a  single  plane.  It  would  require  assumption  of  a  fault 
plane  striking  HHZ  and  dipping  steeply  eastward.  This  orientation  is 
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incompatible  with  nodal  planes  inferred  from  the  broadband  data  (Choy  et. 
al.,  1983)  and  with  direct  evidence  provided  by  aftershock  patterns  and 
by  the  north-trending,  westward-dipping  bedrock  crack  found  in  May,  1982, 
at  the  southern  boundary  of  the  January  aftershock  zone  (Wetmiller 
et.al.,  1982).  Furthermore,  it  cannot  account  for  the  hypocentral  depth 
of  the  March  31  aftershock,  which  is  inferred  to  be  less  than  5  km  (from 
Wetmiller  et.al.,  1983). 

In  probably  the  most  comprehensive  analysis  of  the  New  Brunswick 
events  to  date,  Wetmiller,  et.al.  (1983)  present  body  and  surface  wave 
evidence  in  support  of  their  choice  of  focal  mechanism  for  the  mainshock 
and  a  detailed  accounting  of  the  spatial  distribution  and  composite 
mechanisms  of  the  aftershocks  as  determined  from  three  field  surveys 
conducted  in  January,  April,  and  June,  1982,  following  the  largest 
events.  Their  source  parameters  differ  from  those  of  Choy  et.al.  (1983) 
in  dip  and  rake  and  they  account  for  a  discrepancy  (in  dip,  at  least)  as 
being  due  to  Interpretation  of  a  double  event  where  movement  progressed 
up  a  steepening  fault  plane.  Aftershock  activity  convincingly  supports 
conjugate  rupturing  on  intersecting  fault  planes  of  steepening  dip 
toward  the  surface.  They  propose  that  the  January  9  and  March  31  after¬ 
shocks  occurred  on  the  eastern  limb  and  the  January  11  aftershock  on  the 
western  limb.  Figure  14  illustrates  the  described  fault  orientation. 
Hasegawa  (1983)  used  well-defined  short-period  P-wave  first  motions  and 
analysis  of  fundamental  mode  Love  and  Rayleigh  surface  waves  to  deter¬ 
mine  P-nodal  solutions  for  the  mainshock  consistent  with  those  based  on 
the  aftershock  activity  and  encompassing  the  previously  reported  results 
based  on  moment  tensor  Inversions.  He  reported  that  study  was  continuing 
on  the  comparatively  sparse  and  indistinct  P-wave  first  motions  and  sur- 
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face  wave  records  of  the  larger  aftershocks  with  the  specific  objective 
of  establishing  their  differences  in  source  mechanism  as  compared  with 
the  January  9  mainshock.  Thus,  despite  the  abundance  of  data,  there 
remain  substantial  disagreements  with  regard  to  the  fault  mechanism  and 
focal  depth  of  the  mainshock  and  all  of  the  large  aftershocks.  Focal 
depths  for  the  aftershocks  have  not  been  determined  based  on  recordings 
at  teleseismic  distances. 

The  relevant  source  parameters  for  the  individual  proposed 
mechanisms  are  summarized  in  Table  5.  Using  the  values  of  strike,  dip, 
and  rake  reported  in  these  studies,  theoretical  surface  wave  modal  exci¬ 
tations  were  calculated  assuming  the  azimuth  appropriate  for  propagation 
to  SCP.  Absolute  spectral  amplitudes  were  calculated  for  fundamental 
and  four  higher  order  modes  and  corresponding  Love-to-Rayleigh  ratios 
were  formed  for  assumed  depths  from  1-10  km.  Individual  plots  by 
depth  of  these  theoretical  spectra  for  each  of  the  source  mechanisms  in 
Table  5  are  included  in  Appendix  5;  Love-to-Rayleigh  fundamental  mode 
spectral  ratios  for  each  depth  are  shown  in  Figures  15  -  18.  The  test 
of  whether  surface  wave  data  from  a  single  station  could  corroborate 
a  preferred  mechanism  rested  in  the  degree  of  fit  between  the 
theoretical  spectra  and  the  results  of  the  band-pass  filtering. 


Mainshock 


Strike 

Dip  (NW) 

Rake 

Depth 

Choy  et. al . 

195* 

(318*)* 

65* 

70* 

9.0  +  1.3  km 
(from  depth  phases) 

Nabelek  et.al. 

3*  ( 183  * ) 
(306*)* 

34* 

95* 

7  km 

Nabelek  and  Toksoz 
(in  preparation) 

175  +  5* 
(29D* 

55  +  2* 

87  +  10* 

7  +  1  km 

Wetmiller  et.al . 

195* 

(318*)* 

50* 

120* 

7  km 

(from  aftershock 
and  SW  analysis) 

January  11  Aftershock 

SE 

Wetmiller  et.al . 

332* 

48* 

59* 

6  -  7  km 

(95*)* 


*  These  are  the  "azimuths  from  the  hanging  wall,"  as  required  by 
Harkrider's  convention.  See  Appendix  1  for  details. 
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CHAPTER  5 

INTERPRETATION  OF  RESULTS 

Comparison  of  theoretical  dispersion  curves  and  empirical  band¬ 
pass  filter  results  (Figures  12  and  13  and  Appendix  4)  permitted  ready 
Identification  of  the  long  period,  fundamental  mode  Rayleigh  and  Love 
waves.  Thus,  conclusions  regarding  source  mechanism  and  depth  were 
based  entirely  on  their  Love-to-Rayleigh  amplitude  ratios.  Although  the 
presence  of  higher  modes  was  supported  by  the  data,  accurate  identifi¬ 
cation  of  individual  modes  was  out  of  the  question,  eliminating  these 
as  sources  of  independent  evidence  for  detailed  analysis  of  the  differ¬ 
ences  among  events.  Some  interesting  observations  were  gleaned 
from  the  high  frequency,  Lg  arrivals,  however:  some  consistent  with 
higher  mode  excitation,  some  providing  useful  m^g  estimates.  These 
are  presented  as  empirical  evidence,  worthy  of  further  attention  in 
subsequent  investigations. 

Fundamental  Mode  Analysis  for  Depth  and  Source  Mechanism 

In  comparing  theoretical  Love  wave  spectra  (Appendix  5,  Figures  A5-la 
through  20a)  for  different  mechanisms,  it  becomes  Imnediately  apparent  that 
the  excitations  are  almost  Identical  for  the  fundamental  mode.  The 
Nabelek  and  Choy  mechanisms  (Table  5)  show  nearly  identical  features  for 
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Nabelek  and  Choy  mechanisms  (Table  5)  show  nearly  identical  features  for 
fundamental  mode  Rayleigh  waves;  where  differences  occur,  they  reflect 
similar  shape  but  higher  amplitude  trends  for  the  Nabelek  mechanism 
(compare  Figures  A5-lb  through  5b  with  A5-6b  through  10b).  As  depth  in¬ 
creases,  the  discrepancies  in  amplitude  decrease  and  are  more  restricted 
to  increasingly  longer  periods.  Wetmiller  et.al.'s  January  9  mechanism 
yields  Rayleigh  wave  amplitude  spectra  which  differ  substantially  from 
the  Choy  spectra  for  shallow  source  depths  (compare  Figures  A5-lb  through 
5b  with  A5  -lib  through  15b).  Remarkable  agreement  is  evident  in  the 
short  period  range  of  the  spectrum  as  depth  increases,  while  very 
different  amplitudes  are  seen  for  periods  greater  than  10  seconds  for 
all  depths.  Wetmiller  et.al.'s  January  11  mechanism  is  substantially 
different  in  Rayleigh  excitation  at  all  depths  (See  Figures  A5-16b 
through  20b),  with  a  prominent  spectral  hole  which  moves  to  longer 
periods  with  increasing  focal  depth. 

From  these  observations,  it  was  clear  that  the  spectral  differences 
that  were  most  likely  to  be  diagnostic  of  source  mechanism  and  depth 
were  the  Love-to-Rayleigh  ratios  for  periods  greater  than  10  seconds. 

This  was  especially  appealing  in  light  of  the  fact  that  long  period 
Rayleigh  and  Love  fundamental  modes  can  be  unambiguously  identified 
and  that  all  band-pass  filters  In  this  period  range  were  the  same. 
Theoretical  Love-to-Rayleigh  spectral  ratios  were  then  formed  for 
fundamental  mode  periods  greater  than  ten  seconds,  superimposing 
ratios  for  all  depths  on  a  single  plot  for  each  mechanism.  These 
results  are  displayed  in  Figures  15  to  18.  The  differences  In  the 
values  of  Love-to-Rayleigh  ratios  between  the  Nabelek  and  Choy 
mechanisms  (Figures  15  and  16)  are  Immediately  apparent  when  viewed  in 
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this  way,  though  peak  amplitudes  with  depth  occur  at  almost  identical 
frequencies  for  the  two  mechanisms.  Wetmiller  et.al's  January  9  and 
January  11  proposed  mechanisms  exhibit  distinctive  features  in  their 
spectral  ratios  (see  Figures  17  and  18). 

Love  and  Rayleigh  fundamental  mode  arrivals  were  selected  from  the 
band-pass  filter  results  (on  transverse  and  vertical  components,  respec¬ 
tively).  Love-to-Rayleigh  amplitude  ratios  were  computed  and  plotted, 
by  event,  on  the  same  amplitude  vs.  period  scales  as  used  for  the 
theoretical  spectral  ratios.  The  mainshock  amplitude  ratios  exhibit 
a  smooth  relationship  in  period,  most  consistent  with  Choy's  mechanism. 
Superimposing  the  data  points  on  the  theoretical  curves  shows  spectral 
ratios  consistent  with  depths  ranging  from  5-8  km,  tending  toward 
6  km  at  periods  greater  than  20  seconds  and  toward  8  km  for  shorter 
periods  (Figure  19).  Equally  convincing  are  the  results  for  the 
January  11  aftershock  (Figure  20)  which  show  evidence  of  a  distinct 
peak  in  the  spectral  ratio  entirely  consistent  with  Wetmiller  et.al.'s 
January  11  mechanism,  with  all  the  data  points  lying  very  close  to 
the  curve  corresponding  to  6  km  depth.  Again,  there  is  a  tendency 
for  the  short  period  data  points  to  indicate  a  slightly  greater 
depth. 

The  results  for  the  January  9  aftershock  (Figure  21),  assuming  the 
same  mechanism  as  the  mainshock,  showed  a  great  deal  of  scatter  but,  in 
general,  had  higher  magnitude  Love-to-Rayleigh  ratios  than  those  for  the 
mainshock  at  periods  less  than  20  seconds,  with  an  Indication  of  a  peak 
ratio  located  between  15  and  20  seconds  period,  that  is,  at  a  lower 
period  than  that  of  the  mainshock.  This  is  consistent  with  trends 
associated  with  shallower  depths  on  the  theoretical  Love-to-Rayleigh 
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curves;  a  shallower  depth  than  the  mainshock  is  thus  suggested  for  this 
event.  However,  the  nunerical  values  of  the  ratio  are  inconsistent  with 
either  Choy’s  or  Nabelek's  mechanisms. 

In  an  attempt  to  modify  source  parameters  to  improve  the  fit  for 
this  event,  it  might  seem  that  Love-to-Rayl eigh  ratios  would  be  more 
sensitive  to  changes  in  strike  than  either  dip  or  rake.  However,  evalua¬ 
tion  of  the  rate  of  change  of  the  Love  and  Rayleigh  radiation  pattern 
functions  with  each  parameter,  over  ranges  encompassing  all  the  proposed 
mechanisms,  indicated  that  the  ratio  is  in  fact  most  sensitive  to  dip. 
Since  the  strike  of  the  fault  plane  seems  so  much  better  constrained  for 
this  event  than  either  dip  or  rake,  it  made  sense  to  vary  these  in  seek¬ 
ing  a  better  match  between  the  data  and  the  synthetic  spectra. 

Taking  into  account  the  differences  in  dip  and  rake  among  the  pro¬ 
posed  mechanisms,  additional  synthetic  spectra  were  generated  for  varia¬ 
tions  of  dip  and  rake  on  Choy's  original  model:  rake  values  of  70*.  90* 
and  110*  were  assumed  with  dip  values  of  65*  and  55’  respectively.  The 
effects  of  these  changes  in  dip  and  rake  can  be  seen  in  Figure  22.  As 
rake  is  increased,  the  curves  tend  to  flatten,  with  peak  ratio  values 
occurring  at  longer  periods,  both  of  which  trends  are  inconsistent  with 
the  data.  As  dip  is  decreased,  the  ratios  take  on  much  higher  values, 
peaking  at  lower  periods  than  those  for  corresponding  rakes  and  the 
larger  dip  value.  The  same  trend  to  peak  at  somewhat  longer  periods 
with  Increasing  rake  Is  observed  with  the  shallower  dip.  Of  particular 
Interest  Is  comparison  of  the  Choy  mechanism  (  9  *  318* »  b  ■  65*, 

A  *70*)  (Figures  15  and  19),  the  modified  Choy  mechanism  (  9  *  318*, 
b  *  55*,  A*  90*)  (Figure  22b)  and  the  Nabelek  mechanism  (  9  •  306*, 
i  *34*,  A  *95*)  (Figure  16).  The  three  sets  of  curves  show  similar 
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trends,  though  the  amplitude  ratio  dependence  on  decreasing  dip  and  in¬ 
creasing  rake  is  not  apparent.  The  closest  match  of  the  data  for  the 
January  9  aftershock  corresponds  to  the  modified  Chcy  mechanism  -  dip 
of  55*  and  rake  of  90’  (Figure  23);  this  result  suggests  that  a  shallow 
depth  on  the  order  of  3  km  might  be  more  consistent  with  the  overall  fit, 
though  scatter  in  the  data  precludes  identification  of  the  period  asso¬ 
ciated  with  the  peak  in  the  Love-to-Rayleigh  ratios. 

Considering  the  dip  and  rake  trends  indicated  by  these  tests, 
increasing  the  dip  and  decreasing  the  strike  beyond  the  ranges  tested 
might  produce  a  similar  effect,  though  this  is  only  conjecture.  It  is, 
however,  the  only  way  to  proceed  if  the  fault  plane  is  indeed  dipping 
more  steeply  near  the  surface  as  Wetmiller  et.al.  (1983)  propose.  The 
data  seem  to  Indicate  a  shallower  depth;  beyond  that,  little  surface 
wave  evidence  is  available  to  constrain  the  mechanism  for  this  event. 

Long  period  data  for  the  June  16  aftershock  is  sparse  due  to  many 
problems  inherent  in  the  data.  While  any  interpretation  suffers  because 
of  these  difficulties,  the  fact  that  the  observed  broad-band  Love-to- 
Rayleigh  ratios  are  very  high  prompted  a  comparison  of  the  shorter  period 
band-pass  results  with  the  theoretical  spectra  for  Wetmiller  et.al.'s 
January  11  mechanism  (Figure  24).  Apart  from  the  two  spurious  points  in 
the  long  period  range,  the  data  seem  to  point  to  a  deeper  event,  on  the 
order  of  8  km  depth.  No  interpretation  beyond  this  was  attempted  for 
this  event. 

None  of  the  fundamental  mode  data  show  patterns  typical  of 
Wetmiller' s  mechanism  for  the  mainshock. 

Both  the  theoretical  and  observational  results  discussed  here  show 
that  the  fundamental  mode  surface  waves  in  the  period  range  of  10  to 
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50  seconds  can  be  used  to  Infer  differences  in  source  mechanism  for 
shallow  (less  than  10  km)  crustal  earthquakes.  In  particular,  depth 
differences  on  the  order  of  one  kilometer  appear  to  be  distinguishable 
with  high-quality  fundamental  mode  data.  This  may  be  significant  for 
improving  depth  determinations  from  regional  network  recordings  where 
accurate  epicenter  location  and  focal  mechanism  determination  typically 
can  be  made,  but  source  depth  is  uncertain. 

Other  Evidence  for  Focal  Depth 

The  usual  problems  inherent  in  using  body  wave  observations  to 
determine  focal  depth  for  shallow  earthquakes  were  encountered  here. 
Examination  of  the  GDSN  short  period  waveforms  for  these  four  events  re¬ 
vealed  few  easily  Identified  pP  arrivals.  Probably  the  clearest  indica¬ 
tion  of  a  useful  depth  phase  appears  on  the  short  period  vertical  signal 
for  the  mainshock  as  recorded  at  station  GRFO  (Germany)  (see  Figure  25a). 
Assuming  Choy's  mechanism  as  the  preferred  focal  mechanism  for  this 
event,  (Figure  26),  the  P  arrival  at  GRFO  (  A  *  47*  in  a  northeasterly 
direction  from  the  source,  with  take-off  angle  of  26.5*)  plots  well 
within  the  compresslonal  quadrant,  while  pP  falls  near  a  node.  Thus, 
the  distinct  arrival  following  the  emergent  P  is  identified  as  sP. 
Measurement  of  the  P-sP  delay  time  and  assuning  P  and  S  velocities  con¬ 
sistent  with  assuned  source  structure  yield  s  depth  for  this  event  of 
9  km.  This  Is  consistent  with  Choy  et.al.’s  (1983)  analysis  of  the  sP 
arrival  at  TOl  (Spain),  for  which  they  also  calculate  a  focal  depth  of 
9  km  for  the  mainshock.  Two  other  stations,  BOCO  (Colombia)  and  ZOBO 
(Bolivia)  (A»  45*  and  65*,  respectively,  just  west  of  due  south  of  the 
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source)  would  be  expected  to  show  strong  pP  arrivals  for  the  mainshock 
and  all  the  aftershocks;  while  very  clearly  present  in  the  waveform, 
a  distinct  onset  of  a  pP  could  not  be  identified  for  the  mainshock. 

The  January  9  aftershock  showed  an  excellent  depth  phase  at  BOCO 
(Figure  25b);  this  is  Identified  as  p P  since  the  mechanism  is 
generally  the  same  as  the  mainshock.  The  measured  P-pP  delay  time  of 
2.7  seconds  translated  to  a  depth  of  7.5  km  for  this  event.  The  signal 
recorded  at  ZOBO  for  the  January  11  aftershock  (Figure  25c)  illustrated 
the  more  typical  arrival  observed  for  all  the  events;  further  phase 
identification  and  depth  calculations  were  not  attempted  in  this 
study.  Choy  et.al.  (1983)  inferred  a  pP  arrival  and  a  corresponding 
focal  depth  of  6  km  for  the  January  11  aftershock  from  the  ZOBO  velocity 
record,  assuming  that  the  mechanism  for  the  event  was  not  significantly 
different  from  that  of  the  mainshock.  Depths  determined  from  these 
teleseismic  depth  phases  give  focal  depths,  in  most  cases,  at  least  2  km 
deeper  than  those  Inferred  from  surface  wave  analysis  and  local  after¬ 
shock  observations.  Only  the  aftershock  on  January  11  gives  exact 
agreement  on  the  result. 

Additional  depth  information  was  provided  by  the  near-source  moni¬ 
toring  of  aftershock  activity  as  reported  in  Wetmiller  et.al .(1983) . 
Figure  27  provides  a  multi-dimensional  view  of  epicentral  and  hypocentral 
locations  of  the  aftershocks,  revealing  that  this  activity  is  restricted 
to  depths  less  than  7.5  km.  These  data  suggest  shallower  depths  than 
those  indicated  by  body  wave  analyses,  however,  since  the  monitoring 
system  was  not  Installed  until  January  10,  much  important  Information 
on  the  aftershock  activity  that  Immediately  followed  the  mainshock  and 
first  large  aftershock  is  missing.  This  evidence,  nonetheless.  Indicates 
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that  the  malnshock  rupture  began  at  depth  and  propagated  upward. 

Summary  of  Results  on  Depth  Determination 

The  depths  for  the  mainshock  and  January  9  aftershock  obtained  in 
this  study  using  fundamental  mode  surface  waves  appear  to  be  somewhat 
shallower  than  depths  inferred  from  body  wave  depth  phases  observed  at 
teleseismic  distances  and  are  more  consistent  with  the  depths  suggested 
by  the  locally  monitored  aftershock  activity.  Assuming  that  all  three 
types  of  observations  are  reliable  in  determining  depth,  and  recalling 
the  shape  of  the  curve  of  the  Love-to-Rayleigh  ratio  for  the  mainshock 
data,  the  apparent  inconsistencies  might  be  accounted  for  by  assuming 
the  upward  propagating  fault  rupture  suggested  by  Choy  et.al.  (1983). 

The  high  frequency  body  wave  first  arrivals  would  reflect  the  point  of 
initiation  of  the  rupture  while  the  longer  period  body  and  surface  wave 
arrivals  would  reflect  an  average  depth  over  which  the  fault  moved. 

This  is  also  consistent  with  the  trend  in  Love-to-Rayleigh  ratios  to 
give  greater  depths  at  shorter  periods  observed  in  the  band-pass  ratios 
for  the  mainshock  and  January  11  aftershock.  Furthermore,  it  logically 
follows  that  the  aftershocks  would  be  largely  restricted  to  depths 
shallower  than  the  initial  rupture  if  it  indeed  propagated  upward  toward 
the  surface.  To  test  the  idea  of  an  upward  propagating  rupture  for  the 
mainshock,  theoretical  spectral  amplitudes  for  different  combinations  of 
depths  were  summed  separately  for  Love  and  Rayleigh  contributions  and  the 
ratios  calculated  and  plotted  for  comparison  with  the  same  data  analyzed 
In  Figure  19.  Results  of  preliminary  testing  are  presented  in  Figure  28. 
The  trend  for  Love-to-Rayleigh  ratios  to  peak  at  shorter  periods  for 
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shallower  depth  ranges  is  clear  in  Figure  28a;  the  higher  ratio  values 
associated  with  deeper  events  is  illustrated  in  Figure  28b.  After  a 
number  of  attempts  at  fitting  the  data  to  a  combination  of  various  depth 
contributions,  the  best  match  resulted  from  summing  the  spectra  for 
depths  of  4,  5,  and  8  km,  as  shown  in  Figure  29.  With  the  initial  rup¬ 
ture  set  at  8  km,  Choy’s  body  wave  depth  determination  of  9  km  is  more 
acceptable  in  light  of  this  surface  wave  evidence.  Focal  depth  for  the 
January  11  aftershock  is  consistently  determined  at  6  km  by  all  three 
methods;  while  the  evidence  is  scanty  for  the  June  16  event,  a  depth 
on  the  order  of  7  -  8  km  is  Indicated  by  both  surface  wave  and  local 
aftershock  results.  Table  6  provides  a  concise  summary  of  depths 
determined  by  the  various  methods. 


Table  6 

Summary  of  Depths  Determined  by  Various  Methods 


Event 

Body  Wave 
Evidence 

Aftershock 

Evidence 

Surface  Wave  Evidence 
of  Present  Study 

1/9/82 

mainshock 

9  km 

7  km 

7  _+  1.5  km 

1/9/82 

aftershock 

7.5  km 

7  km 

3  -  5  +  2  km 

1/11/82 

aftershock 

6  km 

6  -  7  km 

6  +  1  km 

6/16/82 


7  km 


8  +  2  km 
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lg  Arrivals  and  m^g  Estimates 


Since  higher  order  modes  were  not  resolvable  In  the  data,  comparison 
with  theoretical  spectra  was  unfounded.  However,  three-dimensional  plots 
of  peak  amplitudes  vs.  frequency  in  the  higher  frequency  range  associated 
with  Lg  arrivals  did  indicate  some  consistency  with  values  predicted 
theoretically.  This  evidence  is  discussed  in  greater  detail  in  Appendix  6. 

Some  additional  observations  were  made  in  the  short  period  region  of 
the  spectra  by  considering  the  relative  peak  amplitudes  of  the  original 
waveforms  on  the  short  and  intermediate  period  vertical  signals  and  the 
band-pass  traces  and  relating  these  to  the  published  magnitudes  of  the 
events.  Table  7  is  a  compilation  of  the  available  magnitude  data  for 
these  events. 


Table  7 


Month/Day/Event 

mjj 

mbLg* 

Ms* 

_  ** 
mbLg 

u  ★  ★★ 
Ms 

1/9  mainshock 

5.7 

5.8 

5.2 

5.9 

4.9 

1/9  aftershock 

5.1 

5.3 

3.9 

5.2 

3.6 

1/11  aftershock 

5.4 

5.5 

4.5 

5.5 

4.25 

6/16  aftershock 

4.7 

4.6 

- 

4.3 

- 

*  Reported  In  Wetmiller  et.al.  (1983).  m^g  values  are  from  the 
Canadian  Network;  Ms  are  from  NE IS 

**  As  measured  from  the  digital  records  at  SCP.  m^g  for  the 

clipped  mainshock  signal  was  calculated  from  a  short  period 
record  simulated  from  the  Intermediate  period  record 

***  As  measured  from  the  photographic  records  at  SCP 
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Nutt 1 i ‘ s  formula  (taken  from  von  Seggern  and  Alexander,  1982)  for  deter¬ 
mining  Lg  magnitude  very  specifically  requires  measurements  of  sustained, 
high  amplitude,  vertical  component  arrivals  of  one  second  period;  these 
tend  to  be  the  highest  amplitude  arrivals  on  the  short  and  Intermediate 
period  signals. 


flifcLg  *  .30  +  1.66[log(  A*)]  +  log(A/T) 

Since  the  distance,  A  ,  is  the  same  for  these  events,  the  change 
in  mfcLg  between  events  is  due  entirely  to  amplitude  differences  of  the 
one-second  arrivals.  Specifically, 

A  mbLg  *  (mbLg)2  ‘  (mbLg)l 
A  "bLg  =  log  (A2/T)  -  log  (A]/T) 

A  mbLg  *  log  (A2/A!) 

Table  8  shows  mb  values  predicted  from  a  reference  mbLg  (the 
January  11  aftershock  in  the  case  of  short  period  data  and  the  January  9 
mainshock  in  the  case  of  the  intermediate  period  data).  Evidently,  the 
ensemble  of  higher  mode,  Lg  arrivals,  propagating  along  this  path  in 
the  northeastern  United  States,  effectively  transmits  this  information 
on  source  strength.  If  Lg  can  be  shown  to  be  the  reliable  measure  of 
magnitude  that  Is  suggested,  the  results  in  Table  9  indicate  that  the 
magnitude  for  the  aftershock  on  June  16  should  be  revised  downward  to  a 
value  closer  to  4.3. 

It  Is  also  Interesting  to  notice  the  trends  in  amplitude  ratios  ob¬ 
tained  from  the  peak  arrivals  at  each  central  frequency  on  the  band-pass 
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Table  8 


Short  Period  (All  events  compared  to  1/11  aftershock  because  of 
clipping  for  the  mainshock): 

Ratio  of  Peak  Associated  Adjusted  Reported 

Signal  Amplitudes  =  log  (Ratio)  by  ^m^g  m^ 

Z  -  Component 


1/11  to  1/9 
aftershocks 

1.9 

.29 

5.1 

5.1 

1/11  to  1/11 

1.0 

0.00 

5.4 

5.4 

1/11  to  6/16 
aftershocks 

10.5 

1.02 

4.4 

4.7 

Intermedi ate 

Period 

_( AT  1  events  compared  with 

1/9  mainshock) 

1/9  main  to 

1/9  aftershock 

8.2 

.9 

4.8 

5.1 

1/9  main  to 

1/11  aftershock 

1.9 

.3 

5.4 

5.4 

1/9  main  to 

6/16  aftershock 

32.5 

1.5 

4.2 

4.7 

filters  where  all  ratios  are  comparisons  with  the  associated  mainshock 
peak  (see  Table  9).  The  expected  differences  due  to  variations  in  excita¬ 
tion  at  the  source  caused  by  the  different  source  mechanisms  are 
apparent,  but  It  is  also  of  Interest  to  note  the  frequencies  at  which 
the  ratios  are  consistent  with  the  differences  in  the  magnitudes  of  the 
events.  Except  for  the  January  11  aftershock,  that  information  is 
clearly  present  in  the  frequencies  around  1  Hz. 

In  spite  of  the  noted  differences,  trfiich  can  no  doubt  be  attributed 
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to  variations  in  modal  excitation  at  the  source,  visual  inspection  of 
the  higher  mode  contributions  to  the  theoretical  Rayleigh  amplitude 
spectra  (Figures  A5-lb  through  20b)  show  a  consistent  maximun  level  of 
excitation  for  the  modes,  taken  together,  regardless  of  depth  in  the  "gra¬ 
nitic  layer"  for  the  different  mechanisms.  This  conforms  to  the  expecta¬ 
tion  that,  for  a  given  mechanism,  the  same  fraction  of  the  total  source 
energy  remains  trapped  as  surface  wave  energy  in  the  same  constant  velocity 
layer.  Therefore,  magnitudes  determined  by  peak  amplitudes  rf  the  re¬ 
corded  signal,  representing  the  ensemble  of  Lg  modes,  should  give 
accurate  measures  of  the  strength  of  the  source  for  shallow  earthquakes. 
m^Lg  estimates  from  network  averages  and  from  the  single-station  SCP 
agree  closely  with  one  another  and  with  teleseismic  magnitude  estimates. 


53 


CHAPTER  6 

CONCLUSIONS  AND  DIRECTIONS  FOR  FURTHER  STUDY 

Single-station  surface  wave  data  from  the  1982  New  Brunswick 
sequence  of  earthquakes  were  used  to  help  determine  the  preferred 
focal  mechanisms  including  focal  depths  for  the  four  events  studied.  By 
applying  band-pass  filters  over  the  complete  frequency  range  provided  by 
the  DWWSSN  short.  Intermediate,  and  long  period  digital  data  recorded  at 
SCP,  group  velocities  and  amplitudes  were  empirically  determined  for 
fundamental  and  higher  mode  arrivals.  Theoretical  Rayleigh  and  Love  dis¬ 
persion  curves  were  generated  for  an  assumed  path  structure  which  ade¬ 
quately  fit  the  observed  data.  Synthetic  Love  and  Rayleigh  amplitude 
spectra  and  spectral  ratios  were  calculated  for  various  source  mechanisms 
proposed  by  other  investigators  and  compared  to  the  data  derived  from  the 
band-pass  filters.  Evaluation  of  source  mechanism  and  depth  determina¬ 
tion  were  based  on  fundamental  mode  data  of  periods  10  seconds  and  greater 
and  supported  by  qualitative  behavior  of  the  higher  mode  excitation. 

Fundamental  mode  Love  and  Rayleigh  ratios  corroborated  Choy  et.al.'s 
(1983)  choice  of  source  parameters  for  the  mainshock:  fault  plane 
striking  at  195*,  dip  of  65*  to  the  northwest,  and  a  rake  of  70*.  The 
preferred  average  depth  was  found  to  be  7  km  as  opposed  to  the  9  km  value 
calculated  on  the  basis  of  depth  phases  observed  at  teleseismic  distances. 
The  surface  wave  results  for  this  event  can  also  be  Interpreted  as  indi¬ 
cative  of  rupture  which  propagates  updip.  For  the  January  9  aftershock, 
a  mechanism  similar  to  that  of  the  mainshock  is  supported  and  a  shallower 
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depth  is  indicated  by  the  shift  in  the  peak  Love-to-Rayleigh  ratio  to 
shorter  periods.  The  January  11  aftershock  is  well  described  by 
Wetmiller  et.al.'s  (1983)  mechanism  which  indicates  thrusting  on  a  plane 
conjugate  to  the  fault  plane  of  the  mainshock,  with  strike  of  332*,  dip 
48*  to  the  southeast,  and  a  rake  of  59*.  Surface  wave  Love-to-Rayleigh 
ratios  and  body  wave  depth  phase  identification  both  support  a  6  km  depth 
for  this  event.  Although  subject  to  uncertainty  because  of  significant 
long  period  noise  and  delayed  triggering  of  the  intermediate  period  re¬ 
cording  system,  analysis  of  the  June  16  event  reveals  similarities  in 
mechanism  to  the  January  11  aftershock,  with  surface  wave  evidence 
supporting  a  deeper  focus,  on  the  order  of  8  km  depth.  The  chosen 
mechanism  for  the  mainshock  and  the  first  aftershock  is  not  entirely 
consistent  with  the  prevailing  ENE  maximum  compressive  stress  charac¬ 
teristic  of  most  areas  in  eastern  North  America,  although  the  inferred 
mechanisms  for  the  January  11  and  June  16  events  are. 

Conclusive  Identification  of  individual  higher  order  Love  and 
Rayleigh  modes  in  the  short  period  range  that  encompass  the  Lg  arri¬ 
vals  was  not  possible.  However,  it  was  apparent  from  comparisons  of 
the  theoretical  spectra  and  three-dimensional  views  of  the  band-pass 
results  that  differences  in  source  excitation  were  well  preserved 
for  these  waves  even  after  traveling  the  1130  km  path  to  SCP.  This  is 
further  support  of  the  observation  that  crustal  structure  In  the 
eastern  United  States  provides  a  more  uniform  wave  guide  for  Lg 
propagation  than  It  does  In  the  western  United  States  (von  Seggern  and 
Alexander,  1982),  where  source  information  (other  than  magnitude)  is 
unrecognizable  after  only  a  few  hundreds  of  kilometers  of  propagation 
from  the  source.  Furthermore,  mbLg  calculations  based  on  the  ensemble 
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of  Lg  modes  constituting  the  highest  amplitude  arrivals  on  the  short 
period  and  intermediate  period  records  proved  remarkably  consistent  with 
mb  as  a  measure  of  source  strength.  This  proved  to  be  the  case  even 
though  the  amount  of  excitation  of  individual  modes  varied  from  frequency 
to  frequency. 

Constraint  of  source  parameters  with  surface  wave  observations  can 
be  optimized  by  using  data  from  many  stations  with  good  azimuthal 
coverage  of  the  event.  Turnbull  (1976)  was  able  to  demonstrate  that  only 
six  stations  were  needed  to  fit  observed  Love  and  Rayleigh  amplitude 
spectra  via  an  Iterative,  least-squares  regression  to  well -constrained 
values  of  dip,  slip,  strike  of  the  fault  plane,  and  moment.  The  present 
study  has  shown  that,  in  the  case  where  good  body  wave  evidence  gives  a 
reasonably  accurate  starting  mechanism  and  where  path  and  source  struc¬ 
ture  can  also  be  assumed  with  confidence,  surface  wave  analysis  at  a 
single  station  can  provide  accurate  estimates  of  focal  depth  and  further 
constraints  on  mechanism,  particularly  in  the  case  of  shallow  earth¬ 
quakes.  Resolution  of  the  individual  higher  modes  is  still  highly 
desirable,  not  only  for  constraining  the  source  mechanism,  but  also  for 
inverting  for  path  structure,  one  of  the  more  interesting  possible  exten¬ 
sions  of  this  study.  How  this  can  best  be  accomplished  is  not  immediate¬ 
ly  apparent;  however,  the  theoretical  dispersion  curves  in  Appendix  3 
suggest  that  one  possible  approach  is  to  band-pass  filter  windows  of  the 
original  time  signal  that  are  limited  by  the  cut-off  frequencies  of 
successive  modes.  The  short  period  records  for  the  mainshock  and  the 
January  11  aftershock  show  clear  evidence  for  the  presence  of  a  dispers¬ 
ing  higher  order  mode,  almost  certainly  the  first  higher  mode.  Although 
this  method  will  encounter  difficulties  in  the  frequency  ranges  where 
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the  different  dispersion  curves  overlap,  with  many  modes  arriving  at 
the  sane  time,  successful  Identification  of  just  one  higher  order  mode 
would  be  useful  for  both  source  and  structure  studies. 
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FIGURE  CAPTIONS 


Figure  1  (a)  Earthquake  locations  in  northeast  United  States  and  con 
terminous  Canada,  1534-1959.  Epicentral  region  and  SCP 
location  are  indicated  with  stars  (from  Yang  and  Aggar- 
wal ,  1981;  after  Smith,  1962,  1966). 


1  (b)  Epicenters  of  earthquakes  with  >,  2  in  northeast 
United  States  and  adjacent  Canada  as  recorded  by  the 
Northeastern  Seismic  Network  for  the  period  1970-1979 
(after  Yang  and  Aggarwal,  1981). 

Figure  2  Seismicity  of  the  New  Brunswick  study  area  prior  to  the 
1982  sequence.  No  particular  trend  to  the  activity  is 
indicated,  although  two  events  with  magnitude  >  3  did 
occur  within  the  period  1970-1981.  Uncertainty  in  epi¬ 
center  location  is  not  less  than  10  km.  Note  the  starred 
locations  of  the  January  and  March,  1982,  sequence  as  well 
as  the  separate  June  16,  1982,  event  (after  Wetmiller  et. 
al.,  1983). 


Figure  3 


Figure  4 


Figure  5 


Figure  6 


Figure  7 


Geology  of  the  New  Brunswick  study  area.  Note  that  all 
the  events  of  the  earthquake  sequence  occur  within  a 
granitic  Intrusion  of  the  Miramichi  Anticlinoriun  which 
has  not  yet  revealed  any  evidence  of  previous  deformation. 
Existing,  mapped  faults  were  apparently  not  Involved  in 
the  1982  earthquake  activity  (after  Wetmiller  et.al.,  1983). 

Aftershock  activity  of  the  epicentral  region  superimposed 
on  local  geology.  Although  not  shown  on  the  map,  the 
epicenter  of  the  June  16  event  is  located  within  this 
same  granitic  pluton  (after  Wetmiller  et.al.,  1983). 

Effect  of  bandwidth  on  peak  amplitude  and  arrival  time. 

The  differences  in  group  velocity  determinations  due  to 
the  factor  of  ten  bandwidth  difference  amount  to  only 
.1  km/sec.  The  effects  on  peak  amplitudes  are  substantial; 
however,  working  with  spectral  ratios  will  cancel  these 
amplitude  effects. 

Observed  signal  and  suite  of  band-pass  filters  for  central 
frequencies  ranging  from  5  Hz  to  8.5  Hz  applied  to  the 
short  period  data  for  the  January  11  aftershock. 

Observed  signal  and  suite  of  band-pass  filters  for  central 
frequencies  ranging  from  5  Hz  to  8.5  Hz  applied  to  the 
short  period  data  for  the  June  16  event. 
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Figure  8 


Figure  9 


Figure  10 


Figure  11 


Observed  signal  and  suite  of  filters  for  central  frequen¬ 
cies  ranging  from  .05  Hz  to  .25  Hz  applied  to  the  inter¬ 
mediate  period,  vertical  component  data  for  the  January  9 
main shock. 

Great  circle  paths  and  stations  used  to  determine  two  of 
the  velocity  models  selected  as  possible  path  and  source 
structure.  Path  A1  extends  from  OGD  through  WES  to  STJ 
(Ogdensburg,  N.J.;  Weston,  Mass.;  St.  Johns,  Newfoundland). 
Path  G  represents  structure  between  OGD  and  MNT  (Montreal) 
and  OGD  and  OTT  (Ottawa)  and  is  the  model  selected  to  re¬ 
present  path  and  source  structure  in  this  study.  Approxi¬ 
mate  epicenter  and  station  locations  are  indicated  (after 
Taylor  and  Toksoz,  1982b). 

Fundamental  mode  Love  (a)  and  Rayleigh  (b)  group  velocity 
dispersion  curves  for  the  three  structure  models  tested, 
plotted  with  data  points  representing  group  velociities 
associated  with  each  peak  on  the  band-pass  filters.  Of 
the  three.  Path  G  dispersion  fits  the  data  points  best.  On 
this  and  all  subsequent  dispersion  plots,  the  symbols 
used  to  represent  the  peaks  were  assigned  by  amplitude. 

The  amplitudes  of  each  enveloped  peak  were  normalized 
to  the  peak  amplitude  for  that  band-pass  and  then 
assigned  a  symbol  according  to  the  following  scheme: 

10  -  the  point  having  the  peak  amplitude  for  the 
given  central  frequency 

9  -  those  point(s)  representing  peaks  along  the  fil¬ 
tered  signal  with  amplitude(s)  80-99%  of  maximum 
for  that  trace 

8  -  those  point(s)  representing  peaks  along  the  fil¬ 
tered  signal  with  amplitude(s)  60-79%  of  maximum 
for  that  trace 

7  -  those  point(s)  representing  peaks  along  the  fil¬ 
tered  signal  with  amplitudes  less  than  59%  of 
maximum  for  that  trace. 

Fundamental  mode  Love  (a)  and  Rayleigh  (b)  group  velocity 
dispersion  curves  for  Path  G  and  modified  Path  G  struc¬ 
tures,  plotted  with  data  points  representing  group  veloc¬ 
ities  associated  with  each  peak  on  the  band-pass  filters. 
Modifying  the  Path  G  model  was  accomplished  by  changing 
the  3.8  km/sec  shear  velocity  layers  to  3.5  km/sec 
In  order  to  effectively  thicken  the  granitic  layer 
and  sharpen  the  velocity  contrast  at  Its  base.  The 
expected  effects  on  the  theoretical  dispersion  are  evi¬ 
dent;  however,  the  original  Path  G  structure  provides 
a  better  overall  fit  to  the  data,  so  it  was  retained  for 
the  duration  of  the  study. 


64 


Figure  12 


Figure  13 


Figure  14 


Figure  15 


Figure  16 


Figure  17 


Figure  18 


Figure  19 


Theoretical  dispersion  predicted  by  the  Path  G  velocity 
model  for  fundamental  and  first  four  higher  mode  Love 
waves  for  the  period  range  .2  -  80  sec.  The  theoretical 
curves  are  superimposed  on  the  transverse  component, 
band-pass  amplitude  results  for  the  January  9,  1982, 
mainshock.  Fundamental  mode  is  the  solid  line;  suc¬ 
cessively  higher  order  modes  are  dashed  lines  of  de¬ 
creasing  dash-length. 

Theoretical  dispersion  predicted  by  the  Path  3  velocity 
model  for  fundamental  and  first  four  higher  mode  Ray¬ 
leigh  waves  for  the  period  range  .2  -  100  sec.  The 
theoretical  curves  are  superimposed  on  the  vertical 
component,  band-pass  amplitude  results  for  the  January 
9,  1982,  mainshock.  Fundamental  mode  is  the  solid  line; 
successively  higher  order  modes  are  dashed  lines  of  de¬ 
creasing  dash-length. 

Conjugate  fault  planes  proposed  for  the  New  Brunswick 
earthquakes  {from  Wetmiller  et.al.,  1983). 

Fundamental  mode  synthetic  Love-to-Rayleigh  spectral 
amplitude  ratios  for  depths  ranging  from  1  -  10  km 
for  the  Choy  mechanism:  6  *  195*,  &  »  65*, 

A  *  70*.  The  legend  for  depth  is  given  in  kilometers. 

Fundamental  mode  synthetic  Love-to-Rayleigh  spectral 
amplitude  ratios  for  depths  ranging  from  1  -  10  km 
for  the  Nabelek  mechanism:  8  «  3* { 183* ) ,  S  »  34*, 

A  ■  95*.  The  legend  for  depth  is  given  in  kilo¬ 
meters. 

Fundamental  mode  synthetic  Love-to-Rayleigh  spectral 
amplitude  ratios  for  depths  ranging  from  1  -  10  km 
for  the  Wetmiller  et.al.  mechanism  as  proposed  for 
the  January  9,  1982,  mainshock:  9  ■  195’,  i  ■  50*, 

A  *  120*.  The  legend  for  depth  is  given  in  kilo¬ 
meters. 

Fundamental  mode  synthetic  Love-to-Rayleigh  spectral 
amplitude  ratios  for  depths  ranging  from  1  -  10  km 
for  the  Wetmiller  et.al.  mechanism  as  proposed  for 
the  January  11,  1982,  aftershock:  8  »  332*,  S  “48*, 

A  ■  59*.  The  legend  for  depth  Is  given  in  kilo¬ 
meters. 

Observed  ratios,  calculated  from  the  band-pass 
results  for  the  January  9,  1982  mainshock,  super¬ 
imposed  on  the  theoretical  fundamental  mode  Love- 
to-Rayleigh  curves  (by  depth)  for  the  Choy  et.al. (1983) 
mechanism.  A  depth  range  of  4  -  8  km  is  Indicated  by 
these  points.  In  this  and  all  subsequent  figures,  stars 
Indicate  ratios  obtained  by  selecting  fundamental  mode 
arrivals  from  the  dispersion  curves,  finding  their  asso- 
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Figure  20 


Figure  21 


Figure  22 


Figure  23 


Figure  24 


Figure  25 


amplitudes  on  the  transverse  and  vertical  components 
and  calculating  their  ratios.  For  band-passes  with 
central  frequencies  in  the  plotted  period  range,  the 
arrival  corresponding  to  the  peak  amplitude  is  assured 
to  be  the  fundamental  mode,  so  peak  amplitude  ratios 
were  calculated  for  each  transverse-to-vertical  pair 
available.  This  did  not  contribute  much  new  informa¬ 
tion  for  the  mainshock,  but  helped  provide  additional 
data  points  for  the  aftershocks,  where  even  fundamen¬ 
tal  mode  identification  was  difficult. 

Observed  ratios,  calculated  from  the  band-pass  results 
for  the  January  11,  1982,  aftershock,  superimposed 
on  the  theoretical  Love-to-Rayleigh  curves  (by  depth) 
for  the  Wetmiller  et.al.  (1983)  January  11  mechanism. 
Although  constrained  by  a  single  value  near  the  peak, 
the  data  points  are  consistent  with  the  theoretical 
curve  for  the  6  km  depth. 

Observed  ratios,  calculated  from  the  band-pass  results 
for  the  January  9,  1982,  aftershock,  superimposed 
on  the  theoretical  Love-to-Rayleigh  curves  (by  depth) 
for  the  Choy  et.al.  (1983)  mechanism.  Although 
not  really  consistent  with  the  theoretical  curves 
for  this  model,  the  trend  to  peak  at  shorter 
period,  corresponding  to  shallower  depth  than  the 
mainshock.  Is  evident. 

Comparison  of  synthetic  Love-to-Rayleigh  spectral 
amplitude  ratios  for  modified  Choy  mechanisms  for 
depths  ranging  from  4  -  9  km.  In  all  cases, 

9  -  195*;  (a)  S  •  55*.  A  -  70’;  (b)  i  «  55* , 

A  •  90*;  (c)  i  •  55*,  A  -  110*;  (d)  i  «  65*, 

A  -  70*;  (e)  i  »  65*.  A  •  90*;  (f)  S  »  65*. 

A  -  110*. 

Band-pass  results  for  the  January  9,  1982,  aftershock 
superimposed  on  the  synthetic  Love-to-Rayleigh  ratios 
(bv  depth)  for  the  modified  Choy  mechanism  (  9  »  195*, 

*  ■  55’,  A  ■  90*)  for  depths  from  4  to  9  km.  Al¬ 
though  the  data  are  scattered  and  fall  beyond  the 
depth  ranges  covered  by  the  theoretical  curves,  a 
depth  less  than  4  km  is  Indicated. 

Band-pass  results  for  the  June  16,  1982,  aftershock 
superimposed  on  the  synthetic  Love-to-Rayleigh  ratios 
(by  depth)  for  the  Wetmiller  et.al. (1983)  January  11, 
1982,  mechanism.  The  values  suggest  a  range  of  possible 
focal  depths  of  7  -  9  km. 

GDSN  first-motion  data  for  the  January  9,  1982,  New 
Brunswick  mainshock,  recorded  at  GRFO;  for  the  January 
9,  1982,  aftershock,  recorded  at  BOCO;  and  for  the 
January  11,  1982,  aftershock  recorded  at  ZOBO. 
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Figure  26  Choy  et.al.  (1983)  preferred  focal  mechanism  (solid 

line),  8  »  195*,  i  *  65*,  A  ■  70*.  Triangles  are 
takeoff  angles  of  P  and  pP;  squares  are  sP  takeoff 
angles.  All  takeoff  angles  are  plotted  on  a  lower 
hemisphere  plot  (after  Choy  et.al.,  1983). 

Figure  27  Composite  distribution  of  January  and  April  aftershocks 

from  the  local  surveys  (after  Wetmiller  et.al.,  1983) 

Figure  28  Love-to-Rayleigh  spectral  amplitude  ratios  plotted  for 

a  variety  of  depth  ranges,  compared  with  data  for  the 
January  9  mainshock.  Spectral  contributions  for  differ¬ 
ent  depth  ranges  are  calculated  for  Love  and  Rayleigh 
amplitudes  and  ratios  plotted. 

(a)  Depth  ranges  include:  2,  3,  4  km;  5,  6,  7,  km; 

8,  9,  10  km. 

(b)  Depth  range  shown:  7,  8  km. 

Figure  29  Best-fitting  Love-to-Rayleigh  spectral  amplitude  ratios 

for  the  January  9  mainshock.  Spectral  contributions 
for  4,  5,  and  8  km  depths  are  calculated  for  Love  and 
Rayleigh  amplitudes  and  ratios  plotted. 
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TAYLOR/TOKSOZ  ASSUMED  PATH  G  STRUCTURE 
BANDPASS  RESULTS  FOR  1/9/82  MAINSHOCK 
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LOVE/RAYLEIGH  SPECTRAL  AMPLITUDE  RATIO 

NABELEK  MECHANISM 
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LOVE/RAYLEIGH  SPECTRAL  AMPLITUDE  RATIO 

NEW  BRUNSWICK  1/9/82  MAINSHOCK  -  BANDPASS  RESULTS 
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LOVE/RAYLEIGH  SPECTRAL  AMPLITUDE  RATIO 

NEW  BRUNSWICK  1/9/82  AFTERSHOCK  -  BANDPASS  RESULTS 
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Figure  29 


APPENDIX  1 


Harkrlder  Program 
Description 
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The  use  of  Harkrider's  layered  half-space  dispersion  pro¬ 
gram,...  has  been  accomplished  much  like  one  would  handle 
a  folk  tale  -  handed  down  from  one  graduate  student  to  the 
next. 

-  Turnbull  (1976) 

The  purpose  of  this  appendix  is  to  facilitate  future  use  of  the 
Harkrider  programs  by  describing  specifically  the  steps  to  be  followed 
in  preparing  the  data  and  control  parameters  for  input  to  the  main  pro¬ 
gram  which  generates  the  synthetic  seismograms.  The  routines  were 
written  in  single  precision  for  use  on  the  PRIME  computer  and  sent  to  us 
in  that  form.  The  computations  performed  in  this  study  were  carried  out 
on  an  IBM  computer,  so  the  programs  were  modified  as  needed  in  order  to 
accommodate  the  different  system  requirements.  The  original  programs 
and  support  libraries  needed  for  execution  on  the  PRIME  have  been  stored 
in  the  first  17  files  on  tape  HARK1  (standard  labeled,  6250  bpi)  which 
resides  in  the  PSU  IBM  tape  library  (a  listing  of  tape  contents  can  be 
obtained  quickly  by  running  PSUDEBE  on  the  tape).  The  modified  versions 
for  use  on  the  I8M  are  contained  in  the  last  files  on  HARK1  and  these 
are  the  versions  which  will  be  discussed  here.  It  must  be  pointed  out 
that  the  compiler  used  to  execute  the  main  program  was  an  H-extended, 
translated  version  of  VTPLOT,  which  modified  VS  Fortran  as  needed  to 
process  the  data  for  plotting  on  the  Versatec.  After  Fall  semester, 
1983,  this  compiler  was  no  longer  supported,  so  conversions  will  be 
necessary  to  ensure  full  compatibility  with  VS  Fortran  (if  continued 
IBM  use  Is  anticipated). 

The  first  step  in  preparing  the  data  is  to  scan  frequency/phase 
velocity  space  for  approximate  roots  to  the  period  equation.  Input 


98 


parameters  to  SCAN  are  defined  as  they  were  used  in  this  application; 
further  work  with  the  program  may  reveal  subtleties  in  their  meanings 
that  were  missed. 

LMODET  =  NCASE  =  the  number  of  modes  to  process 

IDUM  *  0 

INST  =  0  =  presumed  to  be  related  to  instrinent  correction,  set  equal 
to  zero  since  the  SCP  data  were  already  deconvolved  for  instru¬ 
ment  response 

IDENT  =  80  characters  available  for  commenting  the  trial  run 

N  ■  the  number  of  layers  in  the  structure  model 

NW  -  1  if  there  is  no  water  layer  included  in  the  structure 
(0  if  there  is) 

H,  ALP,  BET,  RHO  =  layer  parameters;  thickness  of  layer,  compres- 
sional  velocity,  shear  velocity,  and  density  should  be  indi¬ 
cated  for  each  of  the  N  layers 

CL,  CS,  NPTC  *  largest,  smallest  and  number  of  values  of  phase 
velocity  to  be  scanned  in  the  region  of  interest;  in  this 
case,  the  largest  was  just  greater  than  the  specified  shear 
velocity  in  the  half-space,  the  smallest  was  just  less  than 
the  shear  velocity  in  the  uppermost  layer;  100  was  felt  to 
be  a  reasonable  nunber  of  values  to  search  through  (130  might 
have  been  simpler  given  the  specified  phase  velocity  range) 

FL,  FS,  NPTT  =  largest,  smallest  and  nunber  of  points  to  be  scanned 
in  frequency;  our  range  was  chosen  as  1  -  .01  Hz  (1  sec  -  100 
sec  period)  on  the  basis  of  the  range  of  values  to  be  studied 
and  with  some  attention  to  time  cost 

NMODES  *  the  number  of  modes  requested;  for  want  of  a  better  idea, 
the  nunber  of  modes  was  held  the  same  (NMODES  ■  LMODET  =  NCASE). 
This  can  be  sensibly  varied  in  later  steps;  the  reasons  for 
retaining  three  different  control  parameters  at  this  stage  is 
unclear. 


The  program  is  set  to  scan  for  both  Rayleigh  and  Love  modes  automatically. 
If  two  sets  of  the  above  input  variables  are  specified.  The  original 
PRIME  version  (RLSCFRQ)  allows  for  one  or  the  other  depending  on  the 
values  specified  for  NLO  and  NRA. 
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Output  of  this  scanning  routine  includes  listings  of  input  and 
computed  layer  parameters,  a  visualization  of  the  sign  of  the  period  equa¬ 
tion  in  frequency/ phase  velocity  space  (the  roots  are  sought  where  the 
function  has  a  zero  value),  and  “INPUT  FILES  FOR  LATER  CALCULATIONS." 

The  output  of  interest  appears  under  this  last  heading  where  starting 
phase  velocities  and  frequencies  are  indicated  by  mode.  The  first  entry 
is  for  the  fundamental  mode;  subsequent  entries  correspond  to  as  many 
higher  order  modes  as  were  requested.  Since  there  were  100  frequency 
points  to  be  tested  in  the  case  illustrated  in  Figure  A1  -  1,  the 
frequencies  for  the  fundamental  mode  cover  the  full  range:  100  fre¬ 
quency  values  spanning  the  range  from  .01  to  1  Hz,  with  a  starting  phase 
velocity  value  of  4.1982  km/sec.  The  first  higher  mode  will  have  its 
first  root  at  a  frequency  value  of  .08  Hz,  with  93  roots  extending  out 
to  1  Hz,  and  the  initial  value  for  phase  velocity  should  be  set  at 
4.5429  km/sec.  These  are  the  starting  values  which  must  be  used  as 
input  to  the  Love  and  Rayleigh  eigenvalue-generating  programs,  which 
search  the  frequency/phase  velocity  space  on  a  finer  grid  for  more 
accurate  values  of  the  roots. 

Input  to  RAYLIT.UPEC  and  LOVIT.UP,  the  programs  for  more  precise 
determination  of  Rayleigh  and  Love  eigenvalues,  respectively,  takes  on 
the  same  form  for  both.  LPAR  is  currently  not  used  in  any  of  the  pro¬ 
grams  -  its  purpose  is  unknown. 

NCASE  «  has  been  assumed  to  be  the  number  of  modes  to  process 

NNN  »  setting  this  not  equal  to  zero  indicates  that  the  same 
structure  model  will  be  applied  for  all  NCASE 

MPRNT  «  controls  the  amount  of  diagnostic  output;  if  equal  to  0, 
maximum  output  is  listed;  if  equal  to  1,  only  the  layer  para¬ 
meters  and  dispersion  calculations  are  listed 


IDENT  *  again,  80  characters  of  identifying  comment 

N  and  NW  *  number  of  layers  and  water  layer  (0  or  1)  (as  specified 
in  SCAN) 


The  following  are  specified  first  for  the  fundamental  mode: 


Cl  =  starting  phase  velocity  (value  from  SCAN) 

DC  =  increment  for  phase  velocity  (.01  works  well) 

LC  =  number  of  frequencies  to  be  calculated  for  this  mode  (100  in 
the  case  of  the  fundamental  mode  for  this  example;  values 
for  higher  modes  are  the  nunber  of  roots  specified  by  SCAN) 

LD  *  controls  the  root-searching  process  and  is  described  within 
the  program 

NLAP  =  appears  to  control  input,  a  choice  of  zero  was  consistent 

NHS  =  0 

FK  *  starting  frequency  (specified  in  SCAN) 

FKD  =  frequency  decrement  =  -.01  Hz  in  this  case 

NLES  =  number  of  layers  to  ignore  in  the  calculation.  This  is  set 
at  zero  to  begin,  the  program  performs  layer  reduction  as 
needed  to  eliminate  instabilities  as  computation  proceeds 

LHF  *  0 

DKS  =  this  specifies  even/odd  modes;  the  value  should  be  set  at 
+1.0  for  fundamental  mode,  -1.0  for  first  higher  mode,  +1.0 
for  second  higher  mode,  etc. 

DELTA  AND  CPERT  =  parameters  which  control  the  perturbation  in 
phase  velocity  and  decide  when  the  root  has  been  achieved 
to  sufficient  accuracy;  they  are  set  within  the  program 
if  not  specified  as  input 


Layer  parameters  must  again  be  specified  as  they  were  for  the  scanning 
routine,  followed  by  starting  phase  velocity  and  frequency  values  for  as 
many  higher  order  modes  as  desired,  adhering  to  the  same  format  as  for 
the  fundamental  mode. 

The  output  from  RAYLIT.UPEC  and  LOVIT.UP  can  include  as  much  detail 
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1***  IBPOT  FILES  FOB  LATEB 

6 

CALCOLATIOMS 

1  1 

•** 

SCAB  FOB  1 ZB 

B80BSHICK  STBOCTOBE  -  PATHG 

-  BATLEIGR 

I  *  9 

MB  « 

1 
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0.010050 
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5. 00 CO 

6.0000 
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2.8000 
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3.  8000 
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10.0000 
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0.0400 
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0.050  0 
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0.080050 

-1.0 

0.00556 
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0.0300 
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Figure  A1  -  1 


as  desired,  controlled  by  the  value  specified  for  MPRNT.  For  example, 
to  ensure  that  the  correct  mode  has  been  found  in  the  root-searching 
process,  layer  displacements  can  be  checked  in  the  detailed  listing. 

The  fundamental  mode  should  show  no  zero  crossings,  while  each  suc¬ 
cessively  higher  order  mode  should  show  one  additional  zero  crossing. 

The  information  needed  for  the  main  program  (EQTR.ZBIG  for  the  PRIME; 
HARKRUNR  and  HARKRUNL  for  the  IBM),  are  the  values  (for  each  frequency 
specified)  of  the  phase  velocity,  group  velocity,  medium  amplitude  res¬ 
ponse,  ellipticity  (In  the  case  of  Rayleigh  waves)  and  an  integer  value 
Indicating  the  nunber  of  layers  used  In  the  calculation.  These  can 
either  be  stored  directly  In  a  file  or  modified  from  the  listed  output 
for  later  use  with  EQTR.ZBIG. 

Commenting  for  the  main  program  is  very  complete,  so  far  as  it  goes. 
The  crucial  input  for  the  case  of  a  double  couple  source  rests  with 


eight  lines  of  input  parameters: 


RE  AO 

( 

) 

NSR 

READ 

( 

) 

LMODO, NTYF, THETA, DELT, FLAM, LMODES 

READ 

( 

) 

LRDP , LEXP, SP, SSV, SSH 

READ 

( 

) 

OIST,TO,DT,LNT,LFFT,LINSY,LATN,LU 

READ 

( 

) 

NO  I SL , STO , DTI , DT2 , DT3 

READ 

( 

) 

N SCALE, XSCALE, SCAF , LPLT , I  BUG , LSAVE , LMOPLT 

READ 

( 

) 

(DEPT(I),  1=  1,  NDEP) 

These  quantities  are  defined  as  follows: 

NSR  =  1  if  only  one  source  structure  is  to  be  considered 

LMODO  =  0  for  the  Love  wave  run;  =  1  for  the  Rayleigh  wave  run 

NTYF  =  2,  specifying  the  douple  couple 

THETA  ■  station  azimuth  measured  counter-clockwise  **  from  the  strike 
of  the  hanging  wall.  This  definition  must  be  considered  care¬ 
fully,  as  it  appears  to  be  the  reverse  of  current  conven¬ 
tions.  Theta  and  other  parameters  of  the  source  geometry  are 
Illustrated  in  Figure  A1  -  2. 

DEIT,  FLAM  *  fault  dip  and  rake  as  defined  in  Figure  A1  -  2 

LMODES  a  the  number  of  modes  to  include  in  the  synthetic  calcula¬ 
tion  from  the  results  of  LQVIT.UP  or  RAYLIT.UPEC 

LRDP  *  4,  indicating  a  double  couple,  in  this  case,  with  no  explo¬ 
sion  component;  therefore,  LEXP  »  0,  and  SP,  SSV  and  SSH 
automatically  become  0 

01  ST  »  the  distance  frm  source  to  receiver  which  may  be  specified 
in  kilometers  or  degrees 

T0,DT,LNT  *  starting  time,  time  interval,  and  number  of  time  points, 
respectively,  with  options  given  as  to  how  they  may  be  speci¬ 
fied 

LFFT  «  1,  calls  for  straight  line  Interpolation  FFT 


**  After  subsequent  work  by  later  Investigators  revealed  inconsistent 

results,  it  was  discovered  that  this  convention  was  reversed  in  final 
revisions  of  the  program.  Therefore,  THETA  should  be  measured  clock¬ 
wise  from  strike. 
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LINS  3  0,  indicates  that  instrument  response  will  not  be 
included  in  the  calculation 

LATN  =  0,  indicates  that  attenuation  values  will  not  be  incor¬ 
porated 


(in  the  present  study,  it  made  sense  to  set  LINS  and  LATN  equal  to  0 
since  amplitude  spectra,  rather  than  synthetic  seismograms,  were  calcu¬ 
lated  for  comparison  with  spectral  values  which  had  already  been  cor¬ 
rected  for  instrument  response  and  for  which  attenuation  was  assumed 
constant.  If,  however,  a  synthetic  seismogram  were  to  be  calculated 
for  comparison  with  the  original  SCP  record,  these  quantities  would  have 
to  be  specified. ) 


LU  3  0,  specifies  that  displacement  is  to  be  calculated  at  the 
surface 

NDISL  =  4,  and  STO  3  0  calling  for  a  step  moment  or  dislocation 
with  both  moments  set  to  10^5  dyne-cm 

DTI,  0T2,  DT3  3  all  equal  to  zero  because  of  the  previous  choice 
of  step  moment 

NSCALE,  XSCALE,  SCAF  =  are  all  scaling  calls  for  the  seismogram  plot 
They  may  be  specified  or  default  values  are  set. 

LPLT  3  0,  calls  for  no  plotting;  if  >  0,  the  numerical  value  gives 
the  number  of  plots  per  page 

IBUG  3  controls  the  amount  of  output  to  be  listed 

LSAVE  3  allows  the  time  series  to  be  saved  in  an  external  file 

LMOPLT  «  1,  specifies  that  each  individual  modal  contribution  is  to 
be  plotted,  as  well  as  the  accumulated  time  series.  If  equal  to 
0,  only  the  final,  accumulated  time  series  is  plotted. 

NDEP  3  the  number  of  source  depths  for  which  complete  calculations 
are  to  be  performed 

LCHG  3  0 

DEPT(I),  I  3  l.NOEP  3  the  values  of  depth  to  use  in  the  calculations 
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Once  all  these  source  parameters  have  been  processed,  the  structure 
parameters  are  called  within  SUBROUTINE  LAYERS  and  calculation  proceeds. 
It  is  a  minor  modification  to  obtain  spectral  amplitudes  rather  than 
the  time  series;  just  prior  to  the  call  to  COOLB,  the  values  can  be 
retrieved. 

With  these  hints  about  the  flow  of  the  program  steps,  accompanied 
by  sample  output  listings  for  immediate  reference,  it  should  be  a  rou¬ 
tine  matter  to  get  the  Harkrider  routines  up  and  running  and  performing 
over  the  wide  range  of  applications  for  which  they  are  designed. 


*1 


Figure  A1  -  2 

Geometry  of  source  parameters  and  relative  position  of 
receiving  station  on  the  free  surface  (after  Ben-Menahem 
and  Harkrider,  1964) 


APPENDIX  2 


Original  Digital  Seismograms 
at  SCP 


Figures  A2  -  1  through  A2  -  12  are  the  original  digital  seismograms 
as  recorded  at  SCP.  Short,  intermediate,  and  long  period  data  are  dis¬ 
played  for  each  event  along  with  information  on  origin  time,  epicentral 
location,  magnitude,  depth  and  record  start  time.  Short  period  data  are 
sampled  at  20  samples  per  second;  intermediate  period  at  10  samples  per 
second;  long  period  at  1  sample  per  second.  The  signal-to-noise  ratio 
is  good  for  both  the  January  9  mainshock  and  the  January  11  aftershock 
(Figures  A2-1  through  A2-3  and  A2-7  through  A2-9).  The  long  period  sig¬ 
nal  for  both  the  January  9  aftershock  and  the  June  16  event  (Figures 
A2-6  and  A2-12)  are  barely  discernible  from  the  noise.  Note  that  the 
intermediate  period  instrument  did  not  begin  recording  data  for  the 
June  16  event  (Figure  A2  -  11)  until  well  after  the  first  arrival.  This 
posed  significant  difficulty  when  the  band-pass  filters  were  applied  to 
this  data,  resulting  in  very  few  usable  group  velocity  data  points. 


Figure  A2 
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Figure  A2 
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Band-Pass  Filtered  Seismograms 
From  Ma inshock 
And  Three  Aftershocks 
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The  figures  in  Appendix  3  include  all  bandpass  filter  results,  as 
obtained  from  the  short,  intermediate,  and  long  period  data  for  the  four 
events  considered  in  this  paper.  Each  figure  includes  the  original  seis¬ 
mogram  (at  the  bottom  of  the  page)  and  results  of  applying  a  suite  of 
band-pass  filters  centered  on  selected  frequencies.  The  time  scale  is 
25  seconds  per  tick  mark  at  the  top  of  the  figure.  On  the  left,  the 
filter  parameters  (lower  cutoff,  lower  corner,  upper  corner,  and  upper 
cutoff  frequencies  as  described  in  Chapter  4)  are  given  for  each  band¬ 
pass  and,  on  the  right,  are  peak  amplitudes,  by  trace,  expressed  in 
microns.  Each  trace  is  normalized  to  this  peak  amplitude.  Filter 
parameters  common  to  each  suite  are  detailed  on  a  separated  page  for 
clarity.  In  some  cases,  only  the  envelopes  of  the  band-passes  are 
displayed  in  an  effort  to  bypass  the  significant  cost  of  plotting 
the  carrier  signal . 

Figures  A3-1  through  A3-12  include  three  suites  of  short  period 
filters  ordered  by  event  date  (January  9  mainshock,  January  9  aftershock, 
January  11  aftershock,  and  June  16  event).  Figures  A3-13  through  A3-42 
are  the  results  of  applying  four  suites  of  intermediate  period  filters, 
again  ordered  by  event  date.  Figures  A3-43  through  A3-54  present  the 
long  period  results  for  all  but  the  June  16  event  where  noise  obscured 
the  filtered  signal  completely.  Short  period  filters  were  applied  to 
the  vertical  component  (SPZ)  only;  intermediate  and  long  period  filters 
were  applied  to  both  vertical  and  transverse  components  (IPZ,  IPT,  LPZ, 
LPT). 
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The  filter  parameters  for  Figures  A3  -  1  through  4  are  as 


indicated: 

FL  =  8.470 

FL  =  7.S7G 

FL  =  7.470 

FL  =  6.970 

FL  =  6.470 

FL  =  5.970 

FL  =  5.470 

FL  =  4.970 


FLC  =  8.475 

FLC  =  7.975 

FLC  =  7.475 

FLC  =  6.975 

FLC  =  6.475 

FLC  =  5.975 

FLC  =  5.475 

FLC  »  4.975 


FHC  =  8.525 

FHC  =  8.025 

FHC  =  7.525 

FHC  =  7.025 

FHC  =  6.525 

FHC  =  6.025 

FHC  *  5.525 

FHC  =  5.025 


FH  =  8.530 

FH  =  8.C30 

FH  =  7.530 

FH  =  7.030 

FH  =  6.530 

FH  =  6.030 

FH  =  5.530 

FH  =  5.030 


OBSERVED  SIGNAL 


JA^ll  lA^Il  2An>>|l  2Am 


Ht*  8«UNS«ICA 


J/9/82 


:  :  S :  36 :  -«2 . 3 


Figure  A3  -  2 
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Figure  A3  -  4 


The  filter  parameters  for  Figures  A3  -  5  through  8  are  as 
indicated: 
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FL  =  4.470  FLC  =  4.475 

FL  =  3.970  FLC  =  3.975 

FL  =  3.470  FLC  =  3.475 

FL  =  2.970  FLC  =  2.975 

FL  =  2.470  FLC  =  2.475 

FL  =  1.970  FLC  =  1.975 

FL  =  1.470  FLC  =  1.475 

FL  =  0.970  FLC  =  0.975 


FHC  =  4.525  FH  =  4.530 

FHC  =  4.025  FH  =  4.030 

FHC  =  3.525  FH  =  3.530 

FHC  =  3.025  FH  =  3.030 

FHC  =  2.525  FH  =  2.530 

FHC  =  2.025  FH  *  2.030 

FHC  =  1.525  FH  =  1.530 

FHC  =  1.025  FH  =  1.030 


OBSERVED  SIGNAL 


Figure  A3  -  5 
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Figure  A3  -  8 
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The  filter  parameters  for  Figures  A3  -  9  through  12  are  as 
indicated: 


FL  =  0.770 

FLC  =  0.775 

FHC  =  0.825 

FH  =  0.830 

FL  =  0.670 

FLC  =  0.675 

FHC  =  0.725 

FH  =  0.730 

FL  =  0.570 

FLC  =  0.575 

FHC  =  0.625 

FH  =  0.630 

FL  =  0.470 

FLC  =  0.475 

FHC  =  0.525 

FH  =  0.530 

FL  =  0.370 

FLC  =  0.375 

FHC  =  0.425 

FH  =  0.430 

FL  =  0.270 

FLC  =  0.375 

FHC  =  0.325 

FH  =  0.330 

FL  =  0.170 

FLC  =  0.175 

FHC  =  0.225 

FH  =  0.230 

FL  =  0.070 

FLC  =  0.075 

FHC  =  0.125 

FH  =  0.130 

OBSERVED  SIGNAL 
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Figure  A3  -  9 
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Figure  A3  -  10 
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The  filter 
indicated: 

FL  =  4.470 

FL  =  3.970 

FL  =  3.470 

FL  =  2.970 

FL  =  2.470 

FL  =  1.970 

FL  =  1.470 

FL  =  0.970 


parameters  for  Figures  A3 


13  through  20  are  as 


FLC  =  4.475 

FLC  =  3.975 

FLC  =  3.475 

FLC  =  2.975 

FLC  «  2.475 

FLC  =  1.975 

FLC  =  1.475 

FLC  =  0.975 


FHC  =  4.525 

FHC  =  4.025 

FHC  =  3.525 

FHC  =  3.025 

FHC  =  2.525 

FHC  =  2.025 

FHC  =  1.525 

FHC  =  1.025 


FH  =  4.530 

FH  =  4.030 

FH  =  3.530 

FH  =  3.030 

FH  =  2.530 

FH  =  2.030 

FH  =  1.530 

FH  =  1.030 


OBSERVED  SIGNAL 
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Figure  A3  -  17 
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Figure  A3  -  19 
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The  filter  parameters  for  Figures  A3  -  21  through  28  are  as 


indicated: 

FL  =  0.970  FLC  =  0.975 

FL  =  0.870  FLC  =  0.875 

FL  =  0.770  FLC  =  0.775 

FL  =  0.670  FLC  =  0.675 

FL  =  0.570  FLC  =  0.575 

FL  *  0.470  FLC  =  0.475 

FL  =  0.370  FLC  =  0.375 

FL  =  0.270  FLC  =  0.275 


FHC  =  1.025  FH  =  1.030 

FHC  =  0.925  FH  =  0.930 

FHC  =  0.825  FH  =  0.830 

FHC  =  0.725  FH  -  0.730 

FHC  *  0.625  FH  =  0.630 

FHC  =  0.525  FH  =  0.530 

FHC  =  0.425  FH  =  0.430 

FHC  =  0.325  FH  =  0.330 


OBSERVED  SIGNAL 


Figure  A3  -  21 
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The  filter  parameters  for  Figures  A3  -  29  through  36  are  as 
indicated: 


FL  =  0.220 

FL  =  0.170 

FL  =  0.120 

FL  =  0.090 

FL  *  0.092 

FL  =  0.082 

FL  =  0.063 

FL  =  0.042 


FLC  =  0.225 

FLC  =  0.175 

FLC  =  0.125 

FLC  =  0.095 

FLC  =  0.097 

FLC  =  0.087 

FLC  =  0.067 

FLC  =  0.047 


FHC  =  0.275 

FHC  =  0.225 

FHC  =  0.175 

FHC  =  0.145 

FHC  =  0.102 

FHC  =  0.092 

FHC  =  0.072 

FHC  0.052 


FH  =  0.280 

FH  =  0.230 

FH  =  0.180 

FH  =  0.150 

FH  =  0.107 

FH  =  0.097 

FH  =  0.077 

FH  =  0.057 


OBSERVED  SIGNAL 


4 


I 


New  8«UNS«{C'»  -  :/$/ 82  :i::3:Sl.S 


STaTlCN:  SC*  1NT£«W  *£RICO  -  2  :*a«T  t;w£:  12:56:57.^8 


aar«  u  aitwia  »>  t 


350 


Figure  A3  -  30 


nEw  B"UN5*IC*  -  1/9/32  3-:  :6:36;«t2.3 

5**M3*««  SCP  INTERN  RE«I30  -  l  l’a*r  r:*E:  I S :  wO:  48.  -.8 

it  HC I'  : 


Figure  A3  -  33 


163 


The  filter  parameters  for  Figures  A3  -  37  through  42  are  as 


indicated: 

FL  «  0.009 

FL  =  0.011 

FL  =  0.012 

FL  =  0.015 

FL  =  0.017 

FL  *  0.021 

FL  *  0.026 

FL  =  0.032 


FLC  *  0.014 

FLC  =  0.016 

FLC  =  0.017 

FLC  =  0.020 

FLC  =  0.022 

FLC  =  0.026 

FLC  =  0.031 

FLC  =  0.037 


FHC  =  0.019 

FHC  =  0.021 

FHC  =  0.022 

FHC  =  0.025 

FHC  ■  0.027 

FHC  =  0.031 

FHC  ■  0.036 

FHC  =  0.042 


FH  =  0.024 

FH  =  0.026 

FH  =  0.027 

FH  =  0.030 

FH  =  0.032 

FH  =  0.036 

FH  =  0.041 

FH  =  0.047 


OBSERVED  SIGNAL 
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Figure  A3  -  42 
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The  filter  parameters  for  Figures  A3  -  43  through  48  are  as 


indicated: 

FL  =  0.082  FLC  =  0.087 

FL  =  0.063  FLC  =  0.067 

FL  =  0.042  FLC  =  0.047 

FL  =  0.032  FLC  =  0.037 

FL  =  0.022  FLC  »  0.027 

FL  ■  0.012  FLC  ■  0.017 

FL  =  0.007  FLC  =  0.012 

FL  =  0.002  FLC  =  0.007 


FHC  =  0.092  FH  =  0.097 

FHC  =  0.072  FH  =  0.077 

FHC  =  0.052  FH  =  0.057 

FHC  =  0.042  FH  =  0.047 

FHC  =  0.032  FH  =  0.037 

FHC  =  0.022  FH  =  0.027 

FHC  =  0.017  FH  =  0.022 

FHC  =  0.012  FH  =  0.017 


OBSERVED  SIGNAL 


Figure  A3  - 
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The  filter  parameters  for  Figures 
indicated: 

FL  =  0.003  F|_C  =  0.008 

FL  =  0.004  FIX  =  0.009 

FL  =  0.004  FLC  =  0.009 

FL  =  0.005  FLC  =  0.010 

FL  *  0.005  FLC  =  0.010 

FL  =  0.009  FLC  =  0.014 

FL  =  0.017  FLC  =  0.022 

FL  ■  0.092  FLC  =  0.097 


A3  -  49  through  54  are  as 

FHC  -  0.013  FH  =  0.018 

FHC  =  0.014  FH  =  0.019 

FHC  =  0.014  FH  =  0.019 

FHC  =  0.015  FH  =  0.020 

FHC  =  0.015  FH  =  0.020 

FHC  =  0.019  FH  =  0.024 

FHC  =  0.027  FH  =  0.032 

FHC  =  0.102  FH  =  0.107 
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Figure  A3  -  54 
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Appendix  4  presents  a  compilation  of  group  velocity  results  for  all 
events  superimposed  and  each  event  separately,  with  overlying  theoretical 
dispersion  curves  appropriate  for  the  Love  and  Rayleigh  modes.  As  in 
Figures  10-13,  the  different  symbols  correspond  to  values  of  peak  ampli¬ 
tudes.  Each  peak  amplitude  is  normalized  to  the  maximum  amplitude  for 
the  given  band-pass  and  then  assigned  a  synbol  according  to  the  following 
scheme: 

10  -  the  point  having  the  peak  amplitude  for  the  given  central 
frequency 

9  -  those  point(s)  representing  peaks  along  the  filtered  signal 
with  ampl itude( s)  80-99%  of  maximin  for  that  trace 

8  -  those  point(s)  representing  peaks  along  the  filtered  signal 
with  amplitude(s)  60-79%  of  maximum  for  that  trace 

7  -  those  point(s)  representing  peaks  along  the  filtered  signal 
with  amplitudes  less  than  59%  of  maximum  for  that  trace 

This  scale  made  possible  rapid  identification  of  significant  arrivals  by 
visual  inspection. 

Figure  A4-1  shows  the  results  for  the  vertical  component  band-pass 
filters  for  all  four  events  with  theoretical  Rayleigh  dispersion  super¬ 
imposed;  Figure  A4-2  shows  this  same  information  separately  for  each 
event.  Similarly,.  Figure  A4-3  shows  the  Love  mode  results  for  all 
events,  while  Figure  A4-4  shows  the  results  separately  by  event. 
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DISPERSION  RESULTS  -  COMPARISON 
TAYLOR/TOKSOZ  ASSUMED  PATH  G  STRUCTURE 
BANDPASS  RESULTS  FOR  ALL  EVENTS 
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DISPERSION  RESULTS  -  COMPARISON 
TAYLOR/TOKSOZ  ASSUMED  PATH  G  STRUCTURE 
BANDPASS  RESULTS  FOR  ALL  EVENTS 
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Appendix  5  presents  theoretical  Love  and  Rayleigh  spectra  and 
Love-to-Rayleigh  spectral  ratios,  by  depth,  for  the  four  source 
mechanisms  studied  in  this  paper.  Fundamental  mode  values  are 
plotted  as  solid  lines.  Successively  higher  order  modes  are  repre¬ 
sented  by  lines  of  increasing  dash-length.  Depths  of  one  to  ten 
kilometers  are  considered  for  each  model. 

In  each  figure,  (a)  shows  the  theoretical  Love  amplitude  spectrum; 
(b)  shows  the  theoretical  Rayleigh  amplitude  spectrum;  and  (c)  shows 
the  Love-to-Rayleigh  spectral  ratio. 

Figures  A5  -  1  through  5  are  results  for  source  parameters  of  the 
Choy  mechanism:  9  =  318°,  S  =  65°,  A=  70°. 

Figures  A5  -  6  through  10  are  results  for  source  parameters  of  the 
Nabelek  mechanism:  9  *  306°,  6  =  34°,  A  =  95°. 

Figures  A5  -  11  through  15  are  results  for  source  parameters  of 
the  Wetmiller  mechanism  for  the  mainshock:  8  =  318°,  i  =  50°, 

A  =  120°. 

Figures  A5  -  16  through  20  are  results  for  source  parameters  of 
the  Wetmiller  mechanism  for  the  January  11  aftershock:  •  =  95°, 

S  =  480,  A  =  590. 
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Three  dimensional  renditions  of  the  peak  envelope  amplitudes  asso¬ 
ciated  with  the  band-pass  filtered  seismograms  were  created  and  are  dis¬ 
played  in  Appendix  6  using  a  Surface  II  plotting  routine  which  generates 
a  complete  grid  of  values  from  irregularly  spaced  points  by  applying  a 
simple  interpolation  scheme.  While  a  more  complete  three-dimensional 
view  would  result  from  similar  processing  of  the  entire  envelope  pro¬ 
duced  for  each  filtered  signal,  the  amount  of  information  anticipated 
from  this  effort  vs.  the  expense  of  re-running  the  band-passes  did  not 
justify  the  undertaking.  Peak  amplitudes  plotted  against  group  velocity 
vs.  period  for  Love  (b)  and  Rayleigh  (a)  waves  for  each  event  are  viewed 
both  from  a  clockwise  45’  and  a  135*  rotation  of  the  group  velocity 
vs.  period  plane,  so  that  behavior  at  higher  group  velocities  can  be 
seen,  especially  when  peaks  at  other  group  velocities  obscure  the  view. 
Discontinuities  in  the  lines  on  the  plot  are  artifacts  of  instabilities 
in  the  interpolation  scheme  when  amplitudes  take  on  extreme  values  in 
close  proximity  to  one  another.  Each  three-dimensional  plot  is  normal¬ 
ized  to  the  highest  amplitude  value  to  be  plotted.  The  long,  more 
coarsely  gridded  axis  represents  period,  with  values  ranging  from  1-10 
seconds;  the  short  axis  represents  group  velocity,  spanning  the  window 
from  2.5  to  4.5  km/sec.  Note  that  both  axes  are  linear.  Although  quan¬ 
titative  Interpretation  of  these  plots  is  not  feasible  for  obvious 
reasons,  there  is  one  very  striking  feature  in  the  Love  wave  results. 
With  the  exception  of  the  June  16  aftershock,  whose  features  are 
obscured  by  the  high  amplitude  arrival  near  ten  seconds  period,  the 
three  events  of  the  January  sequence  show  remarkable  similarities  in 
the  pattern  of  Love  wave  excitation.  A  slightly  lower  amplitude  in 
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the  vicinity  of  one  second  and  a  more  pronounced  peak  at  about  four 
seconds  in  the  plot  for  the  January  11  aftershock  are  the  only  signifi¬ 
cant  differences  noted( compare  Figures  A6-lb  through  4b).  This  is  signi¬ 
ficant  considering  that,  for  all  the  mechanisms  tested  in  this  study,  the 
theoretical  Love  wave  spectra  (Appendix  5)  were  nearly  invariant  -  only 
relative  amplitudes  differed.  The  obvious  inference  is  that,  for  high 
frequency  arrivals.  Love  wave  relative  modal  excitation  at  the  source 
is  nearly  the  same  from  event  to  event  and  that  this  similarity  is  being 
preserved  consistently  for  these  events  after  propagating  the  10.2* 
distance  to  SCP. 

The  Rayleigh  plots  (Figures  A6-la  through  4a)  are  considerably  more 
complicated,  consistent  with  the  character  of  the  theoretical  spectra  for 
Rayleigh  waves.  The  peak  amplitude  that  so  distinctly  shows  up  at  10 
seconds  period  is  most  assuredly  associated  with  the  fundamental  mode, 
but  by  virtue  of  the  controlling  effect  that  a  low  velocity  surficial 
layer  will  have  on  the  higher  frequency  fundamental  mode,  it  is  not 
expected  to  be  present  in  the  shorter  period  signals  for  group  velocities 
greater  than  3  km/sec.  It  is  in  this  shorter  period  higher  group  velocity 
range  that  it  becomes  possible  to  observe  the  higher  mode  behavior.  Com¬ 
parison  of  the  theoretical  spectra  of  the  chosen  mechanisms  and  preferred 
depths  for  the  various  events  shows  the  differences  that  can  be  expected. 
Overlaying  Figure  A5-2b  (for  3  or  4  km  depth)  and  Figure  A5-4b  (for  7 
km  depth)  shows  that  theoretically  this  depth  difference  accounts  for 
the  lower  amplitude  arrivals  at  10  seconds  period  for  the  January  9  after¬ 
shock  compared  with  higher  amplitude  contributions  from  both  fundamental 
and  first  higher  order  mode  at  10  seconds  for  the  mainshock  (compare 
Figures  A6-la  and  A6-2a).  Figures  A5-18b  at  6  km  for  the  January  11 
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aftershock  and  Figure  A5-4b  for  the  mainshock  suggest  that  the  amplitude 
at  10  seconds  for  the  fundamental  relative  to  the  higher  mode  arrivals  is 
smaller  for  the  aftershock  than  for  the  mainshock.  In  fact,  the  higher 
order  modes  have  comparable  amplitudes  to  the  fundamental  mode  amplitude 
for  the  January  11  aftershock. 

In  general  careful  comparison  of  the  three-dimensional  perspective 
plots  in  Figures  A1-A4  and  the  corresponding  band-pass  filters  for  differ¬ 
ent  events  at  common  central  frequencies  (Appendix  3),  shows  that  there 
is  significant  variation  in  relative  excitation  among  the  modes  comprising 
Lg  on  the  original  vertical  component  seismograms.  This  behavior  is  con¬ 
trast  to  that  for  the  Love  modes  discussed  earlier  and  is  consistent  with 
the  variations  in  theoretical  excitation  of  Rayleigh  modes  for  different 
depths  and  focal  mechanisms  as  evidenced  by  the  theoretical  results 
shown  in  Appendix  5.  However,  as  discussed  in  the  text,  these  varia¬ 


tions 


January  9,  1982 
Mainshock 


Normalized  to  Peak  Amplitude 


APPENDIX  7 


Errata 


221 


Since  December,  1983,  when  this  study  was  concluded  and  defended, 
subsequent  investigation  into  the  New  Brunswick  Earthquakes  has  attempted 
to  further  refine  source  parameters  using  the  same  program  but  including  a 
systematic  search  through  incremented  values  of  source  parameters  for  the 
best  match  of  theoretical  and  empirical  Love/Rayleigh  spectral  ratios. 
Problems  arising  from  this  ongoing  study  led  to  pointed  questioning  of  the 
published  conventions  of  source  geometry  utilized  by  the  program.  What  had 
been  thought  of  as  a  reverse  convention  by  most  who  studied  this  paper  was 
found  to  be  just  that;  it  is  to  Li 1 i  Tang's  credit  that  this  suspicion 
was  confirmed  with  Dr.  Harkrider.  The  published  convention  indicated 
counterclockwise  measurement  of  theta  (as  illustrated  in  Appendix  1);  for 
subsequent  programming  reasons,  clockwise  measurement  of  theta  from  strike 
was  preferred. 

Subsequent  work  at  The  Pennsylvania  State  University  (Yan,  1985)  has 
yielded  a  best  fit  to  the  mainshock  with  source  parameters  as  follows: 
theta  equals  42  degrees,  dip  equals  43  degrees,  and  rake  equals  130  degrees, 
most  like  the  results  proposed  by  Wetmiller  et.al.,  1983.  The  results  of 
the  present  paper  preferred  the  mechanism  proposed  by  Choy  (see  Table  5, 
page  38),  with  strike  of  195  degrees  (theta  =  318  degrees,  or  360  -  318 
degrees  by  the  revised  convention),  dip  of  65  degrees,  and  rake  of  70 
degree$(close  to  180  minus  the  130  degree  rake  preferred  in  Yan's  work). 
Considering  the  equation  for  the  radiation  pattern  and  the  dependence  of 
the  various  terns  and  factors  on  theta  and  lambda,  it's  Interesting  to 
note  that,  except  for  the  difference  in  dip  values,  the  synthetic 
Love/Rayleigh  amplitude  ratios  could  be  expected  to  be  nearly  the  same  for 
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the  two  cases.  With  theta  now  360  degrees  minus  its  former  value  and 
lambda  now  roughly  180  degrees  minus  its  former  value,  the  compensating 
effects  alter  the  signs  of  all  the  real  or  imaginary  terms,  keeping 
amplitudes  the  same.  Beyond  that  fortuitous  relationship  in  the  case  of 
the  mainshock,  however,  reasonable  comparison  ceases  for  the  other  events. 

Acknowledgement  of  this  error  in  applying  the  program  necessitates 
retraction  of  the  analyses  involving  comparison  of  synthetic  spectra  and 
empirical  results  obtained  from  the  band-pass  filters.  With  respect  to  the 
earthquakes  in  question,  these  interpretations  have  already  been  revised 
and  improved  through  the  subsequent  work  cited  (Yan,  1985).  Synthetic 
spectral  ratios  are,  thus,  submitted  within  the  body  of  this  paper  and 
Appendix  5  as  illustrative  of  the  diagnostic  nature  of  depth  when  using 
synthetic  surface  wave  spectra  in  comparison  with  empirical  Love/Rayleigh 
spectral  ratios.  They  are,  however,  no  longer  indicative  of  the  source 
parameters  for  the  events  studied. 


